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ABSTRACT : 

An  objective-numerical  scheme  for  locating  and  specifying  the  inten- 
sity of  hyperbaroclinic  zones  (HBZs)  is  presented  and  applied  in  both 
research  and  operational  environments.   The  operational  model,  called 
DG9*,  employs  the  isobaric  gradient  of  virtual  potential  temperature  (8*) 
to  specify  baroc linicity  while  a  directional  second  derivative  of  9*  is 
the  parameter  used  to  locate  the  HBZs  and  indirectly  give  a  measure  of 
their  intensity. 

The  DG9*  model  is  diagnosed  from  an  analysis  of  modeled  temperature 
profiles  and  a  number  of  real-atmosphere  temperature  fields  at  1000,  850, 
700,  and  500  mb ,  as  taken  from  winter  and  summer  situations  in  the  years 
1964-67.   Consideration  is  given  to  location,  intensity,  frequency,  slope, 
and  vertical  and  horizontal  distribution  and  consistency  of  the  HBZs  in 
relation  to  the  numerical  processing  of  data.   A  key  synoptic  time,  0000 
GMT  6  March  1966,  serves  to  exemplify  analytical  and  diagnostic  features 
of  the  HBZs.   The  numerically-analyzed  1000  mb  fronts, produced  operationally 
by  Fleet  Numerical  Weather  Central,  Monterey,  California  since  1965,  are 
statistically  and  graphically  compared   to  subjectively-analyzed  versions 
for  winter  and  summer  periods  up  to  30  days.   The  relation  of  the  HBZ  to 
vertical  motion  and  cloud  is  illustrated  by  a  case  study  from  December  1966. 

Experiments,  using  a  high-resolution  scheme  for  analyzing  the  HBZs,  are 
also  presented.   A  proposed  model,  employing  parameters  of  vertical  sta- 
bility and  isentropic  slope,  utilizes  the  thermal  detail  inherent  in 
vertically-resolved  radiosonde  data  to  give  an  enhanced  field  of  baro- 
clinicity  on  the  constant  pressure  surface.   In  addition,  it  is  shown  how 
the  model  may  be  coupled  with  a  vertical  extrapolation  procedure  to  at 
least  triple  the  amount  of  baroclinicity  data  used  for  HBZ  analysis. 
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HYPERBAROCLINIC   ZONES:      NUMERICAL  ANALYSIS 
AND  APPLICATION  TO  FRONTS* 

R.    J.    Renard 


CHAPTER  I 
HISTORICAL  BACKGROUND  AND  OBJECTIVES  OF  STUDY 

It  is  difficult  to  ascertain  exactly  when  meteorologists  first 
became  aware  of,  or  concerned  with,  atmospheric  baroclinicity .   According 
to  Bergeron  (1959),  Kbppen  was  one  of  the  first  nineteenth  century 
meteorologists  to  perceive  observational  evidence  of  air  masses  and  their 
baroclinic  boundaries,  conventionally  called  fronts .   However,  it  was  not 
until  the  1920s  and  1930s  that  the  so-called  Norwegian  air-mass  and  frontal 
concepts  (Bjerknes,  1918)  were  firmly  established.   These  ideas  live  on 
today  in  the  form  of  frontal  analyses  (e.g.,  warm,  cold  and  occluded 
fronts)  found  on  operational  weather  charts  of  every  major  meteorological 
service  in  both  hemispheres,  without  exception. 

Since  the  1940s,  the  broader  concept  of  a  baroclinic  frontal  zone , 

vice  a  frontal  line,  has  become  an  increasingly  important  theoretical  and 

synoptic  consideration  in  the  life  history  of  many  important  facets  of 

weather  systems.   In  this  connection,  the  term  hyperbaroclinic  zone  (HBZ) , 

after  Godson  (1951),  is  most  appropriate  to  describe  a  restricted  area  of 

relatively  large  values  of  baroclinicity,  without  necessarily  connoting  a 

zone  separating  identifiable  air  masses.   Many  investigations  in  the  pre- 

computer  era  have  demonstrated  the  importance  of  HBZs  to  the  dynamic  and 

thermodynamic  processes  and  structure  of  the  atmosphere  (Bjerknes,  1951; 

McClain  and  Danielsen,  1955;  Newton,  1958;  Palmen,  1951;  Palmen  and  Newton, 

1 
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1951;  Reed  and  Sanders,  1953;  Sanders,  1955).   A  synthesis  of  theories 

and  examples  of  development  and  motion  of  waves,  vortices,  jet  streams,  and 
cloud  and  precipitation  regimes  as  a  function  of  baroc linicity  may  be  found 
in  dynamics  texts  (e.g.,  Haltiner  and  Martin  [1957]  or  Hess  [1959])  and 
applied  manuals  (e.g.,  Riehl  et  al.  [1952,  1954],  the  United  States  Navy 
[1952]  and  the  United  States  Weather  Bureau  [I960]).   However,  quanti- 
tative knowledge  of  the  behavioral  characteristics  of  HBZs  on  a  synoptic 
hemispheric  scale  has  remained  incomplete.   This  is  due  to  the  short- 
comings imposed  by  subjective  analysis  with  its  attendant  ambiguity  in 
locating  and  depicting  the  important  HBZs  and  to  the  derived,  not  observed, 
nature  of  baroc linicity .   The  operational  convention  of  specifying  only 
the  warm-side  boundary  of  the  HBZ  (Taljaard  et  al.,  1961)  has  further 
detracted  from  attention  to  the  zone  rather  than  the  line . 

Attempts  at  overcoming  the  frustrations  of  the  manual  analysis  of 
HBZ  boundaries  have  been  pioneered  by  the  Meteorological  Service  of  Canada 
(Anderson  et  al.,  1955).   Their  frontal-contour  model,  comprising  three 
fronts  and  four  air  masses  has  been  used  to  assist  in  a  quasi-objective 
approach  to  frontal  analysis.   However,  the  model  places  undue  emphasis 
on  the  labeling  of  air  masses  and  fronts  by  means  of  empirically-determined 
frontal  values  of  various  temperature  parameters,  notably  wet-bulb  potential 
temperature.   Moreover,  the  time  involved  in  the  subjective  manipulations 
necessary  for  the  objective  placement  of  frontal  contours  exceeds  the  human 
capacity  for  anything  but  a  regional  effort,  if  operational  deadlines  are 
to  be  met. 

The  dilemna  of  real-time  urgency  in  frontal  analysis  suggests 
numerical -objective  vice  manual -subjective  or  manua 1 - ob j ec t i ve  analysis 


of  HBZs .   Such  an  approach  has  been  possible  since  the  availability  of 
high  speed  digital  computers  in  the  middle  1950s,  but  either  was  not 
attempted  or  was  unsuccessfully  attempted  until  the  development  reported 
on  here . 

The  earliest  and  only  past  effort  in  numerical-objective  HBZ 
analysis  known  to  this  writer  ,  a  statistical  approach,  was  not  success- 
ful, apparently  because  of  an  over  specification  of  front-associated 
parameters  and  criteria  for  front  location  and  intensity.   It  is  a  natural 
tendency  to  presume  that  the  utilization  of  a  multiple  number  of  variables, 
such  as  pressure,  wind,  temperature,  humidity,  and  their  derivatives,  would 
fare  better  than  any  one  of  them  in  objectively  placing  fronts.   But,  the 
diverse  reactions  of  each  of  these  parameters  to  the  synoptic  and  mesoscale 
processes  which  affect  baroclinicity ,  as  modified  by  topography,  compli- 
cates the  standardization  of  locating  HBZs  by  statistical  formulae  or  by 
dynamical-statistical  means.   The  more  recent  numerical-objective  efforts 
by  the  United  States  Naval  Weather  Service  (Hughes,  1967)  and  Air  Weather 
Service  (United  States  Air  Force,  1965-1968),  the  Canadian  Meteorological 
Service  (Creswick,  1967),  and  British  meteorologists  (Kirk,  1965),  have 
taken  advantage  of  the  research  described  in  detail  here  and  reported  on 
earlier  as  preliminary  accounts  to  the  meteorological  community  (Renard 
and  Clarke,  1965a,  1965b,  1966;  Clarke  and  Renard,  1965,  1966a,  1966b; 
Renard,  1967a,  1967b,  1967c). 


Correspondence  with  personnel  at  the  Third  Weather  Wing,  United 
States  Air  Force,  Offutt  Air  Force  Base,  Nebraska,  and  Travelers  Research 
Center,  Hartford,  Connecticut;  also,  a  Report  Bibliography  prepared  for  the 
writer  by  the  Defense  Documentation  Center  (1965). 
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Thus,  the  middle  1960s  found  numerical  techniques  and  high-speed 
computers  used  to  analyze,  diagnose,  and  prognosticate  the  operationally- 
usable  structure  of  the  atmosphere,  from  near  the  surface  to  stratospheric 
levels,  with  one  notable  exception,  namely  frontal  phenomena.   Without 
philosophizing  as  to  the  inherent  merits  and  deficiencies  of  frontal  anal- 
ysis and  prognosis  to  weather  processes  (Sawyer,  1964),  their  importance  is 
evident  immediately  from  the  most  casual  perusal  of  daily  weather  maps  from 
the  United  States  (civilian  and  military)  or  any  other  nation.   Such  charts 
always  include  front  locations  for  at  least  one  level  in  the  troposphere. 

An  example,  clearly  illustrating  the  need  for  unique  placement  of 
fronts  is  shown  in  Fig.  1.   This  illustration  displays  the  0000  GMT  5  March 
1964  surface  fronts  analyzed  by  sixteen  United  States  and  foreign  civilian 
and  military  weather  centers.   The  non-uniformity  of  concept,  definition, 
and  model  are  manifested  by  fronts  which  cross  at  angles  as  much  as  90 
degrees  and  differ  in  position  by  300  miles  or  more  when  apparently  indi- 
cating the  same  HBZ  boundary.   Such  a  situation  is  intolerable  to  a  scien- 
tific discipline  and,  practically  speaking,  could  be  disasterous  to  a 
civilian  or  military  operation  which  depended  on  such  analyses  for  its 
success.   What  are  the  parametric  criteria  for  placing  each  of  the  fronts? 
Are  they  concerned  with  temperature  and/or  temperature  gradient,  cloud  and/or 
precipitation,  vector  wind  and/or  wind  shift,  pressure  gradient  and/or  pres- 
sure tendency?  Quantitative  answers  to  these  questions  do  not  exist  from 
any  of  the  weather  services  whose  frontal  analyses  are  displayed  in  Fig.  1. 
The  situation  in  1970  differs  very  little  from  that  in  1964  or  1960  (Jacobs- 
Haupt,  1963). 

In  order  to  improve  the  matters  just  described,  experiments  were 
begun  in  1964  to  develop  a  model  for  numerically  specifying  hyperbaroc linic 
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zones,  their  boundaries  and  intensities,  with  the  aim  of  objectively 
specifying  these  phenomena  on  a  real-time  operational  basis.   During  the 
initial  phases  of  the  investigation,  the  decision  was  made  to  adopt  a 
definition  for  the  HBZ  which  is  simple  and  acceptable  to  theoreticians 
and  synopticians  alike.   Such  a  definition  is  that  implied  for  the 
frontal  zone  in  the  American  Meteorological  Society's  Glossary  of 
Meteorology  (Brooks  et  al.,  1959).   In  essence,  the  frontal  zone  is 
described  as  a  three-dimensional  zone  or  layer  of  large  horizontal 
density  gradient  bounded  by  surfaces  (i.e.,  frontal  surfaces)  across 
which  this  gradient  is  discontinuous.   The  intersection  of  a  frontal 
surface  with  another  identifiable  surface  or  plane  (e.g.,  a  constant- 
pressure  surface  or  the  x,  p  plane)  generates  a  line  called  the  front. 
Such  a  definition  agrees  with  that  proposed  by  Taljaard  et  al.  (1961)  in 
their  excellent  appraisal  of  the  state  of  frontal  analysis  in  the  early 
1960s.   Theirs  is  still  a  very  timely  estimation.   Since  the  concept 
stated  in  the  Glossary  of  Meteorology  clearly  isolates  density  gradient, 
and  thus  temperature  gradient,  as  the  essential  front- location  parameter, 
the  measure  of  frontal  phenomena  and  baroclinicity  are  of  the  same  nature 
Therefore,  the  hyperbaroclinic  zone,  a  bounded  region  of  relatively  large 
baroclinicity,  is  alternately  referred  to  as  a  frontal  zone  and  the  HBZ 
boundaries  as  fronts.   However,  the  term  front  is  not  limited  only  to  the 
warm-side  boundary  of  HBZs  as  conventionally  analyzed  in  connection  with 
frontal  wave  cyclones  of  the  middle  and  high  latitudes.   Also,  at  this 
point,  the  discontinuity  aspect  of  the  frontal-zone  definition  should  be 
interpreted  as  a  mathematical  specification  not  necessarily  substantiated 
in  the  real  atmosphere  due  to  a  limited  data  network  and  the  scale  of 
analysis . 


Fig.  1.   Surface  fronts  at  0000  GMT  5  March  1964,  as  trans- 
cribed from  analyses  of  the  following  international  and  national 
civilian  and  military  weather-analysis  centers: 

National  Meteorological  Center,  United  States  Weather  Bureau; 
Central  Analysis  Office,  Meteorological  Branch,  Department  of 
Transport,  Canada;  Japanese  Meteorological  Agency;  British 
Meteorological  Office;  Zentralamt,  Deutscher  Wetterdienst  ; 
Icelandic  Weather  Bureau;  United  States  Navy  Weather  Centrals  at 
Alameda,  Guam,  and  Pearl  Harbor;  United  States  Navy  Weather  Facilities 
at  Argentia,  Miami,  Norfolk,  Quonset  Point,  San  Diego,  Sangley  Point, 
and  Yokosuka. 
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The  subsequent  chapters  describe  the  development,  characteristics, 
and  applications  of  a  numerical  model  for  locating  and  measuring  the 
intensity  of  real-atmosphere  HBZs .   Hemispheric,  synoptic-scale  case 
studies  provide  a  limited  four-dimensional  (x,  y,  p,  t)  diagnosis  of 
HBZs  on  a  multi-seasonal  basis.   However,  the  reader  may  first  wish  to 
peruse  Appendix  I  for  a  survey  of  the  writer's  preliminary  experiments 
in  numerical  frontal  analysis.   These  first  probes,  somewhat  air-mass 
oriented  and  based  on  wet  bulb  and  equivalent  potential  temperatures, 
were  not  extended  due  to  the  difficulties  of  adequately  computing  and 
portraying  moisture  on  a  hemispheric  scale  as  well  as  the  artificiality 
of  air-mass  labeling.   Relegating  this  section  to  an  appendix  signifies 
the  terminal  nature  of  these  early  experiments.   Following  these  prelim- 
inary experiments ,  research  was  directed  toward  a  more  dynamic  approach 
to  numerical  HBZ  analysis,  utilizing  virtual  potential  temperature  and 
its  derivatives.   The  model  proposed  for  general  operational  use  and 
diagnostic  testing  is  discussed  at  length  in  Chapter  II.   Chapter  III 
illustrates  the  model  by  examples  derived  from  hemispheric  application 
in  a  framework  of  an  existent  operational  analysis  and  forecast  complex, 
namely  the  Fleet  Numerical  Weather  Central,  Monterey,  California.   A  more 
refined  scheme  for  specifying  HBZs,  featuring  high  resolution  derived 
from  vertical  soundings,  is  presented  next  (Chapter  IV)  as  a  prototype 
for  analysis  on  the  subsynoptic  scale.   Chapters  V-VII  return  to  the 
operational  model  appropriate  to  present-day  synoptic-scale  analysis 
(Chapter  II)  and  focus  on  substudies  concerned  with  the  spatial  and 
temporal  consistency  of  numerical  HBZ  analysis  and  the  relation  of  these 
analyses  to  cloud  and  vertical  motion  patterns.   Finally,  the  treatise 
concludes  with  a  resume'  of  the  investigative  results  (Chapter  VIII). 


CHAPTER   II 

HYPERBAROCLINIC-ZONE  ANALYSIS: 
DEVELOPMENT  AND  CHARACTERISTICS    OF  A   NUMERICAL  MODEL 

The   Model 

The  preliminary  experiments  described  in  Appendix  A,  although 
valuable,  merely  suggested  a  continuation  of  the  search  for  a  more 
universal  approach  to  objective  hyperbaroc linic-zone  analysis,  that 
is,  one  void  of  intimate  association  with  air-mass  identification. 

The  definition  of  the  HBZ ,  basic  to  the  analysis  model  whose 
development  is  outlined  in  this  chapter,  is  described  in  Chapter  I. 
Since  the  essential  parameter  is  atmospheric  baroc linicity ,  the 
numerical  specification  of  the  HBZ  proceeds  from  the  definition  of 
the  magnitude  of  total  baroclinicity  about  a  local  point,  which  is 
given  by  Saucier  (1955)  as 

N  =  g|v  In  9*| .  (1) 

N  represents  the  number  of  unit  solenoids  per  unit  area  in  a  plane 
determined  by  the  near-vertical  three  dimensional  pressure  gradient 
and  a  normal  to  the  isolines  of  virtual  potential  temperature  (0*) 
on  a  constant  pressure  surface;  g  is  gravity.   Since  V9*  has  mathe- 
matical properties  similar  to  V  In  9*  while  offering  greater  facility 
in  handling  of  data,  the  former  is  used  in  the  development  that  follows 


10 
Symbolizing  the  unit  vector  normal  to  the  line  bounding  the  HBZ 
as  n  ,  a  parameter  for  locating  the  line  on  a  constant  pressure  surface 

may  be  defined  by 

i    i  -»      d20* 
DG0*t  =   -   V|V9*|-n  =  -  ^-^.  (2) 

Thus,  DG0*  is  a  directional  second  derivative  of  9*,  constrained  to  lie 
along  n.  . 

With  reference  to  the  identity  in  Eq .  2,  the  direction  implied 

— > 
by  n   is  unknown  initially,  as  the  position  and  orientation  of  the 

front  is  the  quantity  being  sought.   However,  since  the  HBZ  tends  to 

follow  an  isotherm  channel,  a  suitable  approximation  to  n.  is  given  by 

the  unit  vector  perpendicular  to  isotherms.   Therefore,  with  a  view 

toward  quantitative  applications,  the  HBZ-location  parameter  on  a 

constant  pressure  surface  is  redefined  as 

DGe*  =  -  V|V9*|  -nQ^  =  -  VG0*-nQ^  =  -  -^-^  >       (3) 


^n9* 


where    n        is    a   unit   vector   along  V9*, and 


G9*  =    |V9*|    =   d9*/dnQSt.  (4) 

A   One -dimensional   Schematic 
Analysis    of  DG9* 

A   one -dimensional   representation   of   9*,    G9*,    and   DG0*  which 

simulates  the    real   atmosphere    is    an  aid    in  describing   the    synoptic   and 

thermo-physical   properties    of  DG0*.      Fig.    2    presents    continuous 

schematic    profiles    of   9*  and    its    derivatives    oriented   along  V9*  on   an 

isobaric   surface.      The    profiles    are   slightly   asymmetric   about    location 

0.       Points    1   and   2    in   Fig.    2a   are    locations    of   relative   maxima    in   the 
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Fig.  2.   A  continuous  schematic  profile  of  0*  and  the  G0*  and 
DG0*  profiles  derived  from  it. 

a.  0*;  theoretical  profile  with  first-order  discontinuity  also 
shown  in  the  vicinity  of  points  1  and  2  (dashed  line). 

b.  G0*,  derived  from  Fig.  a. 

c.  DG0*,  derived  from  Figs,  a  and  b. 
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curvature  of  the  9*  profile;  theoretically,  such  points  are  first-order 
discontinuities  in  9*. 

The  expression  G0*  is  a  measure  of  the  baroc linicity  (Fig.  2b), 
which  is  maximized  at  location  0  in  the  9*  profile.   The  inflection 
points  (1,  2)  in  G0*  are  suitable  boundaries  of  the  zone  containing 
relatively  large  magnitudes  of  baroclinicity ,  i.e.,  the  HBZ .   A 
higher-order  derivative  of  9*  is  required  to  give  a  clear  and  quanti- 
tative indication  of  the  boundaries  of  the  HBZ.   This  is  shown  by  the 
DG9*  profile  in  Fig.  2c_.   The  maximum  and  minimum  values  of  DG0* 
indicate  warm-side  and  cold-side  boundaries  of  the  HBZ,  respectively, 
with  the  zero  point  representative  of  the  location  of  maximum  baro- 
clinicity (Fig.  2b).   Though  magnitude  of  baroclinicity  is  not 
explicitly  given  by  the  front- location  parameter,  DG9*,  a  relatively 
large  spatial  change  in  V9*  implies  the  existence  of  high  values  of 
baroclinicity  in  the  associated  HBZ.   Statistical  and  synoptic  veri- 
fications of  this  relation  are  cited  in  this  and  later  chapters. 

Physical  significance  of  DG9*  may  be  seen  from  the  relations 
of  geostrophic  thermal  wind,  V  ,  and  the  isobaric  virtual  temperatures, 


T*  and  9*. 


and 


V  I  =  I^YSI     _£_  dT*  _   £  d  In  9* 
T1   '  '&  '      fT*  bnTic~        f  an9* 


^T  i  _   g  a2  In  9* 

^9*     f      an     2  ' 
9* 


(5) 


(6) 


Therefore,  the  shear  component  of  the  thermal  vorticity,  given  by  the 
shear  of  the  thermal  wind  in  Eq .  6,  also  exhibits  a  maximum  in  magnitude 
at  each  of  the  HBZ  boundaries. 
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Analysis  Characteristics  of  the  Model 

In  order  to  demonstrate  the  effects  of  numerical  differen- 
tiation and  interpolation,  and  scale  on  the  location  and  intensity 
of  HBZs ,  three  hypothetical,  but  realistic,  9*  profiles  are  presented 
(Fig.  3)  and  interpreted.   The  profiles  are  assumed  to  be  oriented  in 
the  direction  of  V9*,  here  considered  to  be  colinear  with  the  I  axis 
of  the  grid.   Data  are  specified  for  whole-numbered  grid  points  only 
(Table  1),  consistent  with  numerically-analyzed  observations.   Although 
the  distance  between  grid  points,  d,  is  somewhat  arbitrary  for  the 
arguments  set  forth  in  this  chapter,  it  will  be  of  great  practical 
value  for  what  follows  to  consider  d  =  381  km  . 

The  temperature  scale  in  Fig.  3  is  symbolic  only,  but  as  in 

the  real  atmosphere,  the  warm  air  is  shown  to  be  less  baroclinic  than 

the  cold  air  in  all  three  cases.   The  profiles  are  linear  in  the  warm 

air  but  depart  slightly  from  such  in  the  cold  air.   The  hyperbaroclinic 

zones,  separating  the  coldest  from  the  warmest  air,  are  bounded  by 

break  points  typical  of  nature's  version  of  first-order  discontinuities 

on  the  synoptic  scale.   These  maximum-curvature  points  in  the  modeled 

profiles  are  labeled  as  F   and  F,  for  boundaries  on  the  warm-  and  cold- 
r  w      k 

air  sides  of  each  of  the  HBZs,  respectively.   The  profiles  are  labeled 
N  (narrow),  M  (medium)  and  W  (wide),  descriptive  of  HBZ  widths  of  one, 
two,  and  three  mesh  lengths,  respectively.   However,  for  all  three  pro- 
files, the  HBZ  AG*  is  10  C .   This  results  in  HBZ  gradients  of  10  C 
(381  km)"1  for  N,  5  C  (381  km)"1  for  M  and  3.3  (381  km)"1  for  W. 


381  km  is  the  mesh  length  at  60  deg  lat  in  the  FNWC  grid 
representation  of  the  Northern  Hemisphere. 
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The   aim   of   this    chapter    is    to   simulate ,    in  a    systematic   but 

simple    fashion,    numerical   differentiation  and    interpolation   of    tempera- 
ture  data   using   computations    of   the    type    common   to   meteorological  analysis 
The   modeled   thermal   data   are    considered   numerically   processed   grid-point 
temperatures.      Cases    of   definite    physical-synoptic    significance   are    sel- 
ected   for    simulation,    as,    for   example,    numerical    processing    in   relation 
to    intensity   of    frontal   zones,    mesh    length,    closely  adjacent    fronts, 
frontal    slopes,    and    the    use    of    layer  vice    level    in   the   DG0*  computations. 

The    Simulation   Experiments 

Frontal-zone    Intensity 

Variation  in  HBZ  width--First ,  the  9*  data  in  Fig.  3  are  used  to 
generate  G0*  values  by  employing  the  Lagrangian  3-  and  5-point  numerical 
differentiation  formulae  (Jennings,  1964)  for  centered  differences.  For 
the   3-point    computation   at    the   central   grid    point,    I    ,    the    formulation   is 

l(e*  ,  -  e*    )| 

while    the    5-point    G0*   is    given  by   the    following   expression: 

1 8  [o*      -  0*    ]  -  e*     +  e*    | 

rc.±  .    o+l  o-l o+2  o-2  ,e* 

G0* ^  .  (8) 

G0*  grid-point  data  are  then  used  to  generate  the  DG0*s ,  again  using 

both  the  3-point  and  5-point  numerical  differentiations. 

In  order  to  show  relative  magnitudes  of  frontal  parameters  in 

and  near  the  N,  M,  and  W  frontal  zones,  Table  1  lists  the  grid-point 

9*s  (from  Fig.  3)  and  the  computed  values  of  G0*  and  DG0*.   It  is  to 

be  noted  that  the  magnitudes  of  GQ*   and  DG0*  at  and  between  F  and  F. 

w      k 
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are  similar  to  those  found  in  frontal  zones  on  the  daily  operational 

charts  (Chapter  III). 

The  foremost  question  to  which  the  profiles  and  tabular  data  are 
directed  is  as  follows.   Assuming  the  modeled  positions  of  the  HBZ  bound- 
aries to  be  at  F  and  F,  ,  where  are  the  frontal  positions  as  located  by 
w      k 

numerical  differentiation  and  interpolation  of  the  grid-point  temperature 
data?   To  answer  this  question,  Bessel's  central  difference  polynomial 
interpolation  formula  (Jennings,  1964)  is  utilized  to  interpolate  the 
location  of  DG6*  maxima  and  minima  to  the  nearest  .1  mesh  length  (i.e. 
about  .3  deg  lat  for  d  =  381  km).   These  points  are  representative  of  the 
warm-side  boundary  (conventional  front)  and  the  cold-side  boundary  of  the 
HBZ,  respectively. 

Setting  DG0*  =  y,  the  formulation  for  DG0*  is 

y  "^y 

«*\h,  -  yWs  =  -\^  +  BjCyj-y,,)  *  B^-y^y.^)       (9) 

+  B3[(y2-yi)  -  2(yi-y0)+  (y^y^)], 
where  subscripts  -1,  0,  1,  and  2  indicate  grid-point  locations  and  s  is 
the  fractional  distance  between  1  =  0  and  1  =  1.   The  constants,  B..  ,  B~  , 
and  B„  are  functions  of  s  in  the  form 

Bl  =  S  ~\> 

B2=^f^,  (10) 

s(s-l)(s-|) 
and  B3  =  g . 

Consistent  with  the  number  of  grid  points  used  in  the  5-point  numerical 

differentiation  formula  (Eq.  8)  the  interpolation  formula  (Eq .  9)  is 

truncated  with  the  term  which  is  cubic  in  s  (third  degree  polynomial). 
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The    numerically-interpolated    locations    of   the   HBZ   boundaries 
and   their  attendant   DG9*  values    are   also   shown   in   Table    1.      For 
example,    in   the   case    of   the   N   profile    the    interpolated    frontal-zone 
boundaries   are    interpolated   to  be   at   7.41  and   5.41    (7.31  and   5.51) 
for   the   3-point    (5-point)    computation  with   DG9*  values    of  +0.21   and 
-0.17    (+0.32   and    -0.27),    respectively.      By   comparison,    the   modeled 
N-profile   data    indicate      HBZ   boundaries    at    17   and    16  with  DG0*  values 
of  +0.30   and    -0.23,    respectively. 

The    following   conclusions    are    pertinent    to  Table    1: 

1.  The    computed   narrow    (i.e.    strong)    frontal   zone    is    broadened 
most   relative    to   the   modeled   zone  width,    both   for    the   3-   and   the    5-point 
differentiation  schemes.      In   fact,    the   weak   frontal   zone    is    fore- 
shortened,   not   broadened,    in   the    5-point    case. 

2.  The   3-point   differentiation   scheme   broadens    the   HBZs   more 
than    the    5-point   approximation.      Specifically,    the   computed   HBZ  widths 
are   broadened    by   the    following    percentages    relative    to   the   modeled   zone 
widths  : 

N  M  W 

3-pt  100  15  0 

5-pt  80  10  -20 

3.  The  broadening  effect  is  up  to  20  percent  greater  on  the 
cold-side  boundary  of  the  HBZ.  It  is  to  be  noted  that  the  cold-side 
boundary  borders  on  the  most  baroclinic  of  the  two  air  masses. 

4.  The  maximum  grid-point  magnitudes  of  5-point  G0*  values 
exceed  those  of  the  3-point  scheme  and  may  be  greater  than  the  values 
implied  in  the  modeled  HBZs  of  Fig.  3,  as  noted  below  for  figures 
extracted  from  Table  1. 
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G0*  (C[100  km]"1) 
HBZ         Mode  led  Computed 

5-pt  3-pt 

N  2.63  1.79  1.58 

M  1.31  1.47  1.31 

W  .87  1.00  .92 

Hence,  the  5-point  GQ*   values  are  more  desirable  in  order  to  yield 

frontal  DG0*s  which  clearly  outline  the  HBZs . 

5.   The  maximum  DG0*  is  greater  in  magnitude  than  the  minimum 
DG0*  for  a  given  HBZ,  a  consequence  of  greater  baroclinicity  in  the 
colder  of  the  two  air  masses . 

Further,  note  that  if  the  modeled  HBZ  boundary  in  Fig.  3  is 
located  between  grid  points  the  broadening  and  magnitude  problems  are 
increased.   For  example,  in  the  case  of  the  N  profile,  the  broadening 
is  quantitatively  increased  up  to  nine-tenths  of  a  mesh  length  by 
placing  the  HBZ  boundaries  at  6.01  and  6.11  while  maintaining  the  9* 
value  at  7.01.   Thus,  with  a  HBZ  width  of  0.11  in  the  N  case,  the 
broadening  percentage  increases  to  1700  percent  for  the  5-point  dif- 
ferentiation.  However,  this  circumstance  is  highly  unlikely  since  a 
HBZ  width  <  40  km  is  unrealistic  on  the  scale  of  present-day  obser- 
vations and  analyses.   In  case  of  the  M  and  W  frontal  profiles  the  HBZs 
could  be  reduced  to  a  zone  lying  between  5.91  and  6.11  and  4.91  and 
6.11,  respectively,  without  a  change  in  grid-point  temperature  data. 
This  would  result  in  a  percentage  broadening  of  1000  percent  and  100 
percent,  respectively,  for  the  M  and  W  profiles  in  the  5-point  case. 
However,  nonmeteorological  9*  profiles  would  result  for  the  M  and  W 
cases  if  the  grid-point  temperature  data  remained  unchanged. 
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Variation  In  HBZ  AQ*- -Frontal-zone  intensity  (i.e.  average  A9* 

within  the  transition  zone)  is  dependent  on  the  change  of  9*  in  the 
zone  (i.e.  A9*)  as  well  as  distance  (An   )  through  which  A9*  occurs. 
In  order  to  demonstrate  the  effect  of  variations  of  V9*  in  the  HBZ  due 
to  changes  in  A9*,  the  A9*  is  reduced  to  5C  for  the  N  profile  only  of 
Fig.  3.   This  results  in  a  modeled  transition-zone  gradient  equivalent 
to  that  of  the  M  profile  in  Fig.  3.   Table  2  shows  the  9*,  G0*,  and 
DG0*  tabulations  and  broadening  effects  for  the  3-  and  5-point  differ- 
entiations in  the  same  format  as  Table  1. 

The  broadening  effect  on  the  cold-side  is  greater  than  that  on 
the  warm-side  boundary  of  the  HBZ,  an  effect  similar  in  kind  to  the  N 
profile  (A9*  =  IOC)  even  though  the  computed  G0*s  and  DG9*s  resemble  the 
W  case  of  Table  1.   The  overall  percentage  broadening  is  about  the  same 
as  the  N  case  in  Fig.  3,  namely  110  percent  for  the  3-point  and  70  per- 
cent for  the  5-point  differentiations,  again  showing  the  desirability  of 
the  latter.   In  view  of  these  effects,  it  appears  that  it  is  the  frontal- 
zone  width  more  than  the  magnitude  of  A9*  which  fosters  significant 
computational  broadening  of  the  HBZ.   That  is,  the  real-atmosphere 
narrow  transition  zones  are  broadened  relatively  more  than  the  wide 
zones  by  numerical  processing,  regardless  of  the  magnitude  of  A9*  across 
the  HBZ. 

Both  Tables  1  and  2  suggest  the  insignificance  and  noisiness  of 
DG0*  data  outside  the  HBZ  where  values  generally  become  <  | .05|c(100  km)" 
within  one  mesh  length  of  HBZ  boundaries.  Also,  the  zero  DG0*  values  lie 
half-way  between  the  maximum  and  minimum  values  (to  the  nearest  .1  d)  for 
all  5-point  cases  in  Tables  1  and  2.  These  facts  provide  guidance  in  the 
development  of  an  operational  DG0*  product,  as  discussed  in  Chapter  III. 
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Table  2.   Grid  point  values  of  9*  (C),  and  computed  G9* 
(C[100  km]-1)  and  DG9*  (C[100  km]-2)  at  intervals  of  II  for  the  N 
profile  in  Fig.  3,  for  the  3-  and  5-point  numerical  differentiation 
schemes.   Subsets  of  DG9*  values  indicate  interpolated  HBZ-boundary 
locations  and  their  attendant  DG9*  maxima  and  minima  (circled  numbers). 
HBZ  A9*  =  5C.   F   (F,  )  is  the  position  of  the  warm  (cold)  boundary  of 
the  modeled  HBZsY 
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Variation  in 
Mesh  Length 

In  order  to  demonstrate  the  increase  of  resolution  in  locating 
fronts  with  a  decrease  in  mesh  length,  the  N  and  M  profiles  of  Fig.  3 
were  considered  in  connection  with  a  doubling  of  the  number  of  grid 
points.   Linear  interpolation  was  used  to  obtain  the  0*  values  at  the 
positions  intermediate  to  the  whole  numbered  grid  points  in  Fig.  3  (i.e., 
at  I»  ,.  I,    I  re,  etc.)  with  slight  exception  in  the  M  case.   Here,  a 
deviation  from  linearity  was  employed  to  show  maximum  baroclinicity  at 
the  geometric  center  of  the  HBZ .   The  5-point  differentiation  scheme 
and  AG*  =  IOC  are  common  to  the  two  cases  computed. 

N  Case—Table  3  gives  the  pertinent  information  on  this  case. 
The  positioning  of  the  HBZ  boundaries  is  precise,  hence  the  broadening 
effect  is  zero  compared  to  80  percent  broadening  with  data  at  whole- 
numbered  grid  points  only.   The  maximum  G0*  value  has  increased  by 
almost  70  percent  compared  to  the  N  case  in  Table  1  with  G0*  exceeding 
the  model  value  of  2.63  C  (100  km)   ;  the  associated  DG0*s  have  tripled, 
a  desirable  result  for  quantitative  specification  of  the  HBZ  boundaries. 

M_Cas_e- -Table  3  also  gives  the  details  on  the  M  case.   Much  like 
the  N  case  the  positioning  is  excellent,  with  the  HBZ  foreshortened  by 
two-tenths  mesh  length  (10  percent)  compared  to  the  modeled  width,  mostly 
due  to  conditions  on  the  cold-side  boundary.   Again,  magnitudes  of  G0* 
and  DG0*  have  increased  compared  to  their  counterparts  in  Table  1. 

The  obvious  conclusion  is  that  a  mesh- length  reduction  resolves 
the  computed  position  and  intensity  of  the  HBZ.   Any  nominal  fore- 
shortening due  to  numerical  processing  is  not  only  tolerable  but  desir- 
able due  to  the  smoothing  inherent  in  generating  the  data  at  grid  points 
from  the  reported  observations. 
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Table  3.   Grid  point  values  of  0*  (C),  and  computed  G0* 
(C[100  km]-1)  and  DG0*  (C [ 100  km]"2)  at  intervals  of  .51  for  the 
N  and  M  profiles  in  Fig.  3,  for  the  5-point  numerical  differentiation 
scheme  only.   Subsets  of  DG0*  values  indicate  interpolated  HBZ-boundary 
locations  and  their  attendant  DG0*  maxima  and  minima  (circled  numbers). 
HBZ  A0*  =  IOC.   F   (F  )  is  the  position  of  the  warm  (cold)  boundary  of 
the  modeled  HBZs .W 
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Closely  Adjacent  Fronts 

A  typical  occurrence  is  that  which  may  be  described  as  closely 
adjacent  fronts.   Operationally,  the  fronts  are  referred  to  as  primary 
and  secondary;  the  former  designates  the  most  equatorward  of  the  two 
fronts.   Usually  both  frontal  zones  are  of  the  cold-front  variety  and 
the  HBZ  V9*  associated  with  the  secondary  is  generally  weaker  than  that 
of  the  primary. 

Two  examples  are  given  to  illustrate  this  case.   (1)  A  combin- 
ation of  N  and  M  profiles  with  an  intervening   subbaroclinic  zone  or 
air  mass  of  one-mesh  length  dimension;  appreciable  baroclinicity  occurs 
in  the  air  separating  the  N  and  M  transition  zones.   (2)  Same  as  (1) 
except  an  intervening  air  mass  of  two  mesh  lengths. 

Table  4  illustrates  the  multiple-front  problem  with  values  of 
G9*  and  DGO*. 

One-mesh  air  mass--In  the  one-mesh  air-mass  case  the  broadening 
effects  noted  with  N  and  M  profiles  (Table  1)  are  reduced.   In  other 
words,  the  N  or  primary  frontal  zone  shows  40  percent  vice  80  percent 
broadening  and  the  k-side  is  more  accurately  positioned  than  the  w-side, 
compared  to  the  uncomplicated  profiles  in  Fig.  3.   The  M  or  secondary 
frontal  zone  is  foreshortened  by  40  percent  (vice  10  percent  broadening 
in  Table  1)  with  the  k-side  again  placed  most  accurately.   The  secondary, 
especially  its  w-boundary,  behaves  as  a  very  weak  front  here. 

When  viewing  the  intensities,  the  w-boundary  of  the  primary  is 
similar  to  profile  N  in  Table  1,  while  the  k-boundary  is  reduced  in 
magnitude.   The  w  boundary  of  the  secondary  is  close  to  the  magnitude  of 
the  noise  level  criterion  described  earlier,  hence   the  secondary  is 
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Table  4.   Grid  point  values  of  0*  (C),  and  computed  G9* 

(C[100  km]"1)  and  DG0*  (C [ 100  km]'2)  at  intervals  of  II  for  the 

multiple-front  case  with  air-mass  dimensions  of  one  and  two  mesh 

lengths.   The  primary  (secondary)  front  is  represented  by  the  N  (M) 

HBZ  profile  in  Fig.  3;  5-point  differentiation  scheme  only.   Subsets 

of  DG9*  values  indicate  interpolated  HBZ-boundary  locations  and  their 

attendant  DG9*  maxima  and  minima  (circled  numbers).   HBZ  AG*  =  IOC. 

F   (F.  )  is  the  position  of  the  warm  (cold)  boundary  of  the  modeled  HBZs 
w   k 
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poorly  resolved.   However,  the  secondary's  k-boundary  is  well  defined 
in  magnitude. 

It  should  be  noted  that  a  one-mesh  air-mass  distance  is  not 
very  realistic.   Further,  as  in  the  simple  structures  in  Fig.  3,  the 
numerical  processing  effect  causes  all  HBZ  widths  to  approach  a  medial 
value.   The  relatively  narrow  frontal  zones  broaden,  the  wide  ones  tend 
to  foreshorten. 

Two-mesh  air  mass--When  the  air-mass  zone  is  widened  to  two 
mesh  lengths  both  the  primary  and  secondary  zones  behave  as  if  the 
other  were  not  present.   In  the  two  mesh  case,  the  primary  HBZ  is 
broadened  by  60  percent  vice  80  percent  while  magnitudes  are  approxi- 
mately the  same  as  those  of  the  N  profile  in  Table  1.   The  secondary 
front  is  precisely  located  and  magnitudes  agree  closely  with  the  simple 
M  profile  in  Table  1. 

It  is  obvious  from  this  experiment  that  the  position  and  intensity 
of  each  of  the  HBZs  become  increasingly  independent  as  the  subbaroclinic 
zone  separating  them  widens. 

Vertical  Slope  of 
HBZ  Boundaries 

What  is  the  relationship  of  the  slope  of  the  HBZ  boundaries  and 

the  vertical  variation  of  HBZ  widths  as  caused  by  numerical  processing? 

There  are  two  cases  of  interest:   transition -zone  baroclinicity  weakens 

with  increasing  height;  and  transition-zone  baroclinicity  intensifies 

with  increasing  height.   These  are  applied  to  the  simple  profiles  (N,  M, 

W)  shown  in  Fig.  3.   The  frontal-zone  intensity  is  a  function  of  the 

width  only. 
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HBZ  weakens  with  increasing  height--Fig.  4a  illustrates  the 

situation  for  a  modeled  warm-boundary  slope  of  approximately  1/240, 
with  frontal-zone  characteristics  changing  from  N  to  M  to  W  from  1000 
to  850  to  700  mb,  respectively.   This  results  in  an  associated  cold- 
boundary  slope  of  1/480.   The  1000  to  700  mb  warm-boundary  slope, 
determined  from  numerically-computed  data,  is  1/312,  which  is  a  23 
percent  reduction  from  the  modeled  slope.   At  the  cold  boundary  the 
average  slope  increases  relative  to  the  hypothetical  slope  by  26  per- 
cent (from  1/480  to  1/384)  in  the  layer  1000  to  700  mb . 

The  computed  transition-zone  widths,  in  comparison  with  the 
modeled  values,  range  from  greater  (i.e.  broadening)  at  1000  mb  to 
less  (i.e.  foreshortening)  at  700  mb.   The  overall  effect  is  that  of 
reducing  the  range  of  HBZ  widths  relative  to  the  model  values.   Boundary 
values  of  DG0*  decrease  with  height  consistent  with  the  vertical  gradient 
of  transition-zone  baroc linicity . 

HBZ  intensifies  with  increasing  height--Fig.  4_b  illustrates  the 
case  of  the  frontal  intensity  changing  from  W  to  M  to  N  with  increasing 
height  in  the  layer  1000  to  700  mb .   Again,  1/240  is  taken  as  the  slope 
of  the  warm-air  edge  of  the  HBZ;  the  computed  value  is  1/168.   The 
modeled  cold-side  boundary  of  the  HBZ  is  constrained  to  be  vertical  with 
the  computed  slopes  ranging  from  1/120  in  the  low  layer  to  1/72  in  the 
upper  layer.   Thus,  both  the  broadening  effect  and  magnitude  of  DG0* 
maxima  increase  with  increasing  height. 

Use  of  Layer-mean  Temperature, 
9~*,  in  the  DG9*  Analysis 

Low- troposphere  mean  temperature  (i.e.  is obaric- layer  thickness) 

fields  enjoy  widespread  use  for  locating  fronts  operationally.   Here, 
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700 


Fig.  4.   Numerically-computed  HBZ  boundaries  (dashed  lines)  in 
relation  to  the  modeled  boundaries  (solid  lines);  HBZ  AG*  =  10  C. 

a.  Slope  of  modeled  warm-side  boundary:   1/240;  width  of  HBZ 
increases  with  height 

b.  Slope  as  in  a;  width  of  HBZ  decreases  with  height 
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typical  1000/700  mb  mean- layer  temperature  (0*)  profiles  are  generated 
from  information  in  Fig.  3  and  reasonable  vertical  slopes.   Only  the 
5-point  derivative  scheme  is  used.   The  following  cases  are  considered 
for  HBZ  slopes  of  1/120  and  1/240:   (1)  1000  and  700  mb  --N-type  HBZs 
both  levels;  (2)  1000  mb  --N-type  HBZ,  700  mb  --M-type  HBZ;  and  (3) 
1000  mb  --M-type  HBZ,  700  mb  --N-type  HBZ. 

9*,  at  each  grid  point,  is  obtained  as  an  arithmetic  average  of 
the  1000-  and  700-mb  information  where  the  appropriate  700-mb  profile 
(N  or  M)  is  displaced  one  (two)  mesh  length (s)  toward  the  cold  air  for 
the  1/120  (1/240)  slope  case,  respectively.   In  addition,  3C  to  8C  are 
added  to  the  9*  values  at  700  mb  in  order  to  simulate  the  normal  increase 
of  9*  with  height,  with  the  higher  (lower)  temperatures  used  in  the  warm 
(cold)  air  mass.   The  addition  is  made  in  such  a  way  as  to  preserve  the 
transition  zone  V9*  at  each  isobaric  level.   Such  a  manipulation  of  the 
data  agrees  with  the  observable  result  that  the  cold  air  is  both  more 
unstable  and  baroclinic  than  the  warm  air  in  the  lower  part  of  the 
troposphere.   Table  5  shows  9*,  GQ* ,    DG0*,  and  interpolated  HBZ  boundaries 
similar  to  previous  sections. 

N-type,  both  at  1000  and  700  mb ;  slope:   1/120--This  case  simu- 
lates a  relatively  strong  cold  front  with  steep  slope.   The  numerical 
positioning  of  the  HBZ  boundary  is  exact,  resembling  the  M  profile, 
single  level,  while  the  magnitudes  of  DG9*  are  intermediate  to  those 
for  N  and  M  at  a  single  level  (Table  1).   The  HBZ  width  for  9*  is  the 
distance  from  the  1000  mb  warm-side  to  the  700  mb  cold-side  HBZ  boundary, 
thus  essentially  doubling  the  hypothetical  zone  width  compared  to  either 
1000  or  700  mb. 
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Table  5.  Grid  point  values  of  1000/700  mb  9*,  and  computed  G0* 
(C[100  km]"1)  and  DG0*  (C [  100  km]"2)  at  intervals  of  II  for  N,  M,  and  W 
profile  data  derived  from  Fig.  3  and  frontal  slopes_of  1/120  and  1/240; 
5-point  differentiation  scheme  only.  Subsets  of  DG0*  values  indicate 
interpolated  HBZ-boundary  locations  and  their  attendant  DG0*  maxima  and 
minima  (circled  numbers).  HBZ  A9*  =  10C.  F  (F,)  is  the  position  of  the 
warm    (cold)    boundary   of   the   modeled   HBZs . 
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The  nature  of  the  vertical  slope  is  such  as  to  give  the  effect 

of  an  M  case  in  both  frontal  intensity  and  width  since  the  superposition 

of  F,  1000  and  F  700  tends  to  obliterate  the  identity  of  each  of  these 
k  w 

HBZ  boundaries  when  using  the  mean  temperature  field. 

N-type  both  at  1000  and  700  mb ;  slope:   1/240--Positioning  is 
poor,  with  the  numerical  1000  mb  front  displaced  by  one-third  mesh  length 
toward  the  warm  air  and  the  computed  700  mb  cold-side  HBZ  boundary  ex- 
tended back  into  the  cold  air  by  eight-tenths  mesh  length,  a  36  percent 
increase  over  the  modeled  HBZ  width.   Magnitudes  of  DG9*  are  reduced  to 
those  of  the  W  category,  but  the  broadening  is  typical  of  the  N  category 
(Table  1),  partly  due  to  the  greater  geometric  separation  between  the 
1000  mb  warm-side  and  700  mb  cold-side  of  the  HBZ  compared  to  the  1/120 
slope  case. 

For  further  slope  reduction  the  minor  DG9*  minima  and  maxima,  only 

evident  as  noise  level  values  at  16  and  15  in  Table  5,  become  stronger, 

until  two  separate  and  significant  HBZs ,  F  1000  to  F.  1000  and  F  700  to 

w         k         w 

F.700,  result.   Limiting  such  occurrences  is  the  fact  that  in  the  real 

atmosphere,  the  N-type  fronts  have  large  slope  (i.e.  >  1/120). 

1000  mb  --N-type;  700  mb  -- M-type  ;  slope:   l/120--Position  and 

intensity  of  the  warm-side  boundary  of  the  HBZ  behaves  as  the  M  type  but 

the  cold-side  boundary,  F  700,  appears  as  a  doublet  with  small  magnitudes 

in  each  case.   The  cold  boundary  complexity  results  from  the  superposition 

of  a  strong  F,  1000  and  a  moderate  F  700  front, 
k  w 

1000  mb-- N-type;  700  mb  --M-type;  slope:   l/240--In  this  case, 
due  to  increased  separation  of  the  HBZs  at  each  level,  the  1000  mb  and 
700  mb  frontal  zones  manifest  themselves  as  separate  entities.   With 
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further  reduction  in  slope  the  F,  1000  and  F  700  boundaries  become 

r       k  w 

increasingly  distinct.   In  this  case,  the  F  700  value  of  +0.03  may 
be  disregarded  as  noise. 

1000  mb  --M-type;  700  mb  -- N-type  ;  slope:  l/120--Some  rather 
unexpected  results  show  up  for  this  case.  The  warm  boundary  is  fore- 
shortened, whereas  a  broadening  effect  occurs  at  the  cold  boundary. 

Magnitudes  are  typical  of  the  M  case  (Table  1).   Apparently,  the  F  700 

w 

position  is  the  cause  of  the  peculiar  positioning  of  the  1000  mb  front 
giving  the  properties  of  the  W  profile  to  the  warm  boundary,  while  N 
properties  are  noted  at  the  rear  boundary  where  F  1000  and  F  700  are 

K  K 

superimposed . 

1000  mb  -~ M-type;  700  mb  -- N-type;  slope:   1/240--Comment  is 
similar  to  that  immediately  above  with  the  added  feature  of  a  more 
intense  cold-side  HBZ  boundary  compared  to  the  warm-side. 

Summary  of  the  1000/700  mb  layer  simulation—Use  of  9*  vice  0* 
appears  to  offer  advantages  for  the  front- location  model.   The  broaden- 
ing effect  at  the  HBZ  warm-air  edge  is  reduced,  hence  this  entity  is  more 
accurately  located  by  numerical  means.   Although  DG6*  values  are  less  than 
DG0*  for  similar  transition  zones,  the  magnitudes  remain  well  above  noise- 
level  values.   With  a  view  toward  operational  analysis  of  DG0*,  it  is  not 
too  serious  that  the  cold-side  boundary  locations  are  complex  or  poorly 
located  due  to  the  smearing  inherent  in  the  thickness  model  since,  by 
convention,  only  the  warm  edges  of  frontal  zones  are  portrayed  on  the 
surface  or  near-surface  isobaric  charts. 

The  occurrence  of  fronts  intensifying  with  height  tends  to  produce 
foreshortening  (Table  5)  of  the  HBZ,  which  helps  to  balance  the  broadening 
effect  due  to  sparse  data. 
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General  Conclusions 

The  general  conclusions  concerning  the  accuracy  of  numerical 
HBZ  analysis  as  derived  from  experimentation  with  the  profiles  in  Fig. 
3 ,  are  as  follows : 

1.  Use  of  the  Lagrangian  5-point  vice  the  3-point  numerical 
differentiation  formula  results  in  the  most  accurate  placement  and 
intensity  of  the  HBZ. 

2.  There  is  a  tendency  toward  uniform  HBZ  widths  in  the 
numerical  analysis  of  DG9*,  that  is,  the  narrower  zones  are  widened, 
the  wider  zones  are  narrowed.   Stated  another  way,  broadening  of  the 
HBZ  is  more  of  a  function  of  scale  than  of  degree  of  baroclinicity  in 
the  zone. 

3.  The  most  intense  or  sharply-focused  HBZ  boundary  is  the 
warm-side  (i.e.,  the  conventional  front),  just  as  in  the  actual  atmos- 
phere, a  result  of  the  relatively  smaller  amount  of  baroclinicity  in  the 
warm  compared  to  the  cold  air. 

4.  Although  there  is  a  definite  tendency  toward  an  unreal 
uniformity  of  computed  transition-zone  widths,  there  is  difficulty  in 
devising  an  empirical  correction  as  a  function  of  the  magnitude  of  the 
DG0*  maximum  since,  by  itself,  frontal  intensity  (V9*)  is  not  a  good 
indicator  of  the  numerically-derived  width.   Rather,  it  is  the  true 
frontal-zone  width  which  is  a  good  indicator  of  the  broadening  effect, 
but  true  width  is  not  known.   However,  it  is  apparent  that  the  frontal 
zones  with  relatively  large  thermal  gradients  are  most  likely  to  be 
broadened  but  the  weaker-gradient  zones  may  be  broadened  or  fore- 
shortened due  to  numerical  processing. 
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5.  Overall,  the  0*  model  is  more  desirable  than  the  9*  model 
since  the  effect  of  using  a  mean  temperature  biases  the  front- location 
properties  toward  those  of  the  wide  and  weak-zone  variety  and  this 
reduces  the  broadening  effect.   However,  when  using  9*,  the  cold  boundary 
of  the  lower-level  frontal  zone  cannot  be  located.   This  is  also  true, 

to  a  large  degree,  for  the  cold  boundary  of  frontal  zones  at  the  upper 
level  of  the  layer. 

6.  Closely  adjacent  frontal  zones  are  likely  to  be  misplaced. 
In  that  case,  foreshortening  and  weakening  of  the  conventional  boundary 
of  the  secondary  is  a  likely  effect. 

7.  Reduction  of  mesh  size  aids  in  resolving  the  true  location 
and  intensity  of  HBZs .   Any  tendency  toward  foreshortening  of  the  HBZ, 
as  mesh  size  is  reduced,  is  also  regarded  as  favorable  considering  the 
smoothing  inherent  in  data  processing  and  the  relative  sparsity  of  data. 

8.  The  vertical  variation  of  HBZ  intensity,  in  terms  of  the 
zone's  width,  contributes  toward  uniformity  of  widths  and  vertical  slopes 
of  the  HBZs . 

It  should  be  emphasized  that  the  general  analysis  characteristics 
of  the  DG0*  model  presented  in  this  chapter  were  reduced  from  a  rela- 
tively simple  set  of  temperature  profiles  and  a  controlled  simulation 
experiment.   The  real  atmosphere  presents  considerably  greater  com- 
plexity in  temperature  distributions  with  a  wide  range  of  concurrent 
conditions  of  frontal-zone  intensity,  width,  slope,  and  contiguity. 
Thus,  it  is  logical  to  consider  next  a  two-dimensional  application  of 
the  DG0*  numerical  analysis  scheme  to  a  hemispheric  portion  of  the  real 
atmosphere . 


CHAPTER  III 

AN  OPERATIONAL  APPLICATION  OF  THE 
HYPERBAROCLINIC-ZONE  ANALYSIS  MODEL 

Chapter  II  summarized  the  development  of  a  numerical  analysis 
model  which  has  potential  use  in  identifying,  locating,  and  diagnosing 
HBZs  in  real  time.   Insight  was  also  given  on  the  effects  of  scale  and 
numerical  differentiation  and  interpolation  on  the  computed  values  of 
baroclinicity  and  zones  of  hyperbaroclinicity .   To  this  point  in  the 
presentation,  neither  operational  merits  and  deficiencies  of  the  DG0* 
analysis  nor  the  effects  of  numerical  data  processing  on  the  scheme 
have  been  considered.   Thus,  it  became  desirable  to  program  the  DG0* 
analysis  scheme  for  use  in  a  computer-oriented  weather  central  with 
international  responsibility.   Such  an  environment,  and  the  one  most 
accessible  to  the  physical  location  of  the  subject  research,  is  that  of 
the  United  States  Navy's  Fleet  Numerical  Weather  Central,  Monterey, 
California  (FNWC) .   There  follows  a  discussion  of  the  FNWC  automatic 
data  processing  program  and  examples  of  two-dimensional  DG9*  analyses 
as  performed  in  that  framework.   Since  field  use  is  directed  toward  the 
analysis  of  the  conventional  front,  numerical-objective  and  manual- 
subjective  versions  of  frontal  analyses  are  discussed  in  this  and  later 
chapters . 
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Data-Processing   Framework   for   the   Operational   DG8*  Analysis 

Germane    to   the    research    on  hyperbaroclinic    zones    is    the   FNWC   data 
processing    framework   in  which    the   DG0*  analysis    is    performed.      Holl,    et 
al.    (1963,    1965)    have   developed   a   so-called   mass    structure   model    for 
numerically    processing  meteorological   data    in  a    physically   consistent 
manner.      The   model    is    used   currently   by  FNWC    (Hughes,    1967).      The    fea- 
tures   of   the   model,   which    relate   directly   to   a   comprehension   of   the 
peculiarities    of   the   analysis    of   baroclinicity ,    are    described   below. 

The   data    used   as    input    to   the   DG0*   program   are    those    objectively 
analyzed   by  FNWC    from  a   Northern  Hemisphere    synoptic-time    coverage    of 
more    than  4000   surface   and    500   radiosonde    observation   stations.      These 
data  are    processed  at    12-hour   intervals    to   produce    objective  analyses   of 
sea-level   pressure   and   height   and   temperature    of   all  mandatory   pressure 
levels    up   to   200  mb ,    as   well   as   dew-point   depressions    at    terrain    level, 
850,    700,    and   500  mb .      Except    for   terrain- level   parameters,    the   analyses 
are    produced   on  a   square    grid,    63   x   63,   wherein   the   equator    is    an 
inscribed   circle. 

The   mesh    length   is    381   km  at   60  deg    lat    on  a   Northern  Hemisphere 

polar    stereographic    projection.      At    other    latitudes,    I,    the   geometric 

mesh   distance    is    381   x  i  - — ; : —  )  km. 

VI  +  sin  60y 

The  mass -structure  processing  of  input  data  is  accomplished  in 
the  following  manner.   A  mean  static  stability  for  each  of  five  layers 
(1000  -  775,  775  -  600,  600  -  400,  400  -  300,  and  300  -  200  mb) ,  at  each 
radiosonde  station,  is  determined  as  a  function  of  the  heights  of  the 
mandatory  pressure  levels  and  sea-level  pressure.   Objective  analyses 
of  the  layer  stabilities,  the  1000-mb  height,  and  the  thickness  of  the 
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layer  1000  to  500  mb  yields  grid  point  values  of  these  key  parameters 

at  each  of  the  3969  grid  points.   By  a  complex  matrix  solution,  grid 
point  values  of  temperature  and  height,  at  each  of  the  mandatory  pres- 
sure surfaces,  are  generated  from  the  fields  of  the  key  parameters.   The 
temperatures  so  produced  are  linear  in  p  at  each  grid  point  (p  = 

pressure;  k  =  R  /c  ;  R,  =  the  gas  constant  for  dry  air;  c  =  the  specific 
dpd  p 

heat  at  constant  pressure)  and  closely  approximate  virtual  temperatures. 
The  objectively-analyzed  fields  of  the  state  parameters,  pressure  (or 
height),  temperature,  and  moisture,  are  then  available  to  compute  other 
meteorological  quantities,  such  as  those  selected  for  numerical  frontal 
ana  lysis  . 

For  computing  derivatives,  the  grid-point  values  are  approximated 
by  employing  the  Lagrangian  5-point  differentiation  formula  for  centered 
differences  along  each  of  the  linear  grid  axes,  I,  J,  with  the  two  com- 
ponents combined  vectorially.   See  Eq .  8,  Chapter  II.   As  a  rule,  the 
first  derivatives  were  processed  with  a  low-pass  filter,  having  a  cut- 
off wavelength  of  four  mesh  lengths  (approximately  850  nautical  miles) 
before  further  differentiation. 

An  Example  of  Two-Dimensional  Analyses  of  DG0* 
Fig.  5  indicates  the  typical  mathematical  relationships 
existing  between  FNWC's  objectively-analyzed  isobaric  fields  of  0*  and 
the  patterns  of  GQ*   and  DG0*  derived  from  it  for  a  5  x  8  mesh- length 
area  (d  =  381  km).   The  illustrations  may  be  regarded  as  simulations  of 
middle-latitude  winter  conditions  in  the  low  troposphere. 

The  plan  views  contain  two  separate  zones  of  concentrated  thermal 
gradients  (i.e.,  HBZs),  the  bounding  lines  of  which  are  shown  in  each  of 
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the  Figs.  5a,  b_,  and  £.   The  field  of  9*  in  Fig.  5a  is  used  to  generate 

the  isolines  of  G0*  in  Fig.  5b.   The  HBZ  associated  with  the  coldest 
temperatures  is  most  baroclinic  with  G9*  exceeding  2.0  C  (100  km) 
along  the  zone's  centrum.   Fig.  5£  shows  the  DG0*  field,  exhibiting 
elliptically  shaped  isolines. 

The  locus  of  points  along  which  V0*  changes  most  rapidly  in  the 
direction  of  n„ ,  defines  the  axes  of  maximal  and  minimal  DGG*.   These 
are  the  warm-  and  cold-side  boundaries  of  the  HBZ,  respectively.   Both 
boundaries  tend  to  parallel  the  isentropes,  especially  the  former  (Fig. 
5a).   It  is  to  be  noted  that  the  axes  of  maximum  and  minimum  V9*  (from 
Fig.  5b)  nearly  coincide  with  the  zero  isolines  in  Fig.  5c_;  these  lines 
represent  the  centrums  of  the  hyper-  and  sub-baroclinic  regions,  res- 
pectively.  The  imperfect  coincidence  of  G0*  axes  and  the  zero  isolines 
of  DG0"  is  a  result  of  the  numerically-computed  finite  difference 
approximations  to  the  analytic  expressions.   Fig.  5d  illustrates  pro- 
files of  9*,  G0*,  and  DG0*  along  line  N  in  Figs.  5a--c_.   The  similarity 
to  the  schematic  profiles  in  Fig.  2  is  obvious. 

Operational  Hemispheric  HBZ  Analysis 
Following  closely  on  the  development  of  the  HBZ  analysis  model, 
FNWC  has  been  producing,  twice  daily,  operational  hemispheric  1000  mb 
numerical  frontal  analyses  for  the  guidance  of  Fleet  Weather  Centrals 
and  Facilities  in  analyzing  the  important  synoptic-scale  low-troposphere 
fronts.   Analysis  criteria,  degree  of  smoothing,  and  other  presentation 
features  have  evolved  during  the  period  1965  to  date. 

In  view  of  the  experiments  in  Chapter  II,  a  mean  virtual  potential 
temperature  for  the  layer  1000/700  mb ,  9*,  was  selected  to  generate  the 
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Fig.  5.   Typical  plan  views  and  profiles  of  9*,  G0*,  and  DG0*. 

a.  0*  isolines,  at  intervals  of  2.5  K.  Positive  and  negative 
axes  of  DG0*  superimposed,  the  latter  as  dashed  lines,  the 
former  as  dashed  lines  with  alternating  circles  and  diamonds 

b.  G0*  isolines,  at  intervals  of  0.2  C  (100  km)"1.   Axes  in  the 
G0*  field  shown  as  dotted  lines.   DG0*  axes  as  in  a. 

c.  DG0*  isolines,  at  intervals  of  .05  C  (100  km)"2,  with  axes 
as  in  a;  G0*  axes  from  b. 

d.  Profiles  of  9*,  G0*,  and  DG0*  along  line  N  in  a,  b,  and  c. 
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front-location  parameter  appropriate  to  the  1000  mb  (or  surface)  chart, 

DG0*.   For  this  purpose,  9*  is  derived  from  the  layer  thickness  and 
assigned  to  the  logarithmic  mean  isobaric  level,  836  mb.   The  next 
subsection  focuses  attention  on  the  merits  and  deficiencies  of  the 
operational  DG0*  analysis  with  illustrations  for  a  typical  winter  day. 
A  statistical  evaluation  of  numerical  and  manual  frontal  analyses  serves 
to  highlight  differences  and  similarities  of  the  two  modes  of  analysis. 
Objective  Frontal  Analysis  for  0000  GMT  6  March  1966 

The  radiosonde  observation  time,  0000  GMT  6  March  1966,  was 
selected  to  illustrate  the  features  of  objective  frontal  analysis.   This 
particular  synoptic  time  was  not  chosen  because  of  the  excellence  of  the 
numerical  frontal  analysis,  but  rather  because  of  availability  of  many  of 
the  charts,  fields,  etc.  (Figs.  6  to  11)  necessary  to  demonstrate  the 
present  status  of  the  analysis  using  the  DG0*  scheme. 

Fig.  6  shows  a  typical  numerically-produced  analysis  of  the 
frontal  parameter  on  a  polar  stereographic  projection  of  the  Northern 
Hemisphere.   Because  layer  temperatures  are  used  and  since  the  operational 

application  is  to  (conventional)  frontal  analysis,  only  isolines  of 

-  -2  -2 

positive  DG0*  are  shown  in  Fig.  6,  starting  with  5  x  10   C  (100  km) 

Values  of  DG0*  lower  than  5  units  fall  in  the  noise  range.   The  axes  of 
relative  maxima  in  the  field  of  the  elliptically-shaped  isolines  repre- 
sent the  positions  of  the  warm-side  boundaries  of  the  1000  mb  HBZs ,  hence 
they  may  be  used  by  the  field  meteorologist  as  a  first  guess  to  locate  the 
significant  surface  fronts.   As  an  example,  the  numerically-analyzed  1000 
mb  fronts  in  the  middle  and  eastern  Atlantic  area  are  indicated  by  the 
superimposed  dashed  lines  which  trace  out  the  DG0*  axes . 
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Fig.  6.   DG9*  (1000/700  mb)  analysis  for  0000  GMT  6  March  1966, 
Isolines  at  interval  of  20  x  10"2  C  (100  km)  ,  starting  with  +5  units. 
DG9*  axes  indicated  by  dashed  lines  for  a  section  of  the  Atlantic  area, 
NMC  fronts  are  shown  with  conventional  symbols.   Alphabetical  desig- 
nators A,  B,  C,  D  refer  to  Fig.  16. 
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What   constitutes   a    significant    front    or  HBZ   boundary  may  be 

determined   by    the    spatial   extent    of   a   DG0*  axis    and    the   attendant   mag- 
nitudes,   as   well   as   by   the    location,   general   orientation,    and   shape    of 
this    field.      The    nearly  meridional  eastern  United   States-North  Atlantic 
and   nearly   zonal  western   Pacific   DG0*    axes    are   examples    of   significant 
fronts    in  Fig.    6.      Extensive    zones   with   DG0*   in  excess    of   25   units 
signify    intense    fronts;    these   are   HBZs    containing    large    temperature 
gradients.      Note    the   northern   Pacif ic--Siberia   and  western   Pacific-- 
China    fronts   with   DG0*  exceeding   25   units    over   a    large   area;   maximum 
values    are   greater   than   65    units. 

The    orientations,    but    not    necessarily   the    positions,    of  axes 
of  DG0*   in  Fig.    6   appear   reasonable    relative    to   the    fronts    usually   seen 
on   the   daily   synoptic   analyses   emanating    from   the   National  Meteorologi- 
cal  Center,    Suitland,   Maryland    (NMC).      However,    air-mass   designations 
(now  mostly   extinct,    anyway)    and   movement    (warm,    cold,    etc.)    are   not 
indicated   explicitly    in   the   DG0*   field;    also,    stage    of  development,   as 
the    occluded   stage,    is    difficult    to    interpret    from   the   numerical   analysis. 

The    problem   of    locating   occlusions    in   the   numerical    frontal 
analysis    is   manifest    over   the   eastern   North  America—western  Atlantic 
region   and    the    central-western   Pacific.      Minor   bulges    in   the   DG0*  iso- 
lines   may  be   associated  with   occluded   structure   but   a    unique   empirical 
solution   for    this   aspect    of   the   analysis    is    difficult    to   program.      The 
occlusion    problem  relates    to    the   direction    in  which   the   virtual   potential 
temperature   derivatives   are   computed    (Eqs .     [2]    and    [3]).      A   typical  view 
of   NMC's    subjectively-analyzed   surface    occlusions    in   relation   to   the 
1000/500  mb    thickness    lines    serves    to    illustrate    the    problem  well    (Fig.    7) 
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The  occluded  fronts  over  the  Pacific  Ocean  west  of  180  deg  and  over  the 

western  and  eastern  United  States  intersect  thickness  lines  (i.e.,  mean 

isotherms)  at  angles  which,  in  many  instances,  are  as  much  as  90  deg. 

Such  frontal  orientations  are  not  likely  to  be  derived  by  computations 

of  DG0*  from  Eq .  3,  as  may  be  noted  in  Fig.  6.   With  reference  to  Eqs . 

(2)  and  (3),  if  the  front -location  derivatives  could  be  calculated  in 

the  n   rather  than  in  the  n. .  direction  the  problem  associated  with  the 

t  9* 

numerical  analysis  of  occlusions  would  not  exist.   A  later  section  of 
this  chapter  discusses  a  modification  to  the  DG0*  model  which  somewhat 
alleviates  the  problem. 

The  multiplicity  and  orientation  of  fronts  implied  in  Fig.  6  has 
more  similarity  to  the  frontal-contour  analyses  performed  by  the  Canadian 
Meteorological  Service  (Anderson  et .  al.,  1955)  than  by  NMC .   Lacking  an 
absolute  standard,  analyses  from  both  weather  services  have  been  taken  as 
the  best  available,  especially  in  the  dense-data  North  American  and  North 
Atlantic  areas,  and,  as  such,  have  been  used  in  statistical  comparisons 
with  the  objective  frontal  positions.   A  later  section  of  this  chapter 
presents  results  from  the  coevaluation  of  NMC  and  objective  fronts. 

Numerically-analyzed  fronts  tend  to  be  stronger  over  land  than 
ocean;  fronts  near  the  land-ocean  interface  of  continents  (especially 
eastern  continents)  are  pronounced.   Part  of  the  reason  for  these 
analysis  features  is  due  to  the  sparsity  of  data  over  oceans  and  areas 
of  rugged  terrain.   However,  it  is  also  true  that  the  relatively  large 
intensity  over  land  is  due  to  the  rapidly-varying  height  and  nature  of 
the  underlying  surface.   Use  of  0*  vice  9*  has  not  eliminated  topographic 
effects  on  the  objectively-derived  frontal  patterns,  as  indicated  by  the 
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analysis  in  the  vicinity  of  the  Rockies,  Sierra  Madres ,  Greenland,  and 
Himalayan  areas  in  Fig.  6.   Studies  of  monthly-mean  DG0*  axes  (Chapter 
VIII)  and  attempts  to  remove  orographic  influence  and  interface  effects 
are  germane  to  this  problem  (Reed,  1961). 

As  indicated  in  Chapter  II,  limiting  the  DG0*  plot  to  5  units  and 
greater  screens  out  the  noise  frequently  associated  with  the  very  low  DG0* 
values  in  the   subbaroclinic  regions.   Such  a  scheme  also  avoids  dis- 
playing those  insignificant  zones  of  hyperbaroclinicity  where  DG0*  values 
are  <  +5  units . 

The  occurrence  of  a  double  frontal  structure  in  NMC's  analysis 
for  the  western  Pacific  but  only  a  single  HBZ  boundary  in  the  numerical 
analysis  and  the  tendency  for  DG0*  positions  to  occur  on  the  warm-side 
of  NMC  positions,  as  shown  both  in  the  west  and  east  Atlantic  (Fig.  6) 
are  manifestations  of  broadening  effects  and  the  close  proximity  of  HBZs . 
Such  problems  were  discussed  in  connection  with  the  simulation  experi- 
ments in  Chapter  II.   Grossness  of  mesh  length,  peculiarities  of  the  data 
processing,  and  the  differentiation  and  interpolation  schemes  used 
play  a  part  in  explaining  the  differences  in  positions  of  manually- 
analyzed  and  numerically-analyzed  fronts. 

Fig.  8  illustrates  FNWC's  operational  version  of  the  DG0*  analysis 
for  0000  GMT  6  March  1966.   The  operational  chart  is  a  function  of  the 
limited  amount  of  data  received  by  observation  time  plus  four  hours.   This 
fact  and  the  retention  of  a  few  small  isolated  DG0*  zones  in  the  opera- 
tional version  (areas  A  and  B)  account  for  the  minor  differences  existing 
between  Figs.  6  and  8. 

Figs.  9  and  10  show  the  operational  1000  mb  numerically-analyzed 
fronts  superimposed  on  the  sea- level  pressure  field  and  the  850  mb 
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Fig.  8.   Operational  DG0*  (1000/700  mb)  analysis  for  0000  GMT 
6  March  1966.   Isolines  at  interval  of  5  x  10~2  C  (100  km)"2 ,  starting 
with  the +5  line.   DG9*  axes  (dashed  lines)  are  manually  positioned  for 
portions  of  the  Atlantic  area,  similar  to  Fig.  6.   A  and  B  indicate 
insignificant  DG9*  features. 
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Fig.  9.   FNWC  sea- level  pressure  field,  0000  GMT  6  March_1966. 
Isobar  interval:   4  mb.   1000  mb  numerically-analyzed  fronts  (DG0*  axes) 
from  Fig.  8  (heavy  solid  lines)  with  movement  characteristics  shown  for 
the  west  and  east  Atlantic  Ocean  fronts  by  conventional  symbols. 
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Fig.  10.   FNWC  850  mb  isotherm  and  contour  field,  0000  GMT 
6  March  1966.   Isotherms  (thin  solid  lines)  at  5  C  interval,  contours 
(thin  dashed  lines)  at  60  met  interval.   1000  mb  numerically-analyzed 
fronts  (DG0*  axes)  from  Fig.  8  (heavy  solid  lines)  with  movement 
characteristics  shown  for  the  west  and  east  Atlantic  Ocean  fronts  by 
conventional  symbols. 
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isotherm-contour  chart,  respectively.   The  symbols  designating  warm, 

cold,  and  stationary  are  derived  from  geostrophic  winds  at  sea-level  and 
thermal  conditions  at  850  mb .   Most  of  the  DG0*  axes  are  located  in  or 
are  close  to  troughs  in  the  pressure  field,  giving  evidence  of  internal 
consistency  in  the  mass  structure  model  used  to  process  the  input  data 
at  FNWC   The  850  mb  isotherm  pattern  agrees  quite  well  with  the  position 
of  the  objective  fronts.   There  is  a  distinct  tendency  for  the  isotherms 
to  parallel  the  HBZ  boundary  except  near  pressure  centers,  similar  to 
manually-analyzed  fronts  (Fig.  7). 

Objective  Frontal  Prognosis  for  1200  GMT  7  March  1966 
The  numerical  model  described  here  may  be  used  equally  well  for 
objectively  locating  fronts  in  the  prognostic  thermal  fields.   Fig.  11 
shows  the  operational  36-hour  1000  mb  prognosis  of  DG0*  and  the  verifying 
chart  for  1200  GMT  7  March  1966.   The  FNWC  prognostic  model  used  to  gen- 
erate the  isotherms  from  which  DG0*  is  computed  is  a  modified  barotropic 
scheme  (Hughes,  1967).   Such  a  model  may  be  expected  to  smooth  or  elim- 
inate the  minor  thermal  waves  but  leave  the  more  intense  and  extensive 
synoptic-scale  features  intact.   It  is  for  this  reason  that  the  level 
(DG0*)  vice  the  layer  (DG0*)  HBZ  analysis  is  applied  to  the  prognostic 
thermal  fields.   It  is  emphasized  that  the  prognostic  frontal  positions 
result  from  an  analysis  of  prognostic  thermal  fields. 

Comparison  of  Subjective  and  Numerical  Frontal  Analyses 

1000  mb  Analyses 

Comparative  statistics  on  the  objectively-  and  subjectively- 
analyzed  frontal  positions  serve  to  quantify  the  observations  on  DG0* 
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made  earlier  in  the  chapter  as  well  as  give  direction  to  the  search  for 

modifications  to  the  present  program.   It  should  be  emphasized  that  an 
absolute  standard  for  locating  HBZs  does  not  exist,  therefore  it  is 
necessary  to  compare  manual  frontal  analyses  of  high  quality  with  the 
numerical  product.   The  NMC  analyses  were  chosen  for  this  purpose. 

The  most  extensive  evaluation  was  made  for  a  late  summer  period 
of  sixteen  observations  times  from  23  August  to  6  September  1965  (Table  6). 
Over  250,000  nautical  miles  of  numerically-located  1000  mb  fronts  were 
compared  with  the  analogous  NMC  surface  fronts.   The  more  obvious  con- 
clusions from  this  study  are  as  follows: 

1.  The  two  independent  frontal  analyses  have  the  closest  relation 
over  dense-data  land  and  ocean  areas  ,  such  as  the  United  States  and  the 
Atlantic  Ocean  while  the  relation  is  poorest  over  the  Pacific  Ocean  and 
North  African  areas.   The  1.2  deg  lat  difference  for  Asia  is  considered 
anomalously  low  and  not  representative  of  that  area. 

2.  The  objective  HBZ  positions  north  of  30N,  and  over  eastern 
portions  of  continents  and  western  portions  of  oceans,  compare  most 
favorably  with  the  NMC  positions. 

3.  Only  66  percent  of  the  NMC  fronts  could  be  associated  with 
the  objective  fronts.   The  remainder  are  occlusions  and/or  weak  fronts 
in  eastern  ocean  areas.   Both  of  these  types  are  frequently  analyzed 
from  time  continuity  considerations  only  and  are  in  the  process  of 
frontolyzing. 

4.  A  similar  percentage  of  significant  objectively-analyzed  fronts 
(72  percent)  could  be  associated  with  NMC  fronts.   The  remainder  represent 
baroclinic  zones  over  North  Africa  and  other  subtropical  ocean  and  land 
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areas    as   well   as    those    fixed    to   orographic    features.      Following    preva- 
lent  air   mass    theory,    the    subjective   analyst   hesitates    to    place   stationary 
fronts    on    the    periphery   of   desert    or  mountainous    areas,   despite    the    nature 
of   thermal   patterns. 

A  more    recent,    but   quite    limited,    evaluation   for    the    four   synoptic 
times,    0000   GMT  and    1200   GMT,    26   and   27    January    1967,    considered  approxi- 
mately  57,000   nautical  miles    of  Northern  Hemisphere    fronts  without   regard 
to   classification   by   region.      An  average    position  difference    of   2.6   deg 
lat   between   NMC   surface    and    objective    (DG9*  model)    1000  mb    fronts    resulted, 
Forty-two    percent    of   the    numerically-computed    fronts   were    found   on   the 
warm-air   side    of   NMC   fronts   with   an   equal    number   on   the   cold-air   side. 
The    remainder  were    coincident,    to  within    1/4   deg    lat. 


Tab_le    6.      Statistical   relations    between   numerically-analyzed 
1000   mb    (DG9*  model)    and  manually-analyzed   NMC   surface    fronts    for   sixteen 
observation   times    in   the    period   23   August    to    6   September    1965. 


Sector  Objective   Fronts    Sampled  Average   Distance 

(n  mi)  Separating  NMC  & 


Objective   Fronts 
(deg   lat) 


North  America  49,250  1.6 

Atlantic  ,                               46,900  1.4 

Pacific  46,300  2.1 

Europe-North   Africa  26,000  1.8 

Asia  18,000  1.2 

Overall   average  1.7 
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850  mb  Analyses 

Although,  at  this  point,  it  is  not  the  intent  to  discuss  verti- 
cal consistency  of  numerically-analyzed  fronts  (see  Chapter  V  )  it  is 
appropriate  to  the  substance  of  this  chapter  to  show  the  850  mb 
DG0*  version  for  0000  GMT  6  March  1966,  in  order  to  demonstrate  the 
relative  merits  of  numerical  frontal  analysis  using  mandatory  level  and 
layer  data . 

From  the  conventional  synoptic  point  of  view,  the  850  mb  fronts 
are  expected  to  be  on  the  cold-air  side  of  their  1000  mb  counterparts. 
But,  as  noted  from  a  comparison  of  Figs.  6  and  12,  the  objective  850  mb 
fronts  are  coincident  with  or  are  on  the  warm-air  side  of  the  1000  mb 
fronts  derived  from  the  DG0*  model.   To  show  that  this  abnormal  situ- 
ation is  not  uncommon,  the  following  is  offered  in  evidence.   With 
reference  to  the  evaluation  of  late  January  1967  frontal  analyses, 
described  above,  the  average  distance  between  850  mb  NMC  and  objectively- 
located  fronts  is  2.6  deg  lat  which  is  the  same  distance  as  that  for  the 
NMC  surface  versus  objective  1000  mb .   However,  only  19,000  nautical 
miles  of  front  were  considered  at  850  mb ,  most  of  which  were  over  the 
dense-data  North  American  region.   Thus,  considering  the  area  sampled, 
the  850  mb  positions  of  numerically-analyzed  and  NMC  fronts  compare  less 
favorably  than  their  counterparts  at  or  near  the  surface.   But,  of  greater 
importance  is  the  fact  that  86  percent  of  the  850  mb  numerically- located 
positions  were  found  on  the  warm  side  of  NMC's  850  mb  fronts,  with  an  alge- 
braic average  distance  of  2.3  deg  lat.   Such  bias  relative  to  NMC  did  not 
exist  for  the  objective  1000  mb  fronts  derived  from  9*. 
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Fig.  12.   850  mb  DG0*  analysis  for  0000  GMT  6  March  1966, 
Isolines  at  interval  of  20  x  10-2  C  (100  km)"2,  starting  with  +5 
units.   The  available  NMC  fronts  at  850  mb  are  shown  with  conven- 
tional symbols. 
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Modification  of  DGO*  to  Enhance  the  Analysis  of  Occlusions 
One  of  the  more  obvious  omissions  in  the  analyses  generated  by 
the  DG9*  model  is  the  conventional  occluded  structure,  as  discussed 
earlier  in  the  chapter.   Occlusions  are  accepted  and  common  to  all 
national  and  international  weather  services,  except  perhaps  that  of  the 
Canadians,  whose  working  model  features  troughs  of  warm  air  aloft 
(trowals)  as  substitutes  for  the  occlusion  (Galloway,  1958).   The  well- 
developed  occlusions  on  the  surface  charts  of  United  States  weather 
services  are  generally  associated  with  a  large-amplitude  prefrontal 
tongue  of  warm  air,  as  may  be  noted  over  the  western  Atlantic—eastern 
North  America  region  and  over  the  western  North  Pacific  Ocean  at  0000 

GMT  6  March  1966  (Fig.  7).   Such  occluded  fronts  intersect  the  isotherms 

— » 
at  considerable  angles  and  nQ.  becomes  a  relatively  poor  approximation  to 

n  .   (See  Eqs .  2  and  3). 

The  close  coincidence  of  the  thermal  axis  of  warm  temperature  and 
the  conventional  occlusion  suggests  application  of  the  Laplacian  operator 
to  temperature  in  order  to  locate  the  "occluded"  HBZ  boundaries.   Using  9* 
as  the  thermal  field  of  interest, 

_V  9*  =  -  V|V9*|  -n-   -  Iv9*|v-n-  s  DG9*  +  GKO* .      (11) 
As  noted  in  Eqs.  (5)  and  (6),  the  DG0*  parameter  is  indicative  of  the 
vorticity  of  the  thermal  wind  due  to  shear  and  thus  GK0*  is  representa- 
tive of  the  curvature  component  of  the  thermal  vorticity. 

Fig.  6  demonstrates  the  poor  relation  of  manually  drawn  occlusions 
to  the  field  of  DG9*.   GK0*  (Fig.  13),  which  isolates  the  effect  of  the 
curvature  term,  is  not  adequate  either  to  depict  occlusions  or  any  other 
frontal  phenomena  for  that  matter.   The  same  comment  is  true  for  the 
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Fig.  13.   GK9*  (1000/700  mb)  analysis  for  0000  GMT  6  March 
1966.   Isolines  at  intervals  of  20  x  10"2  C  (100  km)"2,  starting 
with  +5  units.   NMC  fronts  are  shown  with  conventional  symbols. 
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Laplacian  itself  (not  shown)  since  thermal  tongues  also  occur  in  areas 
other  than  the  location  of  conventional-front  discontinuities.   However, 
it  is  productive  to  combine  the  curvature  and  shear  terms  of  Eq .  (11)  in  a 
selective  manner,  that  is  only  in  the  area  of  negative  DG6*  between  the 
zero  values  within  the  HBZ  and  the  adjacent  subbaroclinic  zone.   This 
version  is  called  Mod  DG0*  and  is  shown  in  Fig.  14.   Such  a  procedure 
preserves  the  continuity  of  the  non-occluded  warm-side  boundaries  of  the 
HBZs  depicted  by  DG0*  alone  since  these  occur  with  positive  DG0*.   A  blow- 
up of  a  section  of  the  numerical  frontal  analysis  over  the  western  Atlantic- 
eastern  United  States  area  in  Fig.  14  exemplifies  the  enhancement  of  the 
occlusion  analysis  by  use  of  Mod  DG0*  (Fig.  15).   Note  that  DG0*  and  Mod 
DG0*  are  the  same  except  on  the  cold-air  side  of  zero  DG0*  in  the  HBZ. 

It  is  of  interest  to  compare  Figs.  6  and  14.   Note  the  differences 
in  orientation  and  structure  of  the  HBZ  boundaries  that  are  associated 
with  occlusions,  as  for  example  over  the  extreme  eastern  and  western 
sections  of  northern  United  States  as  well  as  northern  Russia.   Unfortun- 
ately, the  curvature  modification  afforded  by  Mod  DG0*  is  not  always 
desirable.   An  example  of  this  may  be  noted  over  the  subtropical  area 
of  North  America,  a  region  frequently  associated  with  low-tropospheric 
thermal  tongues.   Here,  an  already  troublesome  area  for  frontal  analysis, 
due  to  topography  and  sparse-data  problems,  suffers  further  from  the 
curvature  modification. 

Overall,  the  major  objection  to  the  Laplacian  adjustment  for 
occluded  structure  is  that  the  curvature-vorticity  term  is  not  directly 
related  to  baroclinicity  and  hence  fictitious  temperature  gradients  are 
implied.   These  objections  have  kept  the  program  from  becoming  operationally 
useful. 


63 


_    Fig._14.   (DG0*  +  GK0*)  for  DG0*  <  0,  GK0*  >  0,  and 
|gK0*|  >  |DG9*| :   0000  GMT  6  March  1966.   Isolines  at  intervals 
of  20  x  10"2  C  (100  km)"2  starting  with  +5  units.   Axis  values  > 
IS  units:   xxxx :  NMC  fronts  shown  with  conventional  symbols. 
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Fig.  15.   Modification  to  the  DG0*  (1000/700  mb)  analysis  for 
the  case  of  a  frontal  occlusion  over  the  western  Atlantic—eastern 
United  States  area,   0000  GMT  6  March  1966.   9*  j.n  K:   light  solid 
lines;  major  axis  of  zero  DG9*:   dotted  line;  DG0*  and  mod  DG0*  in 
10~2  C  (100  km)"* :   heavy  solid  and  dashed  lines,  respectively. 
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Some  work  has  been  done  by  Kirk  (1965,  1966,  and  1967)  on  incor- 
porating the  Laplacian  effect  into  objective  frontal  analyses.   The 

2 

unqualified  favorable  appraisal  of  V  T  as  a  front-location  parameter, 

given  by  Kirk  in  these  references,  completely  misrepresents  the  utility 
of  this  approach. 

Effect  of  Mesh- length  Variation  on  DG9* 
In  Chapter  II  the  effect  of  reduction  in  mesh  length  on  the 
location  and  intensity  of  frontal  zones  was  demonstrated  for  a  limited 
case  involving  the  modeled  N  and  M  temperature  profiles.   The  results 
there  suggest  a  similar  study  using  real-atmosphere  data. 

Consistent  with  both  the  simulation  and  operational  experiments, 
the  mean  virtual  potential  temperature  for  the  layer  1000/700  mb  (9*) 
and  the  synoptic  time,  0000  GMT  6  March  1966,  were  selected  for  this 
substudy  relating  mesh  length  variation  to  objective  frontal  analysis. 
In  particular,  the  frontal  zone  chosen  for  study  is  that  portion  of  the 
significant  HBZ  in  the  western  At lantic--eastern  North  American  region 
contained  in  the  area  identified  by  corners  ABCD  in  Fig.  6.   The  subject 
HBZ  is  considered  as  moderate  (M)  in  terms  of  the  symbolic  profiles  in 
Chapter  II  (Fig.  3);  the  average  baroclinic ity  is  1 . 6C  (100  km)    through 
the  frontal  zone.   Starting  with  the  full  mesh  grid-point  values  of  9*, 
the  temperature  data  were  numerically  interpolated  to  the  .75,  .5,  and 
.25  mesh  intersections  using  Bessel's  central  difference  formula.   Sub- 
sequently, the  fields  of  G0*  and  DG0*  were  produced  for  each  of  the  new 
grids  in  a  manner  consistent  with  that  described  earlier  in  the  Chapter. 

Fig.  16  illustrates  the  important  aspects  of  the  DG0*  analysis 
on  the  full,  half,  and  quarter  mesh  grids.   Results  for  the  three-quarter 
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Fig.  16.   DG9*  analysis  as  a  function  of  mesh  length   for 
area  identified  by  corners  ABCD  _in  Fig.  6,  0000  GMT  6  March  1966. 
Positive  and  negative  axes  of  DG0*:   from  full  mesh  (heavy  solid 
lines),  half  mesh  (heavy^  long  dashed  lines),  and  quarter  mesh  (heavy 
short  dashed  lines).   DGG*  isolines  intersecting  positive  (negative) 
axes,  starting  with  +.35  (-.25)  x  10~2  C  (100  km)"^  at  intervals  of 
10  units:   from  full  mesh  (light  solid  lines),  half  mesh  (light  dashed 
lines),  and  quarter  mesh  (dotted  lines).   9*  isolines  (K)  indicated  by 
edges  of  shaded  areas. 
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mesh  were  intermediate  to  those  for  the  full  and  half  mesh  and  are  not 

shown  here.   A  background  of  full-grid  9*  (shaded)  is  shown  to  point  out 
an  important  departure  from  the  modeled  profiles  in  Chapter  II,  namely, 
that  the  warm  air  is  much  less  barotropic  and  the  cold  air  somewhat  more 
baroclinic  than  the  counterpart  sections  in  the  M  profile  (Fig.  3  and 
Table  1).   However,  at  J  <  15  similar  conditions  in  the  modeled  and  real- 
atmosphere  temperatures  do  prevail. 

Analysis  of  the  positions  of  warm-  and  cold-side  boundaries  of 
the  HBZs  and  the  frontal-zone  intensities,  as  implied  by  the  numerical 
DG0*  analysis,  indicates  the  following: 

1.  Full  versus  half  mesh:   Overall,  both  full  and  half  mesh 
patterns  are  similar  in  position  and  extent  of  the  HBZ .   Clear  evidence 
of  the  positive  correlation  of  mesh  size  and  frontal-zone  width  is  con- 
fined to  the  area  below  J14  where  the  latter  is  reduced  up  to  10  percent 
(40  kms)  in  the  half- compared  to  full-mesh  case.   Poleward  of  J14  the 
half-mesh  HBZ  widths  show  both  increases  and  decreases  compared  to  full 
mesh,  varying  in  amount  from  near  zero  to  25  percent.   A  large  segment  of 
the  northern  section  shows  both  the  warm-  and  cold-side  HBZ  boundaries 
displaced  toward  the  warm  air  in  the  half-mesh  compared  to  the  full-mesh 
analysis . 

The  reduction  from  full-  to  half-mesh  relates  well  to  the  increase 
of  intensity  values,  both  along  the  warm-  and  cold-side  boundaries  of  the 
HBZ,  as  expected  from  Chapter  II. 

2.  Full  or  half  mesh  versus  quarter  mesh:   Fracturing  of  HBZ 
boundaries  is  very  evident.   In  addition,  the  HBZ  intensity  does  not 
uniformly  increase  from  the  .5  to  .25  grid.   Apparently,  grid-size 
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reduction  much  below  .5  is  inconsistent  with  portraying  synoptic -scale 
fronts  from  data  originally  fitted  to  the  380  km  full-mesh  grid.   With 
an  initial  gross  temperature  resolution,  the  reduction  to  .25  of  the 
full  mesh  (95  km)  highlights  mesoscale- like  variations  in  the  derivative 
fields.   Such  detail,  partly  derived  from  numerical-finite  differencing 
and  interpolation  procedures,  may  be  considered  noise  here  since  the  data 
at  full  mesh,  as  a  function  of  the  mass  structure  model,  do  not  contain 
such  thermal  resolution. 

It  should  be  noted  that  some  error  in  the  relative  positioning 
and  magnitude  of  the  numerically-analyzed  fronts  in  Fig.  16  is  due  to 
the  difficulty  of  locating  the  axes  with  an  accuracy  excelling  .2  mesh 
length,  even  with  the  reasonably  small  analysis  interval  of  10  units  of 
DOG*. 

Summary 

In  summary,  the  results  using  real-atmosphere  data  are  in  general, 
but  not  specific,  agreement  with  the  simulated  data  experiment  in  Chapter 
II.   It  is  evident  that  the  resolution  of  the  temperature  data  input  to 
the  DG0*  program  must  be  enhanced  before  reductions  in  grid-size  much 
below  380  km  can  be  effective  in  specifying  the  true  width  and  intensity 
of  frontal  zones.   Numerically- processed  and  analyzed  thermal  data,  of 
the  type  utilized  here,  presently  do  not  contain  the  detail  inherent  in 
data  extracted  from  high  quality  subjective  analyses. 

In  order  to  focus  attention  on  the  problem  of  appropriate  scale 
and  resolution  for  acceptable  numerical  frontal  analysis,  the  next 
chapter  considers  procedures  for  increasing  the  detail  of  the  baroclinicity 
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data  used  in  the  DG0*  computation.   Highly  resolved  data,  containing 
information  on  the  mesoscale  structure  of  the  atmosphere, is  readily 
extractable  from  vertical  soundings  of  the  atmosphere. 


CHAPTER  IV 
EXPERIMENTS  TO  ENHANCE  THE  RESOLUTION  OF  BAROCLINICITY 

Most  of  the  deficiencies  of  the  operational  DG0*  model  are  due  to 
the  poor  resolution  of  the  temperature  and  baroclinicity  fields  from  which 
DG0*  is  computed.   Frontal  phenomenon,  by  its  nature,  tends  toward  a  meso- 
scale  structure  while  numerica J ly-processed  observations  describe,  at 
best,  the  synoptic -scale  features.   Therefore,  some  experimental  effort 
was  devoted  to  the  problem  of  enhancing  the  resolution  of  HBZ  analysis  to 
a  degree  commensurate  with  the  entity  being  analyzed. 

Modifications  to  DG9*,  such  as  Mod  DG9*  in  the  last  chapter,  may 
improve  the  meteorological  consistency  of  the  present  operational  frontal 
model,  but  its  merits  cannot  exceed  the  amount  or  scale  of  the  information 
content  inherent  in  the  data  tuned  to  a  mesh  size  of  381  km.   Thus,  it  is 
logical  to  look  to  the  highly  resolved  significant-level  radiosonde  data 
for  improving  the  HBZ  analysis.   One  of  the  more  promising  experiments 
using  detailed  sounding  data  has  its  foundation  in  the  analytical  develop- 
ment which  follows : 

6  =  f  (nQ,p),  (12) 

where  9  is  potential  temperature;  p  is  pressure;  and  n_,  defined  by 

— *  i         i 

nfl  =  V9/|V9| ,    represents   an   axis    oriented    in  a   direction   perpendicular 
to   an   isentrope    on  an   isobaric   surface.      Therefore, 

dG   "   an~6n9+   c^6p'  (13) 

y 
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Along  an  isentropic  surface,  d0  =  0,  and  from  Eq .  13 

Ges^  =  .ae^. 

dng     dp  6n0 
Similar  to  Eq.  (3) , 

a2e 

DGe  s  -  VG0-nQ  =  -  ^-^  .  (15) 

an0 

Computation  of  G0  and  DG0 
The  scheme  proposed  here  involves  a  unique  computation  of  G0  from 
the  two  factors,  -d0/dp  and  6p/5n  This  approach  represents  a  consider- 
able deviation  from  that  used  to  obtain  G0*  as  a  function  of  0*,  since  G0 
is  not  derived  directly  from  a  field  of  0.  However,  DG0  is  computed  from 
G0,  similar  to  DG0*  from  G0*. 

The  Computation 
of  d0/dp 

The  derivative,  d0/dp,  is  determined  for  a  mandatory  pressure  sur- 
face by  a  quadratic  fit  of  potential  temperature  data  at  the  mandatory 
level  and  the  first  significant  levels  above  and  below.   Plotted  soundings 
on  a  thermodynamic  chart  give  the  impression  of  a  sequence  of  layers,  each 
linear  in  the  logarithmic  pressure  ordinate,  with  separating  first-order 
thermal  discontinuities.   In  reality,  the  discontinuities  bounding  these 
significant  layers  are  manufactured  by  criteria  used  for  processing, 
coding,  and  plotting  the  observed  data  (Danielsen,  1959).   Therefore, 
fitting  a  continuous  curve  to  the  transmitted  sounding  for  the  purpose 
of  obtaining  d0/dp  is  reasonable  as  well  as  necessary.   Also,  using  a 
quadratic  rather  than  a  linear  fit  of  temperature  data  agrees  with  the 
nature  of  the  sounding  in  pressure  and  temperature  coordinates  and  ade- 
quately handles  the  situations  in  which  the  mandatory  level  is  also  a 
significant  level. 


73 
The   mandatory- level   values    of   S9/Sp  may  be    obtained   through   the 

solution   for   constants    a,    b,    and   c    in    the    system  of   three    simultaneous 

linear   equations 

2 

0     =a  +  bp     +cp 
u  u  u 

9     =  a  +  b   p     +  cp  2  (16) 

m  mm 

9t    =  a  +  b   Pl  +  cpt    , 

and   from  f-~r—  J     =  b  +  2cp    ,  (17) 

^m 

where  u(l)  is  the  significant  level  geometrically  above  (below)  the  man- 
datory pressure  level,  m.  The  9s  in  Eq .  (16)  are  computed  from  sounding 
values    of   pressure   and   temperature    by         R 

|  -  {^f)>  (18) 

In  order  to  keep  the  program  as  simple  as  possible  in  the  initial 
feasibility  experiment  for  computing  G9,  temperatures  to  the  nearest  whole 
deg  C  (accurate  to  .5C)  were  employed  in  computing  cS9/dp.   This  approach  was 
partly  dictated  by  the  state  of  the  hemispheric  preprocessed  radiosonde 
data  made  available  by  FNWC  . 

Significant  layers  of  shallow  depth  (i.e.,  less  than  10  mb)  were 
not  considered  at  this  point  in  the  experiment  partly  due  to  the  gross 
temperature  resolution  and  also  because  such  layers  may  contribute  intol- 
erable error  in  the  bQ/bp   computation.   Further,  as  an  example  of  the 
inconsequence  of  a  thin  layer,  consider  a  stable  significant  layer  near 


Since  the  data  input  to  FNWC's  data-processing  system  excludes  signifi- 
cant level  information,  it  was  necessary  to  write  a  complex  program  to 
couple  significant  and  mandatory  level  data  at  the  preprocessed  stage. 
The  program  was  written  by  Lieutenant  R.  C.  Schiffner,  Department  of 
Meteorology  and  Oceanography,  Naval  Postgraduate  School,  Monterey, 
California . 
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700  mb  of  10  mb  thickness  (approximately  400  feet)  and  100  miles  in 
horizontal  width.   The  slope  of  this  layer,  6z/6n,  is  1/1320.   Hori- 
zontally wider  and  vertically  thinner  stable  layers  result  in  still 
smaller  slopes  and  concomitant  negligible  baroc linicities  along  the 
pressure  surface.   Thus,  HBZs  of  such  dimensions  may  be  profitably 
omitted  from  the  frontal  analysis,  at  least  in  this  first  attempt  at 
computing  G0. 


The  Computation 
of  6p/6n 

Unlike  c39/dp,  the  computation  of  the  slope  of  the  isentropic 

surface,  6p/6nQ,  at  each  station  was  accomplished  within  the  framework 

of  FNWC's  data  processing  and  analysis  system.   Thus,  9*  was  substituted 

for  9  in  Eq.  0-4).  Available  to  the  project  were  FNWC's  numerically-computed 

fields  of  pressure  for  each  9*  surface  at  2 . 5C  intervals.   From  these 

fields,  slopes  of  the  isentropic  surfaces,  |(6p/6n)  n  |=  6p/6n  ,  were  com- 

P        P 
— » 
puted  by  manual  means.   Here,  n   is  the  unit  vector  normal  to  isobars  on 

P 

l        — *   i 
the  isentropic  surface.   However ,  | (6p/6n)  n  .  |  =  6p/6n   , the  form  of  the  slope 

factor  in  Eq.  (14),  represents  only  that  component  of  6p/6n  which  lies  in 

the  direction  of  V9*  on  the  isobaric  surface.   Since  the  angles  between  the 

— *  — ¥ 

vectors  n  and  n   were  found  to  be  <  15  deg  for  nearly  all  of  the  radio- 
sonde data  used,  the  difference  in  the  two  quantities  is  negligible  (i.e., 
generally  <  3.5  percent  in  this  particular  experiment).   In  fact,  over- 
estimating 6p/6n   by  using  6p/6n   in  Eq.  (14)  is  considered  desirable  in 
view  of  the  smoothing  inherent  in  the  numerically-processed  products. 
Moreover,  6p/6nQ.  was  computed  for  that  isentropic  surface  whose  magnitude 
is  divisible  by  2.5  and  greater  than  9*  at  the  mandatory  level  for  which  the 
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computations  are  being  made.   For  example,  6p/6n   ,  for  a  station  with 
0*  =  301K  at  850  mb ,  was  computed  at  the  0*  surface  of  302. 5K.   Undoubt- 
edly some  fiction  was  introduced  into  the  computed  G©  by  this  simplified 
and  hybrid  combination  of  cSO/Sp  and  6p/6n   .   However,  in  the  realm  of  a 
feasibility  study  such  lack  of  homogeneity  and  detail  may  be  tolerated, 
and,  overall,  may  be  desirable  as  will  be  discussed  in  a  later  section. 

Although  the  ultimate  goal  is  to  improve  upon  the  operational 
DG0*  model  for  frontal  analyses,  the  emphasis  here  is  on  its  predecessor, 
the  baroclinicity ,  whose  properties  largely  determine  the  validity  and 
usefulness  of  the  HBZ  analysis.   Thus,  the  discussion  at  hand  will  con- 
centrate on  the  enhancement  of  the  field  of  baroclinicity  as  generated 
by  the  introduction  of  the  vertical  resolution  inherent  in  significant- 
level  radiosonde  data. 

The  gross  scheme  just  outlined  for  computing  G9  and  its  attendant 
field  of  DG0  will  be  symbolized  as  G0'  and  DG0'  to  distinguish  it  from 
subsequent  experiments  which  more  closely  approximate  the  true  resolution 
of  G9  and  DG0.   The  test  date  is  the  same  one  used  to  exemplify  the  G0* 
and  DG0*  computations,  namely  0000  GMT  6  March  1966. 

G91  and  its 
Comparison  to  G9* 

In  order  to  establish  continuity  with  the  discussion  in  Chapter 
III  and  present  a  base  to  which  G01  may  be  compared,  the  numerically- 
analyzed  850  mb  G0*  for  0000  GMT  6  March  1966  (Fig.  17a)  is  discussed 
first.  The  analysis  area  is  confined  to  North  America  and  the  middle 
and  high  latitude  portions  of  the  adjacent  ocean  areas  from  Europe  to 
180  deg  long.   The  positive  DG0*  axes  closely  parallel  baroclinicity 
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magnitudes  of  1C  (100  km)   ,  being  on  the  weaker-  (stronger-)  gradient 

side  for  the  less  (more)  intense  HBZs .   An  example  of  the  former  (latter) 

may  be  noted  in  the  middle  (western)  North  Atlantic  area.   Baroclinic 

regions  in  which  9*  gradients  are  less  than  1C  (100  km)    are  of  little 

interest  but  such  gradients  do  occur  along  the  DG0*  axes.   This  is 

especially  true  for  low- latitude  ocean  areas.   Note  the  North  Atlantic 

region  south  of  30N.   The  halfwidth  of  HBZ  zones  may  be  noted  on  Fig.  17a 

by  the  transverse  distance  from  the  DG0*  to  G0*  axes.   Average  magnitude 

of  the  halfwidth  is  2.5  deg  lat . 

Fig.  17b  shows  the  numerically  computed  G0'  field,  manually- 
analyzed  in  this  version.   About  40  percent  of  the  North  American  area 
shows  G0'  less  than  1C  (100  km)    in  the  same  area.   Seven  radiosonde 
stations  yield  G0'  values  greater  than  3C  (100  km)   ,  with  an  absolute 
magnitude  of  8. 1C  (100  km)    in  southwest  Canada,  while  only  two  G0* 
grid-point  values  exceed  three  units  of  baroclinicity . 

Table  7  statistically  indicates  the  increased  focus  of  baroclini- 
city using  the  significant- leve 1  radiosonde  data  from  118  stations  in  the 
area  of  analysis  on  Fig.  17b.   The  table  shows  the  number  of  stations  in 
each  of  nine  categories  of  AG0  =  G0'  -  G0*.   The  category  "-.19  to  +.19" 
may  be  considered  as  AG©  ~  0.   HBZs  are  arbitrarily  defined  as  the  areas 
having  G0*  values  >  1C  (100  km)   .   For  such  areas,  G0'  >  G0*  in  the  vast 
majority  of  cases  (85  percent),  while  in  sub  baroclinic  regions  (SBZ) , 
where  G0 '  <  1C  (100  km)"  ,  G0'  <  G0*  almost  without  exception  (98  percent). 
This  attests  to  the  greater  range  of  G0'  relative  to  G0*,  and  especially 
to  the  reduction  of  baroclinicity  in  the  s ubbaroclinic  areas  on  the  G0' 
chart.   Stated  another  way,  the  relative  barotropic  areas  are  more  barotropic 
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while  the  relative  baroclinic  areas  are  more  baroclinic  in  the  GO'  as 
compared  to  the  G0*  fields. 


Table  7.   850  mb  AG©  (=  G9*  -  G©*) ,  0000  GMT  6  March  1966, 
summarized  for  118  radiosonde  stations  (North  America  and  adjacent  ocean 
areas).   HBZ  =  stations  where  G©1  >  1C  (100  km)"1;  SBZ  =  stations  where 
G0'<  1C  (100  km)"1.   Stations  with  absolute  instability  omitted. 


AG©  Number  of  radiosonde  stations 
(C  [100  km]"1) 

HBZ  SBZ 

>  +  2.00  4 

+1.00  to  +1.99  5 

+.60  to  +.99  5  1 

+.20  to  +.59  14  0 

-.19  to  +.19  13  10 

-.20  to  -.59  4  20 

-.60  to  -.99  1  13 

-1.00  to  -1.99  0  26 

<  -  2.00  0  2 


The  major  axes  have  similar  orientation  on  both  the  G0*  and  G0' 
charts,  but  those  on  the  latter  are  generally  displaced  to  the  cold  side 
of  their  counterparts  on  the  former.   Note  the  situation  over  eastern 
United  States  and  the  adjacent  Atlantic  and  the  zone  running  northwest- 
southeast  from  Alaska  to  central  United  States.   These  relations  between 
G0'  and  G0*  (Fig.  12  and  Table  7),  suggest  that  the  nature  of  the  G0' 
computations  reduces  the  broadening  effect  as  compared  to  G0*.   However, 
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an  extra  impetus  to  this  effect  may  be  due  to  the  manner  in  which  the 

slopes  of  isentropic  surfaces  (6p/6n   )  were  computed.   That  is,  using 

the  "higher  9*  divisible  by  2 . 5K"  for  computing  6p/5n   has  the  effect 

of  moving  the  hyperbaroc linic  zones  toward  the  cold  air  from  their  true 

positions . 

It  is  also  noted  that  detail  along  the  G0'  axes  (eastern  United 
States  and  southwest  Canada)  is  greater  than  that  of  G0*.   Secondary  HBZs 
are  both  more  evident  and  intense  on  the  G0 '  analysis  as  seen  on  Fig.  17b 
over  east  central  United  States  and  southeast  Canada.   The  mountainous 
regions  are  problem  areas  on  both  the  G0*  and  GQ'  models;  the  two  analyses 
differ  markedly  here.   However,  occluded  structure,  as  shown  by  NMC  over 
northwest  United  States  and  southwest  Canada  (Fig.  12),  is  indicated 
weakly  with  reasonable  orientation  in  the  G0 '  but  not  in  the  G0*  fields. 

There  still  remain  many  problems  in  obtaining  a  suitable  field  of 
baroclinicity  at  any  level,  especially  at  low  tropospheric  levels.   In  the 
present  experiment,  the  G0'  computation  becomes  invalid  for  the  dry  adia- 
batic  case,  that  is  for  d9/c)p  =  0  and  6p/6nQ^  =  OT.   Also,  the  present 
mode  of  computing  G0'  is  unrealistic  for  the  occasional  case  of  super- 
adiabatic  conditions  (c)9/dp  positive)  in  the  observed  sounding  since  this 
stability  condition  cannot  exist  in  the  numerically- processed  data  from 
which  the  6p/6nQ.  values  were  generated.   Four  such  superadiabatic  values 
occurred  over  the  North  American  area  for  the  synoptic  time  under  experi- 
ment.  This  frequency  would  be  greater,  except  that  for  such  anomalous 
cases  temperatures  to  the  nearest  . 2C  (vice  1C)  were  used  in  computing 
dO/dp.   Teletype  data  transmissions  were  the  source  of  such  detailed 
temperature  data. 
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DG0'  and  its 
Comparison  to  DG0* 

To  advance  the  experiment  to  the  point  of  obtaining  a  DG0'  field 
for  comparison  to  its  counterpart,  DG0*,  the  G0'  data  were  objectively 
analyzed.   Several  G0'  guess  fields  were  considered.   Even  over  a  dense 
data  area  like  the  United  States,  a  "flat"  guess,  as  G0'  =  0,  yields  a 
numerical  G0 '  analysis  with  a  tendency  for  circular  isolines  due  to 
localized  points  of  large  baroc linicity .   Some  improvement  in  shaping  the 
field  to  elliptical  rather  than  circular  isolines  is  achieved  by  using 
the  well-behaved  G0*  field  as  first  guess  to  the  G0'  analysis.   Using 
this  method,  Fig.  18  shows  a  manual  trace  of  the  resulting  850  mb  ob- 
jectively-analyzed G0'  field.   The  reduction  of  resolution  in  G0' -objective 
compared  to  G0' -subjective  (Fig.  17b)  is  obvious.   Also,  differences  in 
the  positions  of  the  axes  of  the  two  fields  give  a  further  measure  of  the 
role  of  objective  analysis  in  misplacing  the  associated  HBZ  boundaries, 
a  subject  which  came  under  discussion  in  Chapter  II. 

The  quantity  DG0'  (=  -VG0'-n   )  was  computed  next.   The  0*  field 
was  used  to  generate  the  directional  second  derivative  since  a  suitable 
numerical  analysis  of  9'  was  not  readily  available. 

Fig.  19  shows  the  DG0*  and  DG0 '  numerical  analyses,  similar  to 
Fig.  17  for  G0*  and  G0'.   In  general,  the  comments  concerning  the  baro- 
clinicity  fields  carry  over  to  the  interpretation  of  Fig.  19.   Due  to  the 
close  proximity  of  G0'  baroclinic  zones  over  the  southeastern  United 
States,  the  nearly  meridional  DG0'  axis  along  the  east  coast  is  disrupted. 
This  is  an  undesirable  feature  in  view  of  the  conventional-front  positions 
for  this  time. 
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Also  note  that  the  axis  of  maximum  baroc linicity  in  the  G0'  and 
G0*  fields  do  not  necessarily  lie  along  their  zero-isoline  counterparts 
in  the  DG0'  and  DG0*  patterns,  respectively,  as  expected  from  the  dis- 
cussion in  Chapter  III.   The  phenomenon  is  of  significance  only  for  com- 
plex HBZ  structures  which  occur  near  the  edges  of  dense  data  areas  or  in 
mountainous  regions.   The  detailed  G0' ,  DG0'  analyses  are  more  susceptible 
than  GQ*,    DG0*  to  this  numerical  problem,  as  noted  over  Mexico  and  the 
adjacent  Gulf  of  Mexico  (Fig.  19). 

Fig.  20  gives  a  measure  of  the  positions  of  the  HBZ  boundaries 
from  the  DG0'  and  DG0*  models  relative  to  subjectively-analyzed  850  mb 
NMC  and  Canadian  fronts.   DG0'  excels  in  detail  and  in  apparently  reducing 
HBZ  broadening  by  placing  fronts  on  the  cold-air  side  of  the  DG0*  positions 
However,  the  DG0*  fronts  compare  very  favorably  with  the  hand-analyzed 
versions  in  the  matter  of  spatial  continuity,  especially  over  the  western 
Atlantic  Ocean. 

700  mb  data  and  fields  were  handled  in  the  same  way  as  850  mb  and 
the  results  show  similar  merits  and  deficiencies.   The  vertical  consist- 
ency of  objective  fronts  from  DG0'  excels  that  from  DG0*,  as  noted  in 
Fig.  21.   Improper  frontal  slopes  (i.e.,  700  mb  positions  on  the  warm-air 
side  of  850  mb  positions)  are  more  prevalent  in  the  DG0*  than  in  the  DG0' 
patterns.   The  vertical  consistency  of  DG0*  fronts  is  illustrated  and  dis- 
cussed further  in  Chapter  V. 

Before  advancing  the  G0  experiment  (Eqs .  12-18)  in  order  to  obtain 
a  greater  measure  of  accuracy  and  resolution  in  the  baroclinicity  fields, 
it  is  to  be  emphasized  that  the  difference  in  detail  between  the  G01  and 
G0*  fields  is  derived  from  the  computation  of  cl9/dp  from  signif icant-leve  1 
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radiosonde  date  (Eq.  14).   The  influence  of  the  9*  field  is  retained  by 
substituting  6p/6n  .  for  6p/6n   in  Eq.  (14)  and  n    for  n   in  Eq .  (15). 

Recomputation  of  G9  with  Certain  Approximations  Removed 
The  previous  section  may  be  regarded  as  a  successful  feasibility 
experiment  in  obtaining  a  field  of  baroclinicity  whose  detailed  structure 
excels  that  of  G0*  due  to  the  vertical  resolution  involved  in  the  com- 
putation of  d9/dp.   Utility  of  the  derived  G0'  and  DG0'  is  evident  even 
though  (1)  error-prone  unchecked  real-time  radiosonde  transmissions  were 
used  to  obtain  c)9/dp,  (2)  temperatures  from  such  soundings  were  rounded 
to  the  nearest  whole  degree  Centigrade,  and  (3)  6p/6n   ,  at  the  mandatory 
level,  was  computed  from  numerically-processed  data  at  a  9*  surface  whose 
value  exceeded  that  at  850  mb  by  as  much  as  2.4C. 

Before  proceeding  to  an  extension  of  this  experiment  involving 
the  enhancement  of  baroclinicity  resolution  by  an  increase  in  data 
density,  it  is  well  to  remove  several  of  the  gross  approximations  in- 
herent in  the  computations  to  this  point.   For  this  purpose,  the  data 
for  850  mb,  0000  GMT  6  March  1966, were  further  manipulated.   Due  to  the 
basic  importance  of  G0  to  the  frontal  parameter,  DG9,  only  the  former  is 
in  focus  here.   Concentration  of  effort  is  on  the  dense-data  areas  of  the 
United  States,  Canada,  and  the  surrounding  ocean  areas. 

For  this  second  and  improved  approximation  to  G0  (hereafter 
symbolized  as  G0')  checked  radiosonde  data  were  used.   The  National 
Weather  Records  Center,  Asheville,  North  Carolina  (NWRC,  1966)  recomputes 
most  soundings  using  both  machine  and  manual  processing.   Among  other 
things,  the  radiosonde  data,  as  originally  recorded,  are  checked  for 
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proper  evaluation  and  inconsistencies  of  mandatory- leve 1  data  as  evidenced 
by  analyzed  charts;  also,  derived  height,  temperature,  moisture,  and  wind 
data  are  recomputed  for  50  mb  intervals  as  well  as  the  significant  levels. 
Such  checking  often  results  in  soundings  differing  considerably  from  the 
transmitted  version  in  both  significant  and  mandatory  temperature  infor- 
mation.  And,  of  course,  both  checked  and  transmitted  soundings  differ 
from  the  original  recording  and  any  systematic  smoothing  of  same  (Daniel- 
sen,  1959). 

In  view  of  the  increased  accuracy  of  the  checked  data  compared  to 
the  transmitted  soundings,  the  arbitrarily  imposed  accuracy  of  1C  in 
temperature  and  the  10  mb  criterion  for  significant  layers  were  removed. 
The  following  criteria  were  substituted  in  computing  the  baroclinicity 
parameter,  G©'.   To  compute  d9/dp  at  850  mb ,  temperature  data  were  re- 
corded to  the  nearest  . 1C  from  850  mb  and  the  first  significant  levels 
above  and  below  850  mb.   If  such  significant  levels  extended  more  than 
50  mb  beyond  the  850  mb  level,  then  the  800  mb  and  900  mb  level  data  were 
used  to  formulate  bQ/bp. 

Use  of  the  FNWC  processed -data  fields  was  maintained  for  computing 
the  slope  of  the  isentropic  surfaces  in  Eq .  (14).  However,  the  slopes  at 
each  station  were  obtained  from  a  more  exacting  differentiation  scheme 
made  available  by  Meteorology  International  (Holl,  1966).   As  noted  in 
Chapters  II  and  III,  the  derivatives  used  in  the  DG0*  program  were  com- 
puted from  the  Lagrangian  five-point  differentiation  formula  for  centered 
differences  (Eq.  [8]).  which  does  not  use  the  datum  at  the  grid  point  of 

The  elaborate  set  of  computerized  checking  procedures,  called  the 
Significant-Level  Punched  Card  Audit,  was  made  available  to  the  writer 
by  the  NWRC. 
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computation.   In  Meteorology  International's  program,  derivatives  are 
computed  from  a  16-point  4x4  array  of  grid-point  data  with  the  com- 
putation point  at  its  centrum.   The  program  is  a  scheme  which  minimizes 
smoothing  of  the  gradient  values  and  weights  heavily  the  influence  of  the 
immediate  surrounding  grid  points  as  compared  to  the  Lagrangian  com- 
putations employed  in  the  DG0*  model.   Further,  for  each  radiosonde 
station,  6p/6n  .  was  interpolated  for  the  850  mb  0*,  accurate  to  .1C. 

A  subjective  analysis  of  the  850  mb  G0*1  data  is  shown  in  Fig. 
22.   The  pattern  should  be  compared  with  those  of  G0'  and  G0*  in  Fig.  17. 
In  comparison  with  G0' ,  G01  shows  only  a  slight  (about  10  percent)  reduction 
in  the  subbaroclinic  area  but  the  stations  having  baroclinicity  magnitudes 
>  3C  (100  km)   have  increased  to  twelve  in  number.   These  figures  agree 
in  principle  with  a  sampling  of  data  from  those  radiosonde  stations  for 
which  there  exists  both  G0'  and  G0'  values;  namely,  63  percent  of  the  HBZ 
baroclinicities  have  increased  while  only  36  percent  of  the  SBZ  values 

have  decreased  from  G0'  to  G0' .   The  absolute  high  magnitude  of  G0'  is 

c  c 

21.30  C  (100  km)   at  station  879  in  southwest  Canada,  thus  increasing 
considerably  the  range  of  G0'  relative  to  G0' . 

Similar  to  G0'  versus  G0*  (Table  7),  G01  gives  more  focus  or 
resolution  than  G0*  (Table  8)  with  82  percent  of  the  HBZ  values  increasing 
and  89  percent  of  the  SBZ  values  decreasing  from  the  G0*  to  the  G01  com- 
putations.  It  is  to  be  noted  that  there  is  a  25  percent  increase  in 
comparative  data  points  in  Table  8  over  that  in  Table  7. 

There  are  some  important  differences  in  the  patterns  of  G0'  and 

G0'.   The  most  significant  of  these  occur  over  the  southwest  and  southeast 
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United  States,  the  former  associated  with  orographic  effects,  the  latter 

with  a  synoptic-scale  frontal  structure  whose  complexity  is  exemplified 
by  the  Canadian  analysis  in  Fig.  20.   Repositioning  of  the  axis  of  maxi- 
mum baroclinicity  in  eastern  and  central  Canada  is  also  apparent.   There 
can  be  little  question  from  Fig.  22  that  the  accuracy  and  resolution  of 
the  radiosonde  data  are  important  factors  in  locating  the  axes  of  baro- 
clinicity and  hence  the  hyperbaroclinic  zones. 


Table  8.   850  mb  AG©' (=  G0£  -  G0*),  0000  GMT  6  March  1966,  sum- 
marized for  147  radiosonde  stations  (North  America  and  adjacent  ocean 
areas).   HBZ  =  stations  where  G0<i  >  1C  (100  km)"1;  SBZ  =  stations  where 
G0'  <  1C  (100  km)"1.   Stations  with  absolute  instability  omitted. 


AGO1  Number  of  radiosonde  stations 
(C  [100  km]"1) 

HBZ  SBZ 

>  +  2.00  9  - 

+  1.00  to  +1.99  6 

+.60  to  +.99  9  0 

+.20  to  +.59  17  5 

-.19  to  +.19  16  34 

-.20  to  -.59  7  23 

-.60  to  -.99  1  6 

-1.00  to  -1.99  1  12 

<  -2.00  0  1 

In  general,  G01  differs  from  G01  by  virtue  of  greater  accuracy  and 

detail  in  computed  c*0/dp  and  a  reduction  of  smoothing  and  increase  in 
accuracy  of  6p/5n... 
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Many  problems  still  remain.   Nine  radiosonde  stations,  labeled 

"U"  on  Fig.  22,  indicate  absolute  instability  in  the  layer  within  50  mb 

of  850  mb.   Thus  6p/6n   ,  as  derived  from  FNWC ' s  mass  structure  model, 

is  in  error  at  these  stations.   A  potentially  greater  problem,  which 

encompasses  the  absolute  instability  difficulty,  is  the  validity  of 

6p/6n  .  in  combination  with  d9/dp. 
6 

Cross -section 
Analyses 

A  cross-section  analysis,  extending  from  station  78367  (Guan- 
tanamo  Bay,  Cuba)  to  72655  (St.  Cloud,  Minnesota)  is  used  as  another 
means  of  exemplifying  the  relation  of  baroclinicity  from  numerically- 
processed  and  checked  radiosonde  data.   Fig.  23  shows  the  former,  with 
9*  data  plotted  at  50  mb  intervals,  while  Fig.  24  shows  the  latter,  with 
9  computed  from  NWRC ' s  Northern  Hemisphere  data  tabulations  (1966)  using 
Eq .  (18).   In  the  case  of  station  78367,  significant- leve 1  data  were 
taken  from  the  real-time  transmitted  sounding,  as  checked  data  were  not 
available.   The  mixture  of  checked  and  unchecked  data  is  justified  for 
this  station  since  mandatory  level  data  on  the  two  versions  are  nearly 
identical,  a  fact  not  true  for  all  radiosonde  stations. 

The  subject  cross-sections,  oriented  approximately  in  the 
direction  of  V9*  at  850  mb ,  are  quite  interesting  for  several  reasons. 
The  sections  intersect  a  complex  frontal  structure  over  the  Florida 
area  (see  Figs.  20  and  22),  extend  from  a  near  maximum  to  a  near  minimum 
in  850  mb  temperatures  (Fig.  10),  and  pass  through  a  relatively  dense 
data  area. 

Contrary  to  the  expectation  that  9*  should  exceed  9,  nearly 
40  percent  of  the  area  from  the  surface  to  600  mb  shows  9  >  9*,  in  some 
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cases  (station  74794,  at  1000  mb)  by  as  much  as  5C .   The  situation  is 

similar  at  the  850  mb  level.   Generally,  9*  exceeds  0  in  the  relatively 

warm,  moist  areas  adjacent  to  the  frontal  zones,  while  in  other  areas  9 

is  greater  than  9*.   To  assist  in  this  analysis  the  quantity  9*  -  9  is 

plotted  at  50  mb  intervals  in  Fig.  23. 

Although  the  major  axis  of  warm  and  cold  air  have  similar  locations 

on  the  two  versions,  detail  on  stability  (from  station  78367  at  lower  levels 

to  72327  at  650  mb)  and  positions  of  the  minor  thermal  axes  differ  markedly 

in  pattern  (as  900-650  mb ,  from  78367  to  72206).   It  appears  that  isen- 

tropic  slopes  on  the  numerical  version  are  excessive  in  stable  areas 

(which  would  lead  to  excessive  G0'  values  in  HBZs)  and  deficient  in  or 

c 

near  actual  unstable  areas  (leading  to  deficient  G0 '  values  in  subbaro- 
clinic  and  unstable  layers),  relative  to  slopes  derived  from  the  checked 
radiosonde  data.   The  former  is  to  be  noted  in  the  vicinity  of  the  sloping 
stable  layer(s)  equatorward  of  station  72311  while  the  latter  is  obvious 
at  the  near-surface  levels  throughout  the  cross  section.   The  situation 
at  low  levels  (notably  at  72327  and  72655)  also  serves  as  an  example  of 
the  inability  of  the  numerical-data-processing  scheme  to  specify  the 
location  of  absolutely  unstable  layers  in  the  actual  atmosphere. 

Broadening  effects  due  to  the  numerical  processing  of  data  are 
quite  obvious  in  Fig.  23  within  and  either  side  of  the  major  zones  of 
hyperstability  and  hyperbaroc linicity ,  especially  from  850  mb  to  700  mb 
between  78063  and  72311.   However,  the  reverse  is  somewhat  true  geo- 
metrically below  900  mb  from  just  south  of  78063  to  74794,  where  9* 
exceeds  9  (moisture  effect)  while  9*  is  generally  less  than  9  north  of 
78063.   This  is  a  propitious  circumstance  for  locating  the  HBZ  in  the 
numerical  field  of  data. 
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It  is  of  interest  to  seek  further  quantitative  comparison  of  the 

isentropic  slope  values  since  the  mass -structure  version  of  this  factor, 

6p/6n_.,  is  retained  in  computing  G©1  .   Table  9  shows  the  calculated 
v  c 

850  mb  slopes  at  each  station  as  derived  from  the  two  cross-sections  as  well 
as  those  6p/6n_.  generated  from  the  numerically-processed  data  fields 
and  used  to  compute  G01  in  Fig.  22. 

Measurements  of  5p/6n  and  Baroclinicity--850  mb  point  slopes  of 
isentropic  surfaces  from  the  cross  sections  of  analyzed  processed  data 
(Fig.  23)  and  checked  radiosonde  data  (Fig.  24),  and  average  slopes  used 
in  computing  G©'  (Fig.  17b)  and  G©'  (Fig.  22),  are  compared  for  the  six 
stations  from  78063  to  72532  in  Table  9.   The  range  of  isentropic  slopes 
is  greatest  from  the  point  slope  computations  made  from  the  cross- 
sections.   Beyond  that,  clear  cut  differences  do  not  exist  between  the 
various  values  of  6p/5n. 

The  6p/6n  values  may  be  combined  with  dG/dp  computed  from  the 
checked  data  to  produce  the  baroclinicities ,  G0,  shown  in  Table  9.   The 
pattern  of  G0s  associated  with  the  isentropic  point  slope  measurements 
from  the  cross-section  in  Fig.  24  agree  with  the  complex  doublet  of  HBZs 
between  stations  72206  and  78367.   Such  is  not  indicated  in  the  other 
versions  of  G©  in  Table  9.   Somewhat  disturbing  is  the  relative  maximum 
of  G©  in  the  vicinity  of  72327  for  all  four  G0s .   Since  this  feature  is 
largely  due  to  the  steep  isentropic  slope,  the  G0  version  which  reduces 
its  significance  is  most  desirable.   Overall,  as  expected,  real  data  vice 
processed  data  and  point-slope  vice  average-slope  calculations  yield  a 
G0  picture  closest  to  the  mesoscale  nature  of  HBZs. 
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Table  9.   850  mb  6p/6n  (mb[100  km]") 0000  GMT  6  March  1966: 
Computations  of  point  and  average  values  for  stations  78063  to  72532 
in  cross -sections  of  Figs.  23  and  24.   G9*s(C[100  km]  ;  computed  from 
tabular  values  of  6p/6n  and  d9/dp  computed  from  checked  radiosonde  data. 


Station 

Point   slope 
from  cross-section 

Checked  data        Processed  data 
(Fig.    24)               (Fig.    23) 

Average   sic 
compu 

ipe   numerically 
ited    from 

Met.    Intl. 
Program 

Lagrangian 

5-point 

formula 

6p/6n 

78063 

48.5 

28.5 

41.0 

28.0 

74794 

21.0 

44.0 

49.0 

-- 

72206 

27.0 

62.0 

45.0 

35.0 

72311 

104.0 

90.0 

75.0 

52.0 

72327 

64.0 

150.0 

60.0 

72.0 

72532 

36.0 

25.5 

36.0 

20.0 

G9 


78063 
74794 
72206 
72311 
72327 
72532 


3.1 
1.0 
1.4 
0.9 
1.0 
0.4 


1.8 
2.2 
3.2 
0.8 
2.2 
0.3 


2.6 

1.8 

2.4 

-- 

2.3 

1.8 

0.6 

0.4 

0.9 

1.1 

0.4 

0.2 

The  resolution  of  baroclinicity  is  advanced  one  more  step  in  the 
next  section  where  the  detail  and  accuracy  of  the  analysis  is  increased 
by  a  scheme  of  extrapolating  significant-layer  G0  data  to  the  850  mb  level 
at  points  between  the  locations  of  radiosonde  stations. 
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G9  Extrapolations 

The  increased  resolution  of  the  field  of  G0',  as  compared  to  G0' , 
is  due  to  the  computation  scheme  and  nature  of  the  data;  the  data  density 
is  essentially  the  same  for  both.   The  average  distance  between  radiosonde 
stations  is  close  to  5  deg  lat  over  dense  data  areas  and  much  larger  than 
this  over  ocean  areas,  thus  limiting  the  enhancement  of  detail  in  the 
analysis  of  baroclinicity  for  any  G0  computation  scheme.   A  partial 
solution  to  the  data  problem  is  made  possible  by  an  extrapolation  technique 
which  has  its  basis  in  the  thermodynamics  and  synoptics  of  the  vertical 
structure  of  HBZs .   It  is  well  recognized  that  a  hyperbaroclinic  zone  of 
significance  is  generally  associated  with  a  relatively  stable  layer  of 
air  in  which  d0/dp  is  negative  and  large.   Vertical  cross-sections, 
oriented  approximately  along  V9  on  an  isobaric  surface,  show  such  layers 
extending  over  atmospheric  depths  of  several  thousand  feet  and  distances 
of  hundreds  of  miles  with  boundaries  identified  by  an  isentrope  or  a  small 
range  of  isentropes  (Saucier,  1955).   That  stability  and,  consequently, 
baroclinicity  have  considerable  uniformity  in  a  particular  stable  layer 
may  be  seen  from  the  cross  section  previously  discussed  (Fig.  24). 

Based  on  this  general  background,  selected  layers  above  and  below 
a  mandatory  pressure  surface  have  been  examined  for  stability  (d9/dp)  and 
isentropic  slope  (6p/6n)  for  the  purpose  of  extrapolating  baroclinicity 
data  to  the  mandatory  level  at  locations  between  the  radiosonde  stations. 
In  this  way,  the  vertically-resolved  temperature  data  contribute  further 
to  the  detail  of  the  baroclinicity  analysis.   The  specifics  of  the  extrapo- 
lation procedure  finally  adopted  are  outlined  below.   Application  is  made 
to  the  850  mb  level,  0000  GMT  6  March  1966. 
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d0/dp- -Extrapolations  are  derived  from  radiosonde  data  between 
700  and  950  mb  only.   de/dp  is  computed  at  900,  800,  and  750  mb,  using 
a  quadratic  fit  of  checked  radiosonde  temperature  data.   For  each  of  the 
three  levels,  data  points  are  selected  for  the  level  and  50  mb  below  and 
above  unless  an  intervening  significant  level  is  available.   For  example, 
to  compute  d©/dp  at  900  mb ,  0s  at  950,  900,  and  850  mb  are  used  if  signi- 
ficant levels  are  absent  between  950  and  850  mb .   If  significant  levels 
existed  at  920  and  855  mb,  then  data  at  the  920,  900,  and  855  mb  levels 
are  used  in  the  c)0/c)p  computation.   Such  a  procedure  obtains  three 
extrapolations  for  each  radiosonde  station  (unless  station  pressure 
<  950  mb)  and  controls  the  layer  from  which  the  extrapolations  are 
derived  (i.e.,  from  approximately  2500  feet  below  to  5000  feet  above 
850  mb).   Due  to  the  complexity  of  localized  stable  layers  near  the 
ground,  the  use  of  data  in  the  low  troposphere  below  the  850  mb  level 
must  be  limited.   A  less  arbitrary  system  could  set  the  lower- level 
limit  at  50  mb  above  the  earth's  surface,  vice  a  fixed  950  mb ,  thereby 
handling  data  in  mountain  and  flat  areas  in  similar  fashion. 

6_p_/5n- -Similar  to  the  previous  application  of  Eq.  (14),  this  param- 
eter is  computed  as  6p/6n   at  900,  800,  and  750  mb. 

G9"--The  product  of  de/dp  and  5p/6n   is  symbolized  as  G0"  for 
the  cases  of  data  computed  at  P   850  mb.   The  computed  G0"  values  are 
extrapolated  to  850  mb  along  the  0*  surface  intersecting  the  point  of 
computation  at  pressure  level  P  (i.e.,  900,  800,  and  750  mb)  and  in  a 
direction  along  the  850  mb  V0*. 

It  is  to  be  noted  that  the  system  adopted  here  extrapolates 
baroclinicity  to  850  mb  rather  than  simply  extrapolating  stability, 
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cWdp,  along  the  isentropic  surface.   One  reason  is  that  the  accuracy 
of  G9  ,  using  a  6p/6n   computed  at  P   850  mb  for  the  known  position 

C  H  '*  <C 

of  the  radiosonde  station,  is  better  than  that  obtained  from  use  of 
6p/6n   at  an  estimated  position  of  9*  at  850  mb . 

A  considerable  amount  of  manual  manipulation  of  data  was  necessary 
in  the  extrapolation  phase  of  the  project  since  this  aspect  was  handled 
as  a  feasibility  study.   Although  d©/dp  was  numerically  computed  at 
P   850  mb  from  checked  radiosonde  data,  the  0*  appropriate  to  the 
6p/6n  .  computation  was  determined  manually  from  P,  9*  soundings 
reconstructed  from  numerical  grid  prints  of  P  (0*)  at  2 . 5K  intervals. 
Accuracy  of  interpolated  0*  is  considered  to  be  +  .3K.   For  values  of 
0*  at  900,  800,  and  750  mb,the  isentropic  slopes  were  interpolated. 

Analysis  of  Combined 

G©'  and  G0"  Data 
c c 

Fig.  25  represents  an  analysis  of  the  combined  baroclinicity 
extrapolations  from  900,  800,  and  750  mb  to  850  mb  and  the  G01  compu- 
tations at  850  mb ,  the  latter  analyzed  in  Fig.  22.   The  extrapolations 
are  plotted  to  the  nearest  .1  C  (100  km)  "  to  distinguish  them  from  the 
850  mb  data  plotted  to  the  nearest  .01  C  (100  km)   .   Since  the  atmos- 
phere is  constrained  to  be  statically  stable  in  the  numerically-processed 
data  (Chapter  III)  used  for  6p/6n   ,  the  900  (800  and  750)  mb  extrapo- 
lations  always  appear  on  the  cold-air  (warm-air)  side  of  each  radiosonde 
station.   For  a  few  stations,  as  72365,  72485,  and  72576,  extrapolation 
data  only  exist   since  station  pressures  are  less  than  850  mb .   The 
extrapolation  data  are  particularly  helpful  in  areas  where  the  algebraic 
sign  of  850  mb  c30/dp  is  positive  (absolute  instability  cases,  plotted  as  U)  , 

as  for  stations  72232  and  72374. 
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Fig.  25.   850  mb  G0'  and  GO"  extrapolations,  0000  GMT  6  March 
1966.   Analysis  interval:   geometric,  with  r  =  2;  units:  C  (100  km)~l. 
Axes:   dashed  lines.   850  mb  G9'  data  (radiosonde-station  circles 
shaded)  plotted  to  nearest  .01;  extrapolation  data,  G0",  (location 
circles  not  shaded)  plotted  to  nearest  .1.   U  indicates  d©/dp  is 
positive;  850  mb  and  extrapolated  data  for  each  station  connected 
by  light  dashed  line. 


103 


Slopes  of  the  isentropic  surfaces  (6p/6n  ^)  are  indicated  by  the 
distances  of  the  G0"  extrapolations  from  the  parent  radiosonde  station. 
These  distances  correlate  inversely  with  the  magnitude  of  G0" .   As  an 
example,  compare  extrapolations  for  78016  and  72206;  the  latter  (former) 
indicates  nearly  vertical  (horizontal)  isentropes  and  relatively  high 
(low)  baroclinicities .   Validity  of  a  baroclinicity  computation  decreases 
with  increasing  distance  from  the  radiosonde  station  (i.e.,  decreasing 
isentropic  slope). 

Extrapolations  were  not  plotted  for  stations  where  the  9*  surface 
at  P   850  mb  could  not  be  located  at  the  850  mb  level  or  where  V9*  at 
850  mb  changed  by  180  degrees  between  the  station  and  its  associated 
extrapolation  point.   For  example,  this  situation  occurred  at  72867 
(76255)  for  the  900  (750)  mb  extrapolation.   A  9*  of  268. OK  (309. 8K) 
did  not  exist  at  850  mb.   For  this  study  all  the  extrapolations  (except 
those  indicated  with  a  question  mark)  were  held  to  be  valid  and  used  in 
the  analysis.   The  resulting  mesoscale  detail  in  Fig.  25  may  be  questioned 
Perhaps  those  extrapolation  points  separated  from  the  parent  radiosonde 
station  by  a  baroclinicity  axis  should  be  given  minimal  weight  in  the 
analysis.   The  same  approach  may  be  appropriate  for  data  associated  with 
nearly  horizontal  isentropic  surfaces  and  isolated  points  of  excessive 
high  and  low  G0"  values.   Subjective  computations  and  use  of  the  hybrid 
mixture  of  9*  and  9  information  very  likely  inflates  the  baroclinicity 
magnitudes  in  the  relatively  stable  layers  but  deflates  these  values  in 
the  layers  with  instability,  as  noted  from  a  discussion  of  Figs.  23  and 
24. 
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In  general,  the  number  of  baroclinicity  axes  from  G9"  and  G9' 

data  combined  have  doubled  compared  to  G01  alone.   This  relation  between 

the  two  fields  is  much  like  that  occurring  in  a  comparison  of  flight 

and  radiosonde  data  for  the  purpose  of  locating  the  jet-stream  axis. 

Each  synoptic-scale  jet  axis  may  be  reconstructed  into  a  multiple  number 

of  splinter  axes  by  increased  resolution  of  the  observational  data. 

Although  the  analyses  over  the  eastern  United  States-western 
North  Atlantic  region  appears  complicated  it  is  similar  in  interpretation 
to  the  detailed  Canadian  CAO  frontal  analysis  shown  in  Fig.  20.   For 
instance,  the  CAO  analysis  shows  850  mb  arctic,  maritime,  and  polar  fronts 
distinguishable  over  the  Florida  Peninsula.   These  frontal  entities  are 
in  much  the  same  position  on  Fig.  25.   For  the  remainder  of  eastern 
North  America,  the  comparison  in  number  of  CAO  and  numerical  HBZs  is 
good,  but  the  positions  of  these  zones  differ. 

The  profiles  of  850  mb  baroclinicity  between  stations  78367  and 

72655  on  Figs.  24  and  25  are  quite  comparable,  except  near  78367.   A 

complex  system  of  hyperbaroclinicity  between  stations  78063  and  72206  is 

found  on  both  charts  as  is  the  axis  of  maximal  values  near  72655.   Station 

72532  presents  a  problem  in  correctly  extrapolating  data  to  the  850  mb 

level  since  an  isentropic  minimum  occurs  in  the  layer  800-700  mb  very 

close  to  this  station. 

Over  mountain  areas  the  analysis  from  G9'  and  G9"  combined  differs 

J  c       c 

sharply  from  that  of  G9'  or  G9*  alone.   In  fact,  differences  between  any 
two  of  the  manually  or  numerically  generated  baroclinicity  analyses  are 
maximized  over  such  areas.  A  satisfactory  solution  to  meaningful  patterns 
over  these  regions  has  yet  to  be  achieved. 
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Summary 
The  various  experiments  to  enhance  the  resolution  of  baro- 
clinicity  for  the  purpose  of  HBZ  analysis  suggest   the  following: 

1.  There  is  considerable  baroclinicity  detail  inherent  in  the 
vertically-resolved  signif icant-leve 1  radiosonde  data.   Such  mesoscale 
detail  cannot  be  retrieved  from  numerically-processed  mandatory- leve 1 
data  alone. 

2.  Since  significant-level  radiosonde  data  (transmitted  or 
checked)  are  of  a  meso-microscale  nature,  the  baroclinicity  values 
derived  from  them  may  not  be  compatible  with  the  horizontally-projected 
synoptic-scale  fields  of  other  meteorological  parameters.   However, 
numerically  tuning  the  baroclinicity  field  by  use  of  smoothers  and 
filters  can  solve  this  problem  to  any  desirable  degree. 

3.  The  numerical  methods  and  grids  in  current  operational  use 
for  processing  and  analyzing  meteorological  data  are  not  appropriate  to 
the  scale  of  baroclinicity  shown  in  Fig.  25.   For  this  reason,  the  develop- 
ment described  in  this  chapter  must  be  regarded  as  a  second-generation 
model  for  the  purpose  of  analyzing  hyperbaroclinicity . 

Next,  attention  is  returned  to  the  mode  of  HBZ  analysis  feasible 
for  present-day  field  operations,  which  is  described  as  the  DG0*  model 
in  Chapter  III.   Considerable  insight  on  synoptic-scale  structure  and 
behavior  of  fronts  may  be  gained  from  a  diagnosis  of  this  model  in 
varying  space  and  time  frames  and  in  relation  to  other  meteorological 
parameters.   Thus,  Chapter  V  concentrates  on  the  four-dimensional  con- 
tinuity of  HBZs ,  as  illustrated  by  case  studies  from  winter  and  summer 
periods.   Chapter  VI  relates  the  HBZ  analyses  to  weather  satellite 
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observations  of  cloud  and  vertical  motion;  and,  Chapter  VII  demonstrates 
the  properties  of  time-averaged  intensity  and  location  of  HBZs  for  a 
selected  30-day  period  in  both  winter  and  summer. 

Together,  the  next  three  chapters  give  the  synoptic  meteorologist 
a  four-dimensional  view  of  frontal  analysis  not  heretofore  obtainable 
from  the  standard  manual  analysis  of  surface  fronts  only. 


CHAPTER  V 

THE  DG0*  MODEL:   CASE  STUDIES  OF  THE  SPACE  AND  TIME 
CONTINUITY  OF  HYPERBAROCLINIC  ZONES 

A  Summer-Season  Case  Study 
The  summer-season  day  of  20  August  1965,  with  emphasis  on  the 
0000  GMT  synoptic  time,  was  selected  primarily  for  a  study  of  the  space 
continuity  of  the  tropospheric  DG©*  fields.   Analyses  at  1000,  850,  700, 
and  500  mb  were  considered  in  the  substudy. 

DG9*  Analyses  at  1000, 
850,  700,  and  500  mb 

The  DG0*  operational  model  developed  in  this  treatise  explicitly 
indicates  boundaries  and  implicitly  the  intensities  of  HBZs .   In  general, 
most  meteorologists  are  familiar  with  the  synoptic-scale  boundaries  and 
their  vertical-plane  slopes  from  an  air-mass  point  of  view.   For  instance, 
it  is  well  known  that,  in  the  low- troposphere ,  frontal  zones  slope  toward 
cold  air  with  increasing  height,  with  slopes  generally  steepest  at  the 
lowest  levels.   In  the  middle  troposphere,  these  zones  of  hyperbaroclinicity 
become  more  horizontal  and  diffuse.   However,  not  much  is  known  of  the 
intensity  of  HBZs,  since,  normally,  the  parameter  is  not  displayed  in  a 
quantitative  fashion  on  the  daily  operational  charts. 

It  is  the  purpose  of  this  case  study  to  diagnose  the  numerically- 
computed  boundaries  of  HBZs  as  generated  from  the  processed  virtual 
potential  temperatures  (0*)  and  to  relate  these  analyses  to  the  previous 
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discussion  concerning  the  effect  of  numerical  processing  on  the  location 

and  intensity  of  HBZs  (Chapter  II). 

Figs.  26-30  illustrate  the  DG0*  fields  and  related  phenomena  for 
0000  GMT.  20  August  1965.   Isolines  on  these  figures  represent  a  manual 
trace  of  selected  positive-  and  negative-valued  lines  of  the  numerically- 
produced  front-location  parameter. 

In  order  to  focus  attention  on  the  HBZs  with  spatial  significance 
and  continuity,  certain  criteria  involving  intensity  and  pattern  were 
established  before  locating  the  axes  in  the  DG0*  fields.   The  major 
criteria  are  as  follows.   (1)  Analysis  is  limited  to  the  area  poleward 
of  15N.   (2)  Closed  isolines  are  not  drawn  around  single  grid-point  values, 
whatever  the  magnitude  of  DG0*.   (3)  The  positive  and  negative  axes  are 
drawn  only  for  segments  along  which  both  associated  axes  are  well  defined; 
DG0*  must  exceed  10  units  at  one  or  more  grid  points  along  one  of  the  axes. 
(4)  Axes  are  drawn  only  where  DG0*  magnitudes  exceed  5  units.   However, 
such  axes  are  connected  through  a  DG0*  area  of  less  than  5  units  magnitude 
(but  greater  than  zero)  where  the  interval  of  separation  is  less  than  2 
mesh  lengths,  and  the  associated  positive  or  negative  axis  is  continuous 
using  the  5-units  criterion.   See  DG0*  analysis  over  northern  United 
States  on  Fig.  26.   (5)  Axes  associated  with  areas  of  circular  DG0*  isolines 
are  not  drawn.   See  Kamchatka,  northern  Russia,  and  the  southeast  China 
coastal  areas  in  Fig.  26.   Such  areas  generally  signify  cut-off  pools  of 
cold  air  or  complex  isotherm  patterns.   (6)  Unless  a  HBZ  is  associated 
with  at  least  one  grid-point  G0*  value  in  excess  of  1  C  (100  km)    the 
DG0*  positive  and  negative  axes  are  not  shown.   Such  is  the  case  for  many 
ocean  and  low- latitude  HBZs.   (7)  Furthermore,  axes  obeying  the  criteria 
in  (1)  to  (6)  are  also  required  to  be  approximately  parallel  although  the 
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transverse  distance  separating  the  HBZ  boundaries  is  of  no  particular 
limitation . 

Establishing  criteria  of  this  character  has  as  its  purpose  an 
objective  separation  of  the  significant  from  the  insignificant  synoptic- 
scale  HBZs .   This  selective  screening  may  be  also  accomplished  on  an 
operational  basis.   Frontal  zones  weaker  than  those  in  focus  here  are 
considered  as  noise;  the  complex  zones  omitted  by  application  of 
criteria  (1)  to  (7)  generally  involve  undesirable  fiction  due  to  the 
numerical  processing  in  combination  with  grid  size  and/or  scale  of  the 
thermal  feature. 

1000  mb  (Fig.  26) --The  HBZs  at  1000  mbs  show  less  horizontal 
continuity  and  have  more  complex  structure  than  their  free -atmosphere 
counterparts.   Features  existing  at  1000  mb  but  not  found  at  success- 
ively higher  levels,  as  over  mountainous  terrain  ,  are  of  little  real 
consequence.   In  general,  oceanic  positive  and  negative  axes  are  parti- 
cularly sporadic,  especially  so  over  the  Pacific  Ocean.   The  DG9* 
intensities  are  greatest  over  land.   The  pronounced  HBZ  over  western 
North  America,  a  persistent  summer  feature,  is  a  function  of  the 
differential  large  scale  heating  either  side  of  the  land-sea  interface, 
with  some  contribution  from  the  regional  trade  wind  inversion  (Neiburger 
et  al.,  1961).   Occlusions  over  land  and  ocean  areas  are  poorly  indicated 
as  expected  from  discussions  in  Chapters  II  and  III. 

850  mb  (Fig.  27) --Smoothing  and  reduction  of  segmented  HBZs  are 
evident  when  compared  to  1000  mb .   Intensity  of  baroclinicity  change  at 
the  cold  boundaries  of  the  HBZs  appears  less  than  at  warm  boundaries  for 
the  non-stationary  frontal  zones  over  east  continents-west  oceans,  as 
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Fig.  26.   DG9*,  1000  mb ,  0000  GMT  20  August  1965.   Interval: 
10  x  10-2  C  (100  km)"2,  starting  with  the  5  isoline.   +  (-)  DG6* 
isolines :   solid  (dashed)  lines.   Warm-air  (cold-air)  boundaries  of 
HBZs :   heavy  solid  (dashed)  lines. 
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Fig.    27.      DG0*,    850  mb ,    0000   GMT   20  August    1965.      Interval: 
10  x    10-2   C    (100   km)"2,    starting  with    the    5    isoline.      +    (-)    DG9* 
isolines :      solid    (dashed)    lines.      Warm-air     (cold-air)    boundaries    of 
HBZs :      heavy   solid    (dashed)    lines. 
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expected  from  the  analysis  in  Chapter  II.   The  quasi-stationary  pattern 
due  to  differential  heating  over  western  North  America  is  st.ll  evident 
at  850  mb. 

700  mb  (Fig.  28)--More  extensive  and  continuous  f ron  ;s  'prevail 
at  700  mb  compared  to  lower  levels.   Reorientation  of  baroc 1 inicity  on 
the  North  American  west  coast  indicates  the  surface-induced  thermal 
structure  does  not  exist  at  this  level.   HBZ  boundary  values  generally 
show  an  increase  from  lower  levels,  especially  over  ocean  areas. 

500  mb  (Fig.  29)--Operational  use  of  frontal  zones  at  levels 
geometrically  above  700  mb  is  practically  non-existent  due  to  the  con- 
centration on  low-level  air-mass  structure.   However,  500  mb  frontal 
phenomena  do  exist,  as  evident  from  Fig.  29.   This  level  shows 
latitudinally  extensive  and  intense  HBZs  which  take  on  the  appearance 
of  spiraling  hemispheric-scale  jet  axes  (Riehl,  Alaka,  Jordan,  and 
Renard,  1954).   Many  of  these  zones  do  not  have  apparent  counterparts  at 
the  lower  levels.   Previous  authors  have  noted  a  similar  relation  between 
the  synoptic-scale  HBZs  of  the  lower  and  upper  troposphere,  as  well  as  the 
relative  importance  of  the  latter  to  movement  and  development  of  tropos- 
pheric  weather  systems  (Sanders,  1967). 

Vertical  Continuity 
of  HBZ  Boundaries 

The  positions  of  the  warm-side  boundaries  of  the  HBZs  at  the  five 
levels,  1000,  850,  700,  500,  and  300  mb ,  have  been  included  here  to  empha- 
size the  status  of  vertical  continuity  in  the  operational  DG0*  model  (Fig.  30) 
The  frontal  contours  that  lie  over  the  eastern  United  States,  Canada,  the 
northern  Atlantic,  and  northwestern  Europe  areas  are  internally  consistent. 
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Fig.  28.   DG9*  700  mb ,  0000  GMT  20  August  1965.   Interval: 
10  x  10"2  C  (100  km)-2,  starting  with  the  5  isoline.   +  (-)  DG0* 
isolines  :   solid  (dashed)  lines.   Warm-air  (cold-air)  boundaries  of 
HBZs :   heavy  solid  (dashed)  lines. 
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Fig.  29.   DG9*,  500  mb ,  0000  GMT  20  August  1965.   Interval: 
10  x  10-2  C  (100  km)"2,  starting  with  the  5  isoline.   +  (-)  DG0* 
isolines  :   solid  (dashed)  lines.   Warm-air  (cold-air)  boundaries  of 
«R7o  •   I-iaavv  snlid  (dashed)  lines. 
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Fig.  30.   Vertical  continuity  of  positive  DG0*  axes,  0000  GMT 

20  August  1965.   1000  mb :  ;  850  mb :   ---;  700  mb :   •••;  500  mb : 

///:  300  mb:  xxx . 
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Thus,  slopes  appear  most  logical  for  those  frontal  zones  which  are  not 
tied  to  orographic  features  and  where  data  are  ample.   Removal  of  1000  mb 
fronts  would  result  in  further  improvement  in  vertical  consistency, 
particularly  over  Asia.   The  weak  low-level  fronts  over  oceans  have 
poorest  spatial  continuity,  as  manifested  by  HBZs  tilting  toward  warm 
air  with  increasing  height.   With  reference  to  Chapter  II,  improper 
slopes  may  be  expected  for  HBZs  which  intensify  with  increasing  height, 
since  the  relative  broadening  effect  increases  with  the  frontal-zone 
strength.   However,  distinctly  separate  upper-troposphere  fronts  may 
account  for  the  seemingly  poor  continuity  between  upper  and  lower 
troposphere  over  the  ocean  areas.   Moreover,  the  liklihood  of  tropo- 
pause  intersections  with  the  300  mb  surface  further  complicates  the 
continuity  shown  in  Fig.  30. 

Although  not  shown  here,  the  front- location  parameter  was 
computed  from  the  objectively-analyzed  field  of  unprocessed  observed 
dry-bulb  temperatures  at  1000,  850,  700,  and  500  mb ,  0000  GMT  20  August 
1965.   Maximum  differences  existed  between  DG0*  and  the  dry-bulb  version 
at  1000  mb.   Negligible  differences  existed  at  500  mb  and  700  mb .   It 
was  also  evident  that  the  smoothing  inherent  in  numerically-processed  0* 
somewhat  compensates  for  the  potentially  greater  range  of  0*  compared 
to  9  at  all  levels  . 

A  Winter-Season  Case  Study 
This  section  considers  vertical  consistency  and  time  continuity 
of  numerical  frontal  analysis  for  the  winter  period  1200  GMT  25  January 
to  0000  GMT  28  January  1967.   As  before,  1000,  850,  700,  and  500  mb 
levels  are  subjected  to  the  diagnosis.   In  this  case  the  area  is  limited 
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to  the  North  American  section  of  the  Northern  Hemisphere  equatorward 
of  55N  latitude,  in  order  that  conclusions  be  drawn  from  evidence  in 
dense-data  areas  only. 

HBZ  Boundaries  in 
Vertical  Cross  Section 

The  medium  of  presentation  used  here  is  largely  that  of  the 

vertical  cross  section  depicting  HBZ  boundaries  in  the  layer  1000  to 

500  mb ,  as  derived  from  mandatory  pressure-level  analyses  of  0*  and 

DG0*  information.   For  this  series,  consideration  was  given  only  to 

-2 

those  HBZs  marked  by  six  or  more  grid-point  values  of  DG0*  >  5  x  10 

C  (100  km)   ,  at  least  one  of  which  >  15  units. 

Fig.  31  shows  a  schematic  1000  mb  frontal  contour  as  an  aid  in 
identifying  the  locations  of  the  cross-section  groups  in  Figs.  32-34. 
Nineteen  vertical  cross  sections,  intersecting  areas  adjacent  to  develop- 
ing frontal-wave  cyclones  over  central  and  eastern  United  States,  were 
selected  for  analysis.   The  sections  are  categorized  into  groups;  each 
group  is  homogeneous  relative  to  the  position  of  the  apex  of  a  frontal- 
wave  cyclone  and  the  movement  of  the  surface  front.   Both  cold-  and 
warm-side  boundaries  of  the  HBZs  are  shown  in  each  case.   Boundaries 
tilting  to  the  left  (i.e.,  toward  cold  air)  with  increasing  height  have 
positive  or  normal  slope.   0*  values  at  HBZ  boundaries,  for  each  mandatory 
level  considered,  are  also  shown.   Some  of  the  quantitative  conclusions 
which  may  be  drawn  from  this  sample  of  nineteen  cross  sections  are  shown 
in  Table  10  and  discussed  in  the  following  paragraphs. 

HBZ  Width- -Sawyer  (1955),  interpreting  British  flight  data  taken 
at  levels  ranging  from  1,600  to  20,000  feet,  indicates  that  the  widths  of 
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frontal  zones  average  115  n  mi,  with  a  contrast  of  5C(i.e.,  2.5  C 
[deg  lat]   ),  closely  matching  the  N  profile  in  Fig.  3.   However,  if 
radiosonde  data  are  used  to  determine  the  air-mass  contrasts  for  the 
same  times  as  the  twenty  observation  flights  then  the  average  width 
becomes  525  n  mi,  with  a  contrast  of  8  C  (i.e.,. 9  C  [deg  lat]   ), 
similar  to  the  W  profile.   Expressed  as  a  range,  Sawyer's  frontal-zone 
widths  vary  from  50  to  225  n  mi.   More  recently,  Danielsen  has  shown 
analyses  of  frontal  zones  over  the  western  United  States,  derived  partly 
from  middle-  and  high-troposphere  flight  data  in  support  of  Project 
Springfield  (Danielsen,  1964).   Baroclinicity  magnitudes  computed  from 
Danielsen's  selectively-sampled  stable  zones  above  2,000  feet  fall  in 
the  range  12  C  (80  n  mi)"   (i.e.,  9  C  [deg  lat]"  )  to  1.5  C  (40  n  mi)" 
(i.e.,  2.2  C  [deg  lat]   ).   Such  flight-data  values  generally  exceed  those 
of  the  hypothetical  profiles  of  Fig.  3,  which  were  designed  to  be  typical 
of  baroclinicities  diagnosed  from  a  constant-pressure  analysis  of  radio- 
sonde data, as  indicated  below. 

The  numerically-analyzed  HBZ  widths  in  the  January  1967  cross 
sections  average  close  to  400  n  mi,  increasing  slightly  with  increasing 
height.   Average  HBZ  intensity  decreases  with  increasing  height  from 
3.5  C  (deg  lat)"   at  1000  mb  to  1.6 C  (deg  lat)"   at  500  mb ,  thus  ranging 
from  the  N  to  M  profiles  in  Fig.  3.   As  expected,  most  of  the  associated 
DG0*  values  at  both  the  warm-side  and  cold-side  boundaries  of  the  HBZs 
also  decrease  with  increasing  height,  specifically  in  77  percent  of  the 
cases  up  to  700  mb.   Such  vertical  gradients  of  baroclinicity  in  the  low 
troposphere  agree  with  the  generally-held  concept  of  frontal  zones  over 
continental  areas. 
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Table  10.   Statistics  for  the  numerically-analyzed  HBZs  in 
Figs.  32-34.   All  numbers  represent  averages. 


Frontal  Parameter 


width  of  HBZ,  An  (deg  lat) 


AG*  through  HBZ  (C) 
A9*/An  (C  [deg  lat]  *) 


-1, 


A0*  along  HBZ  boundary  (C)   K: 

W: 

A0*/AZ  along  HBZ   boundary 

(C  [1000  ft]"1)  K: 

W: 


Slope  of  HBZ  boundary 


K: 
W: 


1000 


-Level  (mb)  — 
850       700 


■> 


500 


6.2 

22 

3.5 

6.6       6. 
16        13 
2.4       1. 

8  7.6 
12 

9  1.6 

1000/850 

850/700 

700/500 

11.5 
5.0 

9.6 
6.8 

11.6 
10.0 

2.3 
1.0 

1.9 
1.4 

1.5 
1.3 

1000/850 

850/700 

1000/700 

1/55 
1/15 

1/175 
1/164 

1/115 
1/90 

> 


Thermal  Character  of  HBZ  Boundaries- -HBZ  boundaries  are  known  to 
show  an  associated  increase  of  potential  temperature  with  increasing 
height.   O'Connor  (United  States  Navy,  1952)  observed  rather  small 
changes,  as  298  to  302  K  along  the  wintertime  polar  front  from  sea  level 
to  500  mb,  respectively.   However,  the  Canadians  indicate  larger  changes 
are  likely  along  frontal  zones  (Godson,  1950),  as  296  to  312  K  in  winter 
and  305  to  323  K  in  summer  along  the  polar  front  from  850  to  500  mb , 
respectively.   The  6*  values  used  here  steadily  and  consistently  increase 
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in  magnitude  along  the  HBZ  boundaries  at  a  rate  which  averages  close 

to  1.3  C  (1000  ft)  "  from  850  to  500  mb ,  and  hence  are  quite  compatible 

with  the  Canadian  values.   The  less  conservative  nature  of  the  cold 

boundaries  is  borne  out  by  the  higher  average  rate  of  A9*/AZ  compared  to 

that  of  the  warm  boundaries  in  the  layer  1000  to  700  mb  (i.e.,  2.1  versus 

1.2  C  [1000  ft]"1). 

Frontal  Slope—The  frontal-zone  slopes  are  greater  at  low  levels 
and  for  the  warm-side  boundary  of  HBZs ,  as  expected  from  conventional 
notions  of  fronts.   However,  the  slopes  appear  more  erratic  than  those 
which  would  be  expected  from  manually-drawn  fronts;  negative  slopes  (i.e., 
HBZs  sloping  toward  warm  air  with  increasing  height)  occur  in  20  percent 
of  the  cases.   In  the  layer  below  the  700  mb  level,  such  behavior  is 
partly  due  to  the  poor  numerical  resolution  inherent  in  the  mass  structure 
model  for  the  layer  775  to  1000  mb  (Chapter  III)  and  the  predominance  of 
steeply-sloping  cold  fronts  in  the  sample  considered  here.   In  the  layer 
above  the  700  mb  level,  the  picture  may  be  further  complicated  by  fronts 
which  are  rooted  in  the  high  troposphere  and  have  position  and  intensity 
not  intimately  associated  with  the  low-level  phenomena  (Fig.  35).   A 
similar  feature  was  evident  in  the  DG0*  analyses  for  the  summer-season 
case  study. 

Time  Continuity—The  cross-sections  of  HBZ-boundary  profiles 
shown  in  Figs.  32-34  are  graphed  in  such  a  way  as  to  aid  in  analyzing 
time  continuity.   Table  11  presents  information  on  two  time-continuity 
parameters  derived  from  changes  of  the  warm  (w)-  and  cold  (k)-side 
boundaries  of  the  HBZs  in  the  layers  1000/850  mb  and  850/700  mb. 

1.   The  percentage  of  time  that  the  algebraic  sign  of  the  12-hour 
slope  change  is  the  same  for  the  w  and  k  boundaries  of  the  HBZ. 
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Fig.    35.      Schematic   profiles    of  HBZ  boundaries    for  coexisting 
low-   and   high-troposphere    frontal   zones    (Endlich   &  McLean,    1957). 
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Table    11.      Time -continuity  statistics  for  numerically-analyzed 
HBZs     in     Figs.    32-34.      Column  a--Algebraic    signs  of    12-hour   slope 
change    same    for  w-   and   k-side   boundaries    of  HBZ .      Column  b--Algebraic 
signs  of    12 -hour   change    in  HBZ  width   and   slope    of   positive    and   negative 
DG9*  axes    differ. 


Cross   section  group 


I 

II 

III 

IV 

V 

VI 


Ave 

rage 

fre< 

quency 

(pe 

rcent) 

a 

b 

+ 

axis 

-axis 

100 

67 

67 

50 

67 

100 

50 

33 

33 

90 

64 

53 

100 

70 

33 

100 

50 

50 

2.   The  percentage  of  time  that  the  algebraic  signs  of  the  12-hour 
change  in  HBZ  width  and  slope  of  each  of  the  DG9*  axes  differ. 

The  internal  consistency  of  the  numerical  product  is  shown  by  the 
fact  that  over  80  percent  of  the  time  slope  changes  in  the  positive  and 
associated  negative  axes  of  the  DG0*  field  are  similar  in  algebraic  sign. 
There  is  no  apparent  pattern  in  the  changes  of  slope  and  width  of  HBZs. 
In  other  words,  as  the  frontal-zone  boundaries  steepened  or  became  more 
horizontal  with  time  the  HBZ  widths  may  either  increase  or  decrease. 

The  relation  of  12-hour  changes  in  HBZ  width  and  magnitude  of 
DG0*  is  quite  conclusive  in  the  case  of  increasing  HBZ  widths,  since  73 
and  87  percent  of  the  DG0*  values  along  the  associated  positive  and  nega- 
tive axes,  respectively,  decreased  during  the  same  period.   However,  the 
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situation  for  decreasing  HBZ  width  is  more  nebulous  with  an  average  of 
only  46  percent  of  the  changes  in  |DG0*|  being  positive  for  both  axes. 
Thus,  the  weakening  HBZs  (using  zone  width  as  an  indicator)  correlate 
well  with  changes  in  the  magnitude  of  DG0*,  but  the  cases  of  intensifying 
HBZs  are  indefinite  in  this  respect.   This  situation  is  not  unexpected  as 
a  consequence  of  poor  resolution  in  the  numerically-processed  data 
(Chapter  II). 

The  final  figure  of  this  chapter  (Fig.  36)  depicts  a  two-day  time 
continuity  (0000  GMT  26  January  to  1200  GMT  27  January  1967)  of  850  mb 
numerically-analyzed  frontal  positions.   Only  those  frontal  phenomena 
over  the  ample -data  United  States  area  are  in  focus  here.   The  aim  is 
one  of  illustration,  not  statistical  comparison.   The  time  continuity  of 
the  HBZ  warm-side  boundaries  appears  logical  in  the  area  east  of  the  Rocky 
Mountains.   However,  behavior  over  the  mountainous  southwest  United  States 
is  erratic.   Further,  the  absence  of  occluded  structure  is  to  be  noted. 
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Fig.  36.   Time  continuity  of  850  mb  numerically-analyzed  fronts 
(i.e.,  warm-side  boundaries  of  HBZs):   0000  GMT  26  January  1967  (thin 
solid  lines);  1200  GMT  26  January  1967  (thin  dashed  lines);  0000  GMT 
27  January  1967  (heavy  solid  lines);  1200  GMT  27  January  1967  (heavy 
dashed  lines) . 


CHAPTER  VI 

THE  DG9*  MODEL: 

A  CASE  STUDY  OF  HYPERBAROCLINIC  ZONES 

IN  RELATION  TO  CLOUD  AND  VERTICAL  MOTION 

The  hyperbaroclinic  zones  of  major  interest  to  the  field 
meteorologist  are  those  migratory  transition  regions  of  the  middle 
latitudes  separating  cold  polar  from  warm  tropical  air.   Such  zones 
are  invariably  associated  with  synoptic-scale  cloud  and  vertical  motion 
patterns  which  have  a  common  signature  in  the  large-scale  sense  but  vary 
in  small-scale  details  mostly  due  to  geographic  influences.   The  clouds 
are  the  visible  manifestation  of  saturation  due  to  relative  moisture 
content  and  the  time-integrated  effect  of  vertical  motion.   The 
intensity  and  shape  of  the  vertical  motion  field  is  largely  a  function 
of  the  two  baroclinicity-dependent  parameters  which  may  be  described  as 
the  vertical  rate  of  change  of  vorticity  advection  and  the  horizontal 
Laplacian  of  thickness  advection  (Krishnamurti ,  et  al.,  1966;  Sanders, 
1963).   However,  a  one  to  one  correspondence  between  baroclinicity  and 
cloud  patterns  or  baroclinicity  and  vertical  motion  patterns  does  not 
exist  due  to  the  diverse  combinations  of  the  thermal  and  vorticity  forcing 
functions.   Nevertheless,  the  relation  of  cloud  to  hyperbaroclinicity 
has  been  modeled  for  use  in  synoptic  analysis  and  forecasting,  starting 
with  the  Norwegian  concept  of  polar  front- phenomena  in  the  1920s 
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(Petterssen  1956a, b)  and  more  recently  by  those  observant  meteorologists 
attempting  to  exploit  the  operational  uses  of  the  weather  satellite 
(Boucher  and  Newcomb,  1962  ;  Anderson  et   al.,  1966). 

Fig.  37  qualitatively  implies  a  cloud-HBZ  relation  widely  held 
by  the  current-day  field  meteorologist.   Though  certainly  less  ideal 
than  the  early-century  concepts,  the  model  must  be  described  as  a 
heavily-smoothed  version  of  the  real  atmosphere.   The  extensive  cloud 
masses  and  bands  made  up  of  cirriform,  stratiform,  and  cumuliform  clouds 
near  the  cloud-mass  vortex,  in  the  warm  sector,  and  especially  adjacent 
to  and  within  the  frontal  zone,  are  indicators  of  significant  tropos- 
pheric  vertical  ascent.   For  other  areas,  especially  the  closed-cellular 
and  cloud-free  regions,  the  vertical  motion  is  likely  to  be  negative 
over  at  least  part  of  the  troposphere. 

Until  the  existence  of  a  model  to  quantitatively  specify  the 
location  and  intensity  of  frontal  zones,  an  objective  comparison  between 
hyperbaroclinicity ,  vertical  motion,  and  cloud  has  not  been  possible. 
The  research  described  below  concerned  itself  with  such  a  substudy,  with 
emphasis  on  the  relation  of  satellite-observed  cloud  to  the  HBZ. 

Selection  of  Period  and  Area  of  Study 
An  eleven-day  winter  period,  7-17  December  1966,  was  chosen  to 
investigate  the  relationships  of  satellite-observed  clouds  and  FNWC's 
vertical  motion  fields   (Hughes,  1967)  to  the  HBZs  derived  from  the 


FNWC  calculates  vertical  motion  from  a  diagnostic  a)  equation  in 
the  form  V2(ocD)+(pfrn/R)(dWdp2)  =V2(g/f  J[D,T])-J(T,t1)  -  J(D,g/f  [V2T]), 
where  a  is  a  stability  parameter,  J  is  the  Jacobian,  D  is  the  isobaric- 
height  deviation  from  the  United  States  Standard  Atmosphere,  r\   is  absolute 
vorticity,  and  other  symbols  are  standard. 
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DGG*  model  using  the  9*  appropriate  to  the  layer  1000/700  mb.   The 

weather  satellite  information  used  in  the  study  consisted  of  NMC- 

analyzed  nephanalyses  of  the  AVCS-type  ESSA  III  observations.   In 

order  to  consider  a  significantly  large  area  while  at  the  same  time 

maintaining  reasonable  concurrency  of  cloud  and  HBZ  analyses,  only 

the  nephanalyses  within  four  hours  of  the  0000  GMT  and  1200  GMT  times 

were  utilized. 

Use  of  the  satellite-observation  periods  of  2000  to  0400  GMT 
and  0800  to  1600  GMT  determined  two  major  geographical  areas  of  study. 
The  time  interval  centered  at  0000  GMT  gave  coverage  over  most  of  the 
North  Pacific  Ocean  plus  a  portion  of  North  America  west  of  95W,  while 
the  second  interval  gave  coverage  of  the  eastern  North  Atlantic  Ocean 
east  of  55W  ,  Europe,  and  western  Asia  to  90E.   Since  the  seasonally- 
imposed  northern  limit  to  the  daylight  satellite  view  is  close  to  60N 
and  since  the  numerical  frontal  analysis  is  not  too  reliable  at  latitudes 
less  than  15N  due  to  sparse  data,  the  south  and  north  boundaries  of  the 
areas  of  study  were  chosen  accordingly. 

The  December  1966  weather  situation  was  characterized  by 
abnormally  strong  zonal  flow  over  much  of  the  Northern  Hemisphere 
with  mean  troughs  near  both  coasts  of  North  America  (Dickson,  1967). 
In  the  North  Pacific  area  a  series  of  frontal-wave  cyclones  moved 
through  the  latitude  band  40-55N,  recurving  northeast  into  the  Gulf  of 
Alaska.   In  the  North  Atlantic  region,  the  major  storm  track  was  con- 
fined to  polar  latitudes  north  of  60N  due  to  the  presence  of  a  relatively 
strong,  stationary  subtropical  high  in  the  middle  latitudes.   Because  of 
the  blocking  situation,  only  the  trailing  portions  of  fronts  crossed 
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into  the  North  Atlantic  area  of  study.   Over  the  Eurasian  continent, 
the  weather  in  the  early  part  of  the  period  was  dominated  by  a  high 
pressure  cell  over  Siberia  and  weak  frontal  lows  in  the  Mediterranean. 
In  the  last  half  of  the  period  the  weather  was  generated  mainly  by  low 
cells  in  the  North  Sea  and  over  the  northern  part  of  the  continent  with 
frontal  activity  across  Europe  extending  southward  to  the  Mediterranean. 

In  the  illustrations  and  cloud  statistics  which  follow  only 
cloud  amounts  are  considered.   That  is,  no  attempt  was  made,  on  the 
basis  of  the  nephanalyses ,  to  make  any  distinction  between  high,  middle, 
and  low  cloud  types.   An  exception  to  this  procedure  occurred  for  the 
cases  of  cloud  bands  or  masses  composed  entirely  of  cirriform  cloud. 
Such  clouds  were  considered  to  have  no  direct  relation  to  the  low- 
troposphere  HBZs  and  were  not  included  in  the  statistics.   With  reference 
to  standard  satellite  meteorology  terminology,  the  cloud  amount  covered, 
C,  indicates  >  80  percent  cloud  cover  while  mostly  covered,  MCO,  specifies 
50-80  percent  sky  coverage. 

Patterns  of  Vertical  Motion  and  Cloud  in  Relation 
to  the  Numerically-Analyzed  Fronts 

Some  notion  of  the  empirical  relations  between  the  numerically- 
analyzed  1000  mb  front,  satellite-observed  clouds,  and  850  mb  vertical 
motion  may  be  obtained  by  viewing  a  3-day  sequence  (15-17  December  1966) 
for  the  North  Pacific  area  (Figs.  38-40).   The  area  is  not  necessarily 
the  most  desirable  for  such  a  demonstration  due  to  sparsity  of  data,  but 
the  wide  expanse  of  uniform  surface  aids  in  generating  patterns  of  a 
more  universal  nature. 
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Cloud  banding  appears  to  be  well-related  to  the  1000  mb  fronts 

in  orientation.   However,  as  expected,  the  colocation  is  not  as  ideal 
as  that  of  the  model  in  Fig.  37.   The  low- latitude  fronts  are  associ- 
ated with  long-trailing  cloud  bands  extending  from  middle  latitudes 
into  the  subtropics.   Many  of  the  HBZ  boundaries  are  within  the  area 
covered  by  cloud  masses  which  include  bands  not  evident  to  the  cloud 
analyst.   It  is  to  be  noted  that  not  every  DG0*  axis  is  simply  related 

to  a  cloud  mass  or  band  nor  does  every  such  cloud  entity  correspond  to 

-  -2 

a   synoptic-scale   HBZ.      The   vast  majority   of   DG0*  maxima   >   15  x    10 

-2 
(100  km)  '  are  found  in  the  covered  and  mostly  covered  areas. 

The  relation  of  vertical  motion  to  HBZs  is  more  difficult  to 

analyze.   As  discussed  earlier,  a  given  pattern  of  vertical  motion  is 

related  more  clearly  to  cloud  than  to  the  zones  of  hyperbaroclinicity. 

2 
However,  major  areas  of  850  mb  vertical  motion   (i.e.  centers  >  1.0 

cm  sec   )  are  somewhat  related  to  HBZs,  with  60  percent  of  the  centers 
of  ascent  located  within  such  zones  and  most  of  the  remainder  in  the 
adjacent  warm-air  subbaroc linic  zone,  similar  to  models  suggested  by 
Sanders  (1963)  and  Todson  (1964).   The  orientations  of  low- latitude  verti- 
cal motion  patterns  do  not  appear  to  be  closely  associated  with  either  the 
cloud  bands  or  DG0*  axes.   A  part  of  the  problem  of  interpretation  here  is 
due  to  lack  of  radiosonde  data  in  the  area. 

In  summary,  the  illustrations  shown  here  do  not  suggest  an 
unqualified  set  of  usable  relations  between  HBZs,  vertical  motions, 


2  -1 

This  area  was  contained  within  the  1  cm  sec    isoline  for  500  mb 

vertical  motion  more  than  95  percent  of  the  time.   Patterns  at  the  two 

levels  were  very  similar;  only  the  magnitude  differed. 
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and  clouds.   However,  it  is  obvious  that  HBZs  and  areas  of  850  mb 
vertical  motion  equal  to  and  greater  than  +1.0  cm  sec   are  associ- 
ated with  major  cloud  bands  and  masses.   Further,  there  is  the 
suggestion  that  the  less  significant  (i.e.,  very  narrow)  frontal  cloud 
bands  are  found  adjacent  to  and  oriented  similarly  to  the  positive  DG0* 
axes;  in  these  regions  the  vertical  motion  is  generally  less  than  4-  1 
cm  sec   .   In  order  to  quantitatively  specify  the  relations  between 
cloud  and  HBZ,  the  next  section  presents  statistics  derived  from  the 
7-17  December  1966  period. 

Satellite-Observed  Cloud  Bands  Versus  Hyperbaroclinic  Zones 
Stratification  by  Latitude 

Statistics  on  cloud  bands,  stratified  by  latitude,  are  shown 
in  Table  12.   Average  DG9*  values  are  indicated  for  each  of  the  three 
latitude  bands  in  the  two  major  hemispheric  regions  sampled,  namely, 
the  North  Atlantic  east  to  western  Asia  and  the  North  Pacific  east  to 
western  North  America.   The  relatively  high  DG0*  intensities  along  the 
positive  axes  indicate   that  the  warm-side  boundaries  of  hyperbaroclinic 
zones  represent  closer  approximations  to  a  first-order  discontinuity 
than  the  cold-side  boundaries,  as  concluded  in  Chapter  II.  Also,  the 
average  magnitude  of  the  DG0*  values  increases  with  latitude,  an  in- 
direct measure  of  the  poleward  intensification  of  the  frontal  zones. 

The  statistics  on  clouds  versus  HBZs  demonstrate   that  the 
average  percentage  of  C  and  MCO  cloudiness  in  the  HBZ  is  larger  than 
that  in  the  SBZ  at  corresponding  latitudes,  considering  all  latitudes 
between  15N  and  60N.   Conditions  in  the  low  latitudes  of  the  Pacific 
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area  and   the   high    latitudes    of   the  Atlantic  area  are  most   responsible 

for   the   near    two   to   one    ratio   of   cloud    percentages    in   the   HBZ   versus 
the   SBZ.      Although    the    combined    percentage    of  C   and  MCO  cloud   cover 
inside    (especially)    and   outside    the   HBZ   markedly   increases   with    lati- 
tude, the  actual  area    covered   by  clouds    in  each    latitude   band    increases 
more    slowly   since    the   area   ratios    of    the    two    latitude   bands,    30-35N 
and   45-60N,    to    15-30N  are    .86   and    .65,    respectively. 


Table    12--Relation   of   satellite-observed   clouds    to   hyperbaro- 
clinic   and    s  ub  baroclinic   zones,    stratified   by    latitude,    for    period 
7-17   December    1966. 


Latitude             Average   DG0*        Percentage    of  Percentage  Percentage 
band    (deg)       (C[100   km]"2)         latitude   band  of  HBZ  with  of   SBZ  with 
within  HBZ  cloud   amount  cloud   amount 

Negative    Positive  C  MCO  C  MCO 

Axes  Axes 


Atlantic-Eurasian-African   Region 


15-30N 

.11 

-.08 

39 

13 

12 

8 

12 

30-45N 

.19 

-.09 

30 

41 

18 

19 

15 

45-60N 

.20 

-.19 

37 

65 

9 

30 

7 

Area-weighted 

average 

for 

15-60N 

37 

31 

13 

16 

12 

Pacific-North  American  Region 


15-30N 

.11 

-.05 

24 

30 

16 

11 

8 

30-45N 

.17 

-.13 

37 

51 

26 

25 

16 

45-60N 

.20 

-.15 

33 

63 

15 

38 

14 

Area -weighted 
average    for 
15-60N 

30 

41 

20 

19 

11 
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Stratification  by 
Geographical  Area 

Table  13  considers  each  type  of  cloud  band  (frontal  covered, 
non-frontal  covered,  and  non-frontal  mostly  covered)  in  each  of  the 
following  six  geographical  zones:  A,  North  Pacific;  B,  western  North 
America;  C,  eastern  North  Atlantic;  D,  Europe;  E,  western  Asia  (C,  D,  E 
poleward  of  30N) ;  and  F,  Atlantic-Africa-Asia,  15-30N.   Here,  the  word 
"frontal"  refers  to  the  satellite  cloud  entity  labeled  as  such  by  the 
ESSA  nephanalyst,  apparently  when  the  band  is  categorized  as  C  and  is 
thought  to  be  associated  with  a  conventional  front  as  depicted  on  the 
NMC  surface  chart. 

Pacific  area  A  indicates  that  the  average  angle  between  the  C 
bands  and  their  associated  DG9*  axes  is  larger  (29  deg  versus  15  deg) 
and  has  a  wider  range  (10-90  deg  versus  5-45  deg)  in  the  frontal  than 
in  the  non-frontal  cases.  As  discussed  in  Chapter  III,  the  operational 
numerical  frontal  analysis  model  does  not  adequately  depict  occlusions 
while  the  amplitude  and  orientation  of  frontal  cloud  bands  relates  well 
to  such  fronts;  this  is  a  possible  explanation  for  the  observed 
relation.   Corresponding  values  in  area  C  appear  to  disagree  with  this, 
since  the  average  angle  in  the  frontal-covered  category  is  small. 
However,  as  noted  earlier  in  this  chapter,  the  storm  tracks  and  their 
attendant  occlusions  were  well  within  the  area  studied  in  the  Pacific 
but  not  in  the  Atlantic.   Therefore  the  Atlantic  and  Pacific  statistics 
do  not  conflict. 

Further  examination  of  Table  13  discloses  that  in  the  frontal- 
covered  category  a  greater  portion  of  the  cloud  band  fell  into  the  HBZ 
than  in  the  SBZ  in  both  ocean  areas  A  and  C,  but  the  cloud  amount  in  the 
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Table  13 --Relation  of  satellite-observed  cloud  bands  to  hyper- 
baroclinic  and  s ubbaroclinic  zones  in  the  latitude  band  15-60N, 
stratified  by  geographical  area:   North  Pacific--A;  North  America, 
west  of  95W--B;  North  Atlantic,  poleward  of  30N,  east  of  55W--C; 
Europe,  poleward  of  30N--D;  Asia,  west  of  90E--E;  and  the  Atlantic- 
Africa-  Asia  area,  15-30N  between  55W  and  90E--F. 


Area   Type  cloud   Number   Average  Range  Average  DG0*  Percentage 

band  and    of       angle  of  (C[100  km]"^)  of  cloud 

coverage    bands    between  +  angles  along  bands  in 

sampled  DG9*  axis  (deg)  positive  axis  HBZ 
and  cloud 
band  (deg) 


Frontal 

C 

18 

29 

10   to   90 

.17 

55 

A 

Non- frontal 

C 

9 

15 

5   to  45 

.13 

35 

Non-frontal 

MCO 

8 

8 

0   to   10 

.09 

33 

Frontal 

C 

1 

90 

.20 

40 

B 

Non-frontal 
Non-frontal 

C 
MCO 

0 
0 

Frontal 

C 

8 

14 

5   to  35 

.25 

61 

C 

Non- frontal 

C 

10 

29 

0   to   90 

.16 

51 

Non-frontal 

MCO 

1 

10 

.07 

50 

Frontal 

C 

2 

12 

10   to   15 

.12 

78 

D 

Non-frontal 
Non-frontal 
Frontal 

C 

MCO 

C 

7 
0 
0 

13 

0  to  30 

.16 

63 

E 

Non- frontal 

C 

10 

28 

0   to  90 

.23 

55 

Non- frontal 

MCO 

2 

12 

10  to   15 

.15 

20 

Frontal 

C 

0 

F 

Non- frontal 

C 

0 

Non- frontal 

MCO 

12 

Indefinite 

.11 

50 
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SBZ  is  by  no  means  insignificant.   This  result  agrees  with  the  inter- 
pretation of  the  model  in  Fig.  37,  the  broadening  effect  on  the 
numerically-analyzed  fronts,  and  the  extensive  zone  width  due  to  using 
9*  rather  than  9*. 

In  the  non-frontal  C  and  MCO  classes  combined  the  statistics 
are  inconclusive.   In  areas  A  and  E  most  of  the  cloud  banding  occurs 
in  the  subbaroclinic  region  while  in  D  this  is  true  of  the  hyperbaro- 
clinic  zone;  in  C  and  F  the  bands  are  evenly  split  between  the  HBZs  and 
SBZs .   Non-frontal  cloud  bands,  especially  MCO,  are  frequently  associ- 
ated with  the  stratus  and  stratocumulus  attending  the  stable  regions  of 
anticyclones  and,  as  such,  large  amounts  in  the  SBZs  are  to  be  expected. 

The  satellite  cloud  bands  associated  with  major  positive  DG9* 

axes  were  also  studied.   An  axis  was  classified  as  major  if  at  least 

-  _2  -2 

one  DG9*  value  >  15  x  10   C  (100  km)  '  occurred  along  it.   The  results 

were  similar  to  those  presented  in  Table  13  and  are  not  shown  here. 


CHAPTER  VII 

THE  DG0*  MODEL:   HEMISPHERIC  FREQUENCY  AND  AVERAGE  INTENSITY 

OF  HYPERBAROCLINIC  ZONES 

The  DG0*  analyses  yield  quite  readily  to  an  objective  calcu- 
lation of  the  frequency  and  average  intensity  of  HBZs .   Two  30-day 
periods  were  selected  to  demonstrate  a  synthesis  of  objective  1000  mb 
frontal  analyses  for  the  Northern  Hemisphere;  one,  15  January  through 
14  February  1966,  exhibits  winter-time  conditions,  and  the  other,  the 
month  of  July,  1967,  shows  mean  summer-time  conditions.   For  both  winter 
and  summer  periods  the  frequencies  of  subjectively-analyzed  NMC  surface 
fronts  were  also  calculated.   This  was  done  over  a  reduced  area  com- 
prising the  western  Northern  Hemisphere,  as  appropriate  to  NMC's 
facsimilie  transmissions  of  the  0000  GMT  and  1200  GMT  surface  and 
1000/500  mb  thickness  fields. 

The  frequency  and  intensity  statistics  were  accumulated  manually 
from  the  0000  GMT  and  1200  GMT  numerical  DG0*  analyses  (1000/700  mb  model) 

by  tabulating  the  occurrence  of  and  average  intensity  along  those  positive 

—  -2 

axes  where  DG0*  >  .05  C  (100  km)    for  each  of  3844  grid  squares.   Each 

square  measures  381  km  on  a  side,  true  at  60N.   The  NMC  frontal  positions, 

not  allowing  an  intensity  calculation,  were  tabulated  in  a  similar  way  for 

frequency  only. 
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Winter-season  Sample:   15  January  through  14  February  1966 

Fig.  41  highlights  the  important  features  of  frequency  and  in- 
tensity of  numerically-analyzed  1000  mb  fronts  for  the  period  selected. 
For  frequency,  the  maximum  possible  occurrence  per  grid  square  is  62 
which  is  equivalent  to  the  number  of  0000  GMT  and  1200  GMT  synoptic  times 
involved  in  the  statistic.   Minimal  frequency  is  zero,  indicating  no 

occurrence  of  numerical  fronts  in  the  grid  zone  during  the  period.   Aver- 

_2 

age  intensity  values  have  a  lower  limit  of  .05  C  (100  km)  "  since  only 

those  positive  DG0*  axes  with  values  equal  to  or  greater  than  five  units 
were  counted.   The  upper  limit  of  intensity  is  only  a  practical  one  since 
most  map-time  intensities  are  less  than  60  units  of  DG9*.   In  general, 
analysis  at  latitudes  less  than  20  deg  lat  is  omitted  due  to  the  unreli- 
able nature  of  data  for  these  sparse-data  areas. 

As  background  to  the  interpretation  which  follows,  it  is  to  be 
noted  that  the  last  half  of  January,  1966, was  anomalously  cold  over  North 
America  while  the  first  half  of  February,  1966  (that  is,  at  least  from 
February  5  forward),  represented  a  trend  back  to  normal  winter-time  con- 
ditions (Ludlum,  1966).   This  suggests  that  over  the  ample-data  North 
American  region  frontal  activity  should  show  supernormal  occurrence  at 
relatively  low  latitudes  with  subnormal,  but  not  zero,  frequency  at  polar 
latitudes.   In  fact,  this  same  trend  was  evident  over  all  of  the  western 
North  Hemisphere  (Posey,  1966;  Green,  1966).   The  eastern  section  of  the 
hemisphere,  especially  Europe  and  western  Asia,  was  associated  with  rela- 
tively more  disorganization  and  meridional  flow  as  compared  to  the  western 
section,  particularly  in  January.   The  numerically-derived  frontal  averages 
for  these  months  show  essentially  the  same  features. 
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Numerical  Frequency 

Axes  of  maximum  frequency  are  multiple  with  mostly  zonal 
orientation,  except  along  the  continent-ocean  interface  where  overtones 
of  coastal  orientation  are  evident.   Frequencies  equal  to  or  in  excess 
of  30  percent  are  concentrated  over  land  areas.   Regional  maxima  appear 
on  the  Canadian  west  coast  near  55N,  over  the  extreme  North  Pacific  in 
the  Kamchatka  Peninsula  area,  east  of  Japan  at  35N,  over  the  western  and 
southern  Mexican  coastal  area,  and  along  the  Asiatic  front  in  the  foothills 
of  the  Himalayas. 

There  are  several  areas  of  zero  or  near-zero  frequency,  as  noted 
from  the  3-percent  lines  in  Fig.  41.   The  more  extensive  areas  of  minimal 
frequency,  as  in  various  parts  of  the  Arctic  latitudes,  may  be  regarded 
as  the  traditional  air-mass  source  regions.   However,  certain  localized 
areas  of  this  type  are  more  difficult  to  explain,  especially  in  view  of 
the  high  frequencies  in  the  surrounding  grid  zones.   Examples  may  be  found 
over  the  east  coast  of  the  United  States  near  35N,  southwest  United  States, 
southern  Finland,  and  ocean  areas  in  the  latitudinal  band  45-50N.   Most  of 
the  local  minima  over  land  have  the  common  feature  of  being  in  a  mountain- 
lee  area  while  such  ocean  areas  more  likely  represent  tenuous  boundaries 
between  arctic  air  excursions  toward  the  south  and  maritime- polar  air 
traverses  toward  the  north.   Minimal  zones  south  of  the  most  equatorward 
positive  DG9*  axes  are  not  indicated.   These  zones  are  generally  equator- 
ward  of  20N  and  display  erratic  patterns  due  to  sparse  data. 

Numerical  Intensity 

The  intensity  field  shows  that  the  intense  maritime  fronts  occur 
farther  poleward  than  their  counterparts  over  the  continental  regions. 
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Fig.  41.   Frequency  and  intensity  of  numerically-analyzed  1000 
mb  fronts  in  the  period  15  January  through  14  February  1966.   Frequency 
isolines  for  3$  (2  cases),  20$  f 12  cases),  30$  (19  cases),  and  45$  (28 
cases).   Frequency  axes:   20-29$,  heavy  dashed  lines;  >  30$,  heavy 
solid  lines.   Areas  with  frequencies  <  3$,  shaded.   Axes  of  average 
DG9*  intensity  (units,  10-2  C  [100  km]"2):   15-20  units,  light  dashed 
lines;  >  20  units,  light  solid  lines  with  crosses.   Intensity  areas  > 
30  units  enclosed  by  dashed  line.   Intensity  axes  omitted  through  areas 
of  frequency  <  3$. 
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In  general,  highest  intensities  are  found  over  land  with  absolute 

-2 

maxima  of  .49  and  .47  C  (100  km)    on  the  Alaskan  west  coast  near  60N 

and  55N}  respectively.   Such  figures  are  of  significance  considering 
that  a  total  of  11  and  23  cases,  respectively,  are  associated  with  these 
intensities  . 

The  large-scale  relation  between  intensity  and  frequency  is  one 
of  maximum  intensity  lagging  the  most  frequent  position  of  the  front. 
That  is,  the  former  is  toward  the  coldest  air  from  the  axis  of  greatest 
frequency. 

NMC  Frequency 

It  is  of  interest  to  compare  the  frequencies  of  NMC's  manually- 
analyzed  fronts  (Fig.  42)  to  the  objective  version  in  Fig.  41.   The  format 
in  both  figures  is  the  same  for  easy  comparison.   However,  an  intensity 
comparison  is  not  possible  in  the  case  of  the  NMC  fronts. 

The  major  differences  between  NMC  and  numerical  in  this  short- 
period  mean  are  similar  to  those  found  in  the  single-time  views.   For 
instance,  the  many  occlusions  in  the  eastern  ocean  areas  give  more 
meridional  orientation  to  the  NMC  compared  to  the  numerical  frequency 
axes  since  the  latter  does  not  portray  such  fronts  very  well  (Chapter 
III).   Another  extensive  area  of  difference  is  eastern  Canada  and  the 
North  Atlantic  Ocean,  particularly  Greenland  and  Hudson  Bay,  which  shows 
high  frequency  in  the  DG9*  field,  but  near  zero  occurrence  on  the  NMC 
chart.   Such  a  difference  is  to  be  expected  since  it  is  not  common  practice 
to  depict  orographic  and  other  semi -permanent  frontal  zones  on  NMC's  daily 
operational  weather  charts.   The  NMC  and  numerical  patterns  are  most  com- 
patible over  the  Canadian  Archipelago,  mountainous  sections  of  western 
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Canada,  and  the  United  States  and  the  North  Atlantic  equatorward  of 
40N. 

Summer-season  Sample:   July,  1967 
Data  for  the  summer-season  month,  July   1967  (Fig.  43) ,  are  pre- 
sented in  the  same  fashion  as  the  winter  period  discussed  above. 
Initially,  the  data  were  stratified  by  the  upper-air  observation  times 
0000  GMT  and  1200  GMT  in  order  to  analyze  diurnal  variations,  but  one 
month's  data  proved  to  be  inadequate  for  this  purpose. 

To  aid  in  interpreting  the  July   1967  frontal  statistics,  the 
following  comments  regarding  the  Northern  Hemisphere  weather  are  of 
interest.   The  subject  month  was  much  like  January   1966  in  that  North 
American  temperatures  were  well  below  average  east  of  the  Rockies. 
Record  lows  were  common  and  associated  with  polar-front  penetrations 
well  equatorward  of  normal,  especially  over  southeast  United  States 
(Ludlum,  1967).   In  the  European  and  central  and  west  Atlantic  and 
Pacific  areas,  low-troposphere  high  pressure  systems  were  well  developed 
with  anticyc Ionic  and  supernormal  zonal  flow  at  higher  latitudes  than 
usual  for  July.   Thus,  frontal  positions  may  be  expected  poleward  of 
normal  at  these  longitudes.   In  the  Asiatic  sectors  the  flow  patterns 
were  not  suggestive  of  any  marked  or  unusual  frontal  phenomena  (Dickson, 
1967b). 

Numerical  Frequency 

Compared  to  winter,  there  are  fewer  and  less  extensive  frequency 
axes  with  values  >  20  percent  but  relatively  more  points  with  frequency 
>  30  percent  along  the  axes.   Isolines  of  45  percent  and  60  percent  are 
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Fig.  43.   Frequency  and  intensity  of  numerically-analyzed  1000 
mb  fronts,  July  1967.   Frequency  isolines  for  3$  (2  cases),  20$  (12 
cases),  30$  (19  cases),  45$  (28  cases),  and  60$  (37  cases).   Frequency 
axes:   20-29$,  heavy  dashed  lines;  >  30$,  heavy  solid  lines.   Areas  with 
frequencies  <  3$  shaded.   Axes  of  average  DG9*  intensity  (units,  10~2 
C  [100  km]"2):   12.5-15  units,  light  dashed  lines;  >  15  units,  light 
solid  lines  with  crosses.   Intensity  areas  >  20  units  enclosed  by 
dashed  line.   Intensity  axes  omitted  through  areas  of  frequency  <  3$. 
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to  be  found  over  Asia,  Africa,  Mediterranean  Europe,  and  mountainous 
North  America.   An  absolute  maximum  frequency  of  70  percent  is  reached 
along  the  southern  periphery  of  the  Iranian  mountains  (near  30N) ,  just 
inland  from  the  Persian  Gulf  and  the  hot  Saudi-Arabian  desert  to  the 
southwest . 

Areas  with  near-zero  frequencies  are  of  greater  extent  in  summer 
than  in  winter,  both  in  the  Arctic  and  as  buffer  zones  between  the 
zonally-oriented  major  frontal  axes.   The  multiplicity  of  such  minimal 
areas  points  out  the  reduced  latitudinal  migration  of  fronts  in  the 
former  as  compared  to  the  latter  season. 

Except  for  northern  sections,  frontal  frequencies  of  any  con- 
sequence over  oceans  are  non-existent.   However,  the  relation  of 
latitudinal  positions  of  continental  and  maritime  fronts  is  similar 
to  winter.   The  circulation  patterns  prevailing  in  July  1967  undoubtedly 
emphasize  this  relation.   The  zonal  extent  of  fronts  is  interrupted  by 
the  coastal  orientation  and  orographic  influence  over  the  western  contin- 
ental areas  of  Europe  and  North  America.   The  extensive  Asiatic  front 
has  shifted  to  the  polar-side  desert  basins  of  the  Himalayas  at  most 
longitudes,  a  seasonal  shift  expected  from  the  nature  of  the  thermal 
patterns  in  the  two  seasons. 

The  frontal  ring  over  mountainous  western  North  America  may  be 
explained  by  a  number  of  features  on  the  periphery  of  a  July  maximum- 
temperature  area  normally  centered  near  37N,  113W  (U.  S.  Navy,  1966). 
One  such  feature  was  discussed  in  Chapter  V  (Fig.  26),  namely  the  dif- 
ferential heating  across  the  land-sea  interface.   Other  sections  of  the 
ring- like  front  are  due  to  the  polar  front,  usually  found  near  the 
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Canadian-United  States  boundary  in  the  summer  season,  the  sparse-data 

and  cooler  mountain  and  ocean  sectors  along  the  south-southwest  fringe 
of  the  desert  area,  and  a  combination  of  moisture  and  temperature  differ- 
ences between  land  and  water  areas  to  the  east  and  southeast  of  the  hot 
elevated  basin.   The  last -enumerated  phenomenon  generates  the  dew-point 
front,  occurring  often  through  the  Texas  area  in  summer  (U.S.  Navy,  1952). 
A  less  prominent  but  similar  situation  occurs  over  northwest  Africa. 

Numerical  Intensity 

As  expected,  the  intensity  field  is  more  complex  and  much  reduced 

in  magnitude  compared  to  winter.   The  maximum  average  intensity  of  .21  C 

-2 

(100  km)   ,  found  just  off  the  Atlantic  coast  of  African  Morocco,  is  less 

than  half  of  the  absolute  maximum  noted  in  the  winter  sample.   However, 
intensities  over  mountainous  western  North  America  closely  approach  20 
units.   Most  zonally-oriented  fronts  at  latitudes  less  than  40N  and  most 
ocean  fronts  are  not  associated  with  intensity  axes  of  significance,  using 
>  15  units  as  the  criterion.   The  lag  of  the  maximum  intensity  zones 
relative  to  maximum  frequency  zones,  evident  in  winter  over  the  east  con- 
tinent-west ocean  areas,  is  not  apparent  in  the  summer  picture;  rather, 
the  orientation  of  frequency  and  intensity  axes  tends  toward  coincidence 
in  the  mean. 

NMC  Frequency 

The  pattern  of  axes  on  the  NMC  analyses  are  not  entirely  according 
to  expectation  (Fig.  44).   The  multiplicity  of  major  axes  of  frontal 
frequency  not  only  exceed  NMC's  winter  picture  but  numerical  in  summer  as 
well.   The  relation  to  numerical  is  least  expected  as  a  result  of  vari- 
ations noted  in  the  daily  patterns  (Chapter  III).  An  area  in  point  is 
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eastern  North  America  and  extreme  western  Atlantic,  usually  one  of 
maximum  agreement  in  the  NMC  and  numerical  frontal  analyses.   Since 
the  manual  product  does  not  suffer  greatly  from  sparse-data  boundary 
effects  and  since  the  seasonal  stationary  fronts  are  not  depicted  in 
the  daily  subjective  analyses,  the  frontal  frequency  pattern  over  south- 
western North  America  is  simpler  than  numerical.   The  usually  trouble- 
some extreme  northeast  Pacific  area  yields  different  solutions  on  the 
subjective  and  objective  versions.   In  general,  central  and  western  ocean 
areas  and  the  Arctic  region  show  the  best  comparison  between  NMC  and 
numerical  patterns. 

Prior  to  the  development  of  the  frontal  analysis  model  described 
here,  a  simple,  objective,  quantitative  method  for  measuring  the  location 
and  intensity  of  fronts  did  not  exist.   Thus,  frontal  climatology  has 
remained  qualitative  and  descriptive  over  the  years.   However,  it  is 
obvious  from  the  results  of  the  prototype  experiment  of  this  chapter 
that  quantitative  frontal  climatology  is  now  feasible  and  that  such  a 
synthesis  may  be  expected  to  relate  well  to  geography,  season,  and  the 
daily  weather  sequences. 


CHAPTER  VIII 

RESUME"  AND  PROPOSALS  FOR  FURTHER  STUDY 

Summary  of  Results 

The  following  major  accomplishments  were  achieved  during  the  re- 
search study  reported  on  here. 

1.  The  most  important  result  of  the  investigation  is  the  develop- 
ment of  the  first  operationally  successful  numerical-objective  model  for 
the  synoptic-scale  analysis  of  the  location  and  intensity  of  hyperbaro- 
clinic  zones  (HBZs).   The  scheme  locates  the  constant  pressure  HBZ 
boundaries  (i.e.,  fronts)  along  the  axes  of  maximum  baroclinicity  gradi- 
ent,  simulating,  in  principle,  the  positioning  of  fronts  by  subjective 
analysts  who  adhere  to  the  transition-zone  concept  of  frontal  phenomena. 
The  model  is  easily  adapted  to  a  computer-oriented  national  weather  ser- 
vice operation,  as  evidenced  by  the  1000  mb  Northern  Hemisphere  numerical 
frontal  analyses  and  prognoses  being  produced  twice  daily  by  the  Fleet 
Numerical  Weather  Central,  Monterey,  California  (FNWC) . 

2.  A  four-dimensional  synoptic  and  statistical  diagnosis  of  HBZ 
structure  and  motion  is   made  possible  by  the  development  in  _1.   For 
this  purpose,  modeled  and  real-atmosphere  virtual  potential  temperature 
data  were  employed  to  illustrate  the  relation  of  the  merits  and  deficiencies 
of  the  analysis  model  to  the  numerically-processed  and  computed  data.   Gen- 
eral characteristics  of  the  HBZs  are  shown  to  be  universal,  although  a 
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particular  mode  of  numerical-data  processing  and  computation,  namely 
that  of  FNWC ,  was  used  throughout  the  study.   Grid  size,  frontal- 
zone  intensity  (especially  that  due  to  zone  width),  air-mass 
dimension  (i.e.,  width  of  the  adjacent  subbaroclinic  zone),  and  use  of 
layer  versus  mandatory- leve 1  temperatures  are  shown  to  be  factors  in 
yielding  computed  frontal  structure  which  differs  in  degree  from  the 
real  atmosphere  as  a  function  of  the  numerical  manipulation  of  data.   In 
general, the  HBZs  are  weakened  and  the  widths  broadened  by  numerical  pro- 
cessing, with  these  effects  somewhat  proportional  to  the  original  intensity 
in  the  transition  zone.   For  a  given  mesh  size,  use  of  mean-layer  vice 
single- level  temperatures  reduces  the  broadening  effect  on  the  HBZs.   In 
this  regard,  use  of  the  1000/700  mb  mean  virtual  potential  temperature 
field,  9*,  was  found  to  be  suitable  for  generating  1000  mb  (or  surface) 
fronts.   The  reduction  in  mesh  length,  which  must  be  closely  tied  to  data 
density,  appears  to  offer  the  most  potential  for  improvement  in  accuracy 
and  detail  of  numerical  frontal  analysis. 

3.  Several  case  studies  from  the  years  1964-67  in  extended  time 
and  space  frames  indicate   the  broader  utility  of  the  HBZ  as  compared  to 
the  conventional  air-mass  associated  front.   The  numerical- ana  lysis  model 
has  the  advantage  of  locating  fronts  at  any  level,  whether  moving  or  oro- 
graphically  fixed,  and  not  just  those  fronts  separating  air  masses  of 
recognizable  thermal  characteristics  and  dimensions.   However,  occluded 
structure  is  not  evident  from  use  of  the  model  without  empirical  modi- 
fication.  Further,  the  relation  of  HBZs  to  vertical  motion  and  cloud  is 
shown  to  be  complex  and  not  operationally  usable  for  locating  HBZs. 

4.  The  study  also  proposes  a  more  elegant  model  for  the  purpose 

of  analyzing  the  mesoscale  structure  of  HBZs.   The  high  resolution  inherent 
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in  vertical  soundings,  combined  with  a  vertical  extrapolation  procedure, 

increases  the  resolution  and  amount  of  baroclinicity  data  to  the  extent 
that  the  fine  structure  of  baroclinicity  is  clearly  evident  on  the  con- 
stant pressure  mapping.   The  capacity  of  the  major  national  weather 
centrals  is  not  sufficient  presently  to  take  full  advantage  of  such  a  model 
but  it  is  offered  as  a  prototype  for  future  operational  use. 

Proposed  Areas  for 
Further  Research 

The  following  areas  for  extended  research  on  HBZs  appears  pro- 
ductive at  this  time. 

1.   It  has  been  demonstrated  that  use  of  vertically-resolved 
thermal  data  allows  the  mesoscale  structure  of  HBZs  to  be  analyzed.   In 
reality,  the  synoptic-scale  HBZs  comprise  a  multiple  number  of  these  sub- 
zones  of  hyperbaroclinicity  (Kreitzberg,  1964).   In  this  study  such  res- 
olution was  introduced  at  the  baroclinicity  stage  of  the  HBZ  analysis. 
It  is  proposed  that  the  more  basic  field,  temperature,  be  resolved  by 
extrapolating  the  vertical  temperature  structure  into  the  horizontal  or 
constant  pressure  plane  and  then  applying  the  subject  model  for  the 
numerical  analysis  of  HBZs.   With  a  significant  increase  in  data  density, 
a  judicious  grid-size  reduction  can  be  fully  effective  in  increasing  the 
resolution  of  the  analyzed  HBZs.   For  a  high-resolution  model  of  this 
type,  it  is  presumed  that  considerable  smoothing  of  temperature  and/ or 
its  derivatives  will  be  necessary  in  order  to  locate  the  boundaries  of  the 
synoptic-scale  HBZs,  as  found  by  previous  researchers  working  with  detailed 
frontal  structure  (e.g.,  Beniura,  1962). 
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2.  The  operational  surface-front  analysis  model  proposed  here, 
utilizing  9*,  is  essentially  a  thermal-  as  opposed  to  a  weather -oriented 
scheme.   With  little  guidance  from  national  or  international  meteorologi- 
cal commissions,  the  latter  approach  is  becoming  popular  with  the  field 
meteorologist  as  a  convenient  way  of  forcing  apparent  agreement  with  the 
satellite-observed  cloud  masses  and  bands.   Perhaps  regarding  the  frontal 
zone  as  the  centrum  of  a  circulation  entity,  (Bleeker,  1961;  Eliassen, 
1962;  Todson,  1964;  and  Kirk,  1966)  would  allow  relating  vertical  motion, 
and  hence  cloud,  more  intimately  to  baroclinicity ,  thus  joining  the  thermal 
and  weather  concepts  of  the  frontal  zone.   With  the  advent  of  night  and  day 
coverage  by  a  multiple  number  of  weather  satellites,  empirical  modeling  of 
vertical  circulation,  weather, and  baroclinicity  may  now  be  feasible. 

3.  The  limited  substudy  concerned  with  the  long-period  (as  30 
days)  frequencies  and  average  intensities  of  HBZs  indicated  the  potential 
use  of  these  charts  for  discriminating  between  the  fronts  fixed  to  oro- 
graphy and  those  which  are  progressive.   For  operational  analysis  and  fore- 
casting it  may  be  desirable  to  eliminate  the  orographic-type  fronts  in 

the  daily  analysis  of  HBZs.   Some  preliminary  regional  work  contributing 
to  the  solution  of  this  problem  has  been  done  by  Reed  (1961),  who  deter- 
mined tables  of  average  seasonal  surface  temperatures  for  synoptic  map 
times  over  the  western  United  States.   Reed  suggests,  but  does  not  show, 
that  the  deviation  of  temperature  from  a  long-term  mean  be  used  to  locate 
fronts  in  the  complex  mountain  areas.   The  existence  of  monthly  mean 
radiosonde  temperatures,  recently  updated  by  the  United  States  Navy  (1966), 
makes  such  an  experiment  feasible.   It  is  proposed  that  the  HBZ  model 
developed  here  be  used  in  connection  with  the  temperature  anomalies. 
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4.  Again,  with  reference  to  the  substudy  mentioned  in  3  above, 
a  computation  of  the  average  frequency  and  intensity  of  numerically- 
analyzed  HBZs ,  as  determined  from  a  suitable  period  of  several  years, 
will  allow  the  establishment  of  a  much  needed  quantitative  climatology  of 
Northern  Hemisphere  fronts. 

5.  Operationally,  there  is  need  for  a  computer-plotted  front  in 
addition  to  the  isolines  of  the  front-location  parameter  since  there  is 
some  subjectivity  in  locating  the  implied  boundaries  of  the  HBZ  on 
meteorological  charts  of  small  scale.   Some  work  has  been  done  on  this 
problem,  using  a  directional  third  derivative  of  0*,  but  the  effort  was 
terminated  due  to  the  errors  inherent  in  fields  of  such  high-order 
derivatives . 

Final  Remarks 

Experience  in  applying  the  numerically-analyzed  fronts  emphasizes 
the  need  for  a  man-machine  mix  (Snellman,  1969)  due  to  the  current  gross 
grid  size  and  poor  horizontal  resolution  in  the  numerically-processed 
operational  data.   Using  the  DG0*  model,  the  computed  frontal  positions 
may  be  regarded  as  reasonably  good  first  approximations  to  the  actual 
locations.   Mesoscale  adjustment  of  the  objectively- located  fronts  as  a 
function  of  the  adjacent  and  updated  pressure,  wind,  and  temperature 
observations  is  still  deemed  to  be  necessary.   Such  an  approach  does 
require  plotted  data  as  well  as  a  knowledgable  meteorologist  but  the 
time  and  effort  involved  is  much  less  than  a  totally  manual  approach. 

It  is  of  interest  to  note  the  recent  theoretical  work  by  Williams 
(1967)  which  has  as  its  goal  the  dynamical  modeling  of  mesoscale  fronto- 
genesis  from  initial  large-scale  conditions  of  a  simple,  vertically-bounded, 
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atmosphere.   Williams'  developmental  model  yields  vertically- 
sloping  frontal  zones  whose  boundaries  may  be  located  in  terms  of 
the  characteristics  of  potential  temperature  and  wind  profiles,  quite 
analogous  to  the  operational  frontal-analysis  model  described  here. 


APPENDIX  A 

PRELIMINARY  EXPERIMENTS    IN  HYPERBAROCLINIC-ZONE  ANALYSIS 

Experiments   Using  Wet-bulb   Potential  Temperature 

The    initial    impetus    for   the   study   of   numerical  hyperbaroclinic 

zones    (HBZs)    came    from  some   excellent   statistical  work  on   the   air-mass 

distributions    of  wet-bulb    potential   temperature,    9    ,    published   by 

w 

Harley  (1962).   9  has  been  the  prime  parameter  in  the  daily  opera- 
w 

tional  frontal-contour  analyses  emanating  from  the  Central  Analysis 
Office,  Meteorological  Service  of  Canada  (CAO) ,  since  1953. 

Harley  presents,  among  other  things,  monthly  means,  standard 

deviations,  and  modal  values  of  9  in  three  air  masses  described  by  the 

w  J 

Canadians  as  maritime  tropical,  maritime  polar,  and  maritime  arctic.   His 

data  sample  is  provided  by  a  careful  analysis  of  a  full  year  (1957)  of 

radiosonde  observations  from  89  stations  in  the  area  25-90N  and  15W  to 

170E  across  North  America.   The  method  used  to  determine  an  air-mass  9 

w 

allows  his  statistics  to  be  interpreted  as  usable  climatological  indicators 
of  the  warm-air  edge  of  the  HBZs  separating  two  adjacent  but  dissimilar  air 
masses.   As  an  example,  Table  14  contains  the  monthly  means  and  standard 
deviations  of  frontal  9  s  for  the  polar  front,  which  separates  maritime 
tropical  from  maritime  polar  air.   The  9  values  are  taken  from  Harley, 
but  rounded  to  the  nearest  .1C. 

Some  preliminary  work,  attempting  to  assess  the  utility  of  9  as 
a  front- location  parameter,  is  described  in  the  following  subsection.   In 
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order  to  utilize  the  empirical  results  of  Harley  and  yet  limit  the  initial 
experiments  on  locating  HBZ  boundaries,  the  following  conditions  and  pro- 
cedures were  imposed:   (1)  0  and  V0  are  considered  as  the  prime  front- 

w       w 

location  parameters,  (2)  only  the  poleward  boundary  of  tropical  air  (i.e., 
polar  front)  is  located,  (3)  the  850  mb  level  is  a  suitable  single  test 
level,  and  (4)  since  objectivity  is  the  goal,  any  successful  subjective 
analysis  scheme  should  be  a  prototype  for  a  numerical-objective  approach. 


Table  14.   Monthly  mean  0  and  0  and  their  standard  deviations, 
for  the  polar  front,  as  analyzed  by  CAO. 


Month 

Jan 
Feb 
Mar 
Apr 
May 
Jun 
Jul 
Aug 
Sep 
Oct 
Nov 
Dec 


0 

Standard 

0 

Standard 

W 

(°C) 

deviation 
(°C) 

e 
(°C) 

deviation 
(°C) 

14.0 

1.1 

42.1 

3.1 

14.0 

1.0 

42.2 

3.0 

14.2 

1.0 

42.6 

2.9 

14.9 

1.0 

45.0 

3.0 

15.5 

1.3 

46.7 

4.0 

16.8 

1.3 

50.8 

4.4 

17.7 

1.0 

53.9 

3.7 

17.7 

1.1 

53.9 

3.8 

17.4 

0.9 

52.9 

3.0 

16.5 

1.0 

49.8 

3.3 

15.6 

1.0 

46.8 

3.3 

14.9 

1.1 

44.8 

3.3 
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Relations  of  0  and  V0  to  the  Polar  Front 
w       w 

Summer  Case  Study: 
Manual  Analysis 

The  National  Meteorological  Center  (NMC)  and  CAO  850  mb  polar- 
front  analyses  over  North  America,  for  35  synoptic  times  (0000  GMT  21-26 
May  1963,  0000  GMT  21  June  to  1200  GMT  26  June  1963,  and  0000  GMT  12  July 
to  1200  GMT  16  July  1963),  were  screened  for  their  relationship  to  the 

0   field.   The  0  value  at  each  radiosonde  station  was  obtained  graphi- 
w  w 

cally  from  the  United  States  Naval  Weather  Service's  standard  thermo- 
dynamic chart,  the  AROWAGRAM,  using  transmitted  radiosonde  temperature 
and  dew-point  data. 

Specifically,  answers  to  the  following  equations  were  being 
sought  from  analyses  of  the  850  mb  polar  front  on  the  NMC  and  CAO  charts: 

1.  What  is  the  relation  of  NMC  and  CAO  0  s  along  the  front? 

w 

2.  Are  values  of  0  relatively  uniform  along  either  of  the 

w 

fronts? 


3.  Are  frontal  values  of  0  conservative  with  time? 

w 

4.  In  view  of  answers  to  1,   2,   and  3_,  do   Har ley's  data  assist 
in  objectively  placing  an  air-mass  boundary? 

With  reference  to  these  questions,  analysis  of  the  May,  June,  and 
July  1963  data  indicated  the  following  conclusions.   The  numbered  para- 
graphs refer  to  the  foregoing  questions. 

1.   A  sampling  of  the  polar-front  0  s  at  longitude  intervals  of 

5  deg  over  North  America  generated  the  percentage  frequencies  and  averages 

shown  in  Table  15.   Similarity  of  NMC  and  CAO  frontal  0  s  does  exist  in 

w 

the  mean,  but  apparently  there  are  differences  in  0  values  along  the  two 
versions  of  the  polar  front  on  a  given  day. 
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Table    15.      Percentage   frequency  of  wet-bulb   potential  temperature, 
0    ,   along   the   850  mb   polar   fronts   analyzed  by  NMC  and  CAO. 


e 

w 

Percentage  Frequent 

<°C) 

21-26  May    1963 

21-26  June    1963 

NMC                    CAO 

NMC                     CAO 

12 

1                         3 

12-16  July   1963 
NMC  CAO 


12 

1 

3 

13 

14 

7 

5 

2 

15 

3 

11 

7 

16 

13 

7 

14 

10 

4 

17 

15 

28 

12 

12 

16 

22 

18 

36 

36 

22 

46 

30 

28 

19 

23 

26 

14 

9 

26 

26 

20 

3 

16 

9 

13 

7 

21 

3 

2 

2 

4 

2 

22 

3 

7 

11 

23 

4 

Average 

17.5 

17.5 

17.4 

17.4 

18.6 

18.8 

Number   of 
Cases 

39 

57 

46 

2.      Inspection   of   the   average    frontal  values   of  9     as   a   function 


of   longitude,   by  analysis   center,    shows    the   geographical  variation  of  0 


w 


along  the  front  (Table  16).   The  variation  is  greatest  in  the  case  of  NMC 


fronts.   Apparently  the  polar-front  values  of  0  are  not  uniform  on  the 
average,  as  indicated  in  Table  15  as  well. 
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3.  Some  measure  of  the  conservativeness  of  frontal  0  values 

w 

and  other  features  may  be  obtained  by  noting  a  characteristic  day  in 
the  series.   Fig.  45  for  0000  GMT  and  1200  GMT  21  June  1963  shows  the 
relation  of  NMC  and  CAO  polar  fronts  to  0  data  and  Harley's  statistical 
estimate  of  the  frontal  positions. 

4.  Fig.  45  shows  that  the  orientation  of  manually-analyzed  and 
statistical  frontal  positions  are  similar  but  not  coincident.   In  general, 
the  temporal  non-conservativeness  and  spatial  non-uniformity  of  9  sug- 
gests the  probable  failure  of  trying  to  objectively  locate  fronts  based 
on  a  single  climatological  value  of  9  .   But  the  preponderance  of  NMC  and 
Canadian  positions  within  one  standard  deviation  (a)  of  the  monthly  mean 

9  is  suggestive  of  a  statistical  approach  in  which  the  warm-air  boundary 
w 

of  the  HBZ  is  searched  for  within  the  area  bounded  by  9  +  k,  where 

y  w  — 

k  <  3  C. 


Winter  Case  Study: 
Numerical  Analysis 

The  1963  summer  tests  were  considered  at  least  encouraging  enough 
to  continue  experimenting  with  the  application  of  statistical  values  of  9 
for  frontal  analysis.   However,  since  the  aim  was  to  develop  an  objective 
system  operationally  applicable  to  the  Northern  Hemisphere,  use  of  numeri- 
cally-processed data  and  fields   became  necessary  and  computations  of  9 
were  converted  to  the  numerical  mode.   In  this  connection  an  analytical 
expression  for  9   (Berry  et  al.,  1945)  of  the  following  form  was  programmed 
for  solution  on  FNWC's  Control  Data  Corporation  1604  digital  computer. 


Data  were  processed  in  accordance  with  FNWC's  mass  structure  model  as 
described  in  Chapter  III. 
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A^LK       RdeslOOO^ 

9w  =  9e  eXp   "  \TTC-  +  1000  C    >  (25) 

w     p  p 


where 

,L     m 


9e   =  9  eXp  ^cHv*  (26) 


P     c 

m  is  mixing  ratio,  L  is  latent  heat  of  condensation,  subscripts  s  and  c 
indicate  saturation  and  condensation  levels,  respectively,  e  is  vapor 
pressure,  1000  is  1000  mb,  and  other  notations  are  standard.   Due  to  the 
limited  nature  of  this  experiment  the  various  complex  numerical  approxi- 
mations and  iteration  schemes  used  to  compute  9  will  not  be  discussed 

r     w 

here.   Such  are  on  file  at  FNWC. 

The  synoptic  time  1200  GMT  5  February  1964  was  chosen  for  experi- 
mentation in  objective  frontal  analysis  using  the  numerically-produced 
wet-bulb  potential  temperature  field.   In  this  case  the  field  is  actually 

wet-bulb  virtual  potential  temperature,  9*,  since  FNWC's  processed  tempera- 

w 

tures  approximate  virtual  temperatures.   Fig.  46  shows  a  portion  of  the 
850  mb  9*  field  over  North  America.   The  area  within  9  +  3q,  appropriate 
to  February  (Table  14),  is  regarded  as  the  most  likely  region  within  which 
the  polar  front  is  located.   Using  9*  as  a  suitable  approximation  to  9  , 

r  °   W  rr  W* 

the  first  approximation  to  the  position  of  the  polar  front  is  the  zone 

(called  search  zone)  between  the  9*  isolines  of  11. 0C  and  17. 0C.   Follow- 

w 

ing  the  definition  of  the  front  used  here  (Chapter  I),  the  line  of  maximum 

V9*  within  the  search  zone  became  the  second  approximation  to  the  polar 

front.   In  this  case,  the  orientation  of  the  axis  of  maximum  V9*  is  simi- 

w 

lar  to  the  subjectively-located  fronts,  especially  CAO,  but  locations  differ 
This  and  other  tests  of  a  like  nature,  although  suggestive  of  the  direction 
in  which  to  pursue  the  research,  did  bring  about  the  abandonment  of 
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wet-bulb  potential  temperature  as  a  front  locator  for  the  following 
reasons . 

1.  0*  was  being  considered  solely  because  of  the  existent 

w 

air-mass  statistics  of  Harley  (1962).   By  using  the  Smithsonian 

Meteorological  Tables  (List,  1963)  the  0  statistics  may  be  readily 

w 

converted  to  their  0  counterparts.   Table  14  includes  the  monthly 

mean  frontal  values,  0  ,  so  determined. 

e 

2.  0*  computations  are  difficult  and  costly  in  computer  time 

due  to  the  multiple  numbers  of  scans  necessary  to  approximate  0*  to  an 
r  w 

acceptable  degree  of  reliability.   Consequently,  the  next  testing  phase, 

still  following  the  lead  of  the  Canadian  statistical  work,  utilized 

numerically-computed  0*  vice  0*. 
J  r      e       w 

Experiments  Using  Equivalent  Potential  Temperature 

For  reasons  discussed  above,  further  tests  with  0*  were  aborted 

w 

in  favor  of  the  more  easily  computed  equivalent  virtual  potential  tempera- 
ture, 0*,  using  a  form  of  the  analytical  expression  in  Eq .  (26). 

One  of  the  front-locator  parameters  developed  is  of  the  form: 

F.  =  0*n  |V0*|  ,  (27) 

i    e .  ■   e1 

l 

where  i  refers  to  _p_olar,  maritime,  or  arctic  fronts;  |V0*|  represents 

the  magnitude  of  the  gradient  of  0*;  and  0*  is  a  normalized  value  of  0*, 

e       e  e 

defined  as 

e*n  =0*  -  I 0*  -  9*  I .  (28) 

e.     e .     e    e . 
11  l 

For  the  purpose  of  this  experiment,  9*  +  ka  were  regarded  as  equivalent 

ei 

to  9    "  ka.   The  location  of  the  front  is  taken  as  the  axis  of  maximal 
ei 


values  of  F. . 

l 


178 

For  the  purpose  of  applying  the  scheme  of  Eq.  (27)  as  a  third 

approximation  to  the  850  mb  polar  front  of  1200  GMT  5  February  1964 

(Fig.  46),  the  following  procedure  is  implemented.   First,  the  9*  lines 

are  redesignated  as  9*"  lines.   Thus  9*  of  33. 2C,  42. 2C,  and  51. 2C 

e  e 

become  9*  of  33. 2C,  42. 2C,  and  33. 2C,  respectively.   Next,  a  search 

zone  defined  by  9*  +  3a  is  scanned  for  the  axis  of  maximal  values  of 
e  — 

F. .   Such  a  procedure  forces  location  of  the  front  between  the  maximum 
l 

V9*  and  9*  however  always  within  the  search  zone.   To  assure  this  con- 
e      e 

dition  in  the  machine  computation,  the  F.  values  may  be  arbitrarily  set 

to  small  values  outside  the  search  zone.   In  those  cases  in  which  9*  does 

e 

not  lie  in  a  direction  toward  warm  air  from  the  maximum  V9*  the  F.  com- 

e      l 

putation  of  the  front  location  is  less  accurate  than  that  dependent  on 

the  maximum  V9*  alone, 
e 

Enhancement  factors --In  order  to  ensure  placement  of  the  front  to 

the  warm-air  side  of  the  axis  of  maximum  V9*  in  the  search  zone,  a 

'   e ' 

weighting  factor,  W,  of  dual  form  was  used  in  combination  with  F.  to  gen- 
erate F' . ,  the  final  front-location  parameter  considered  in  this  statistical 
synoptic  approach.   Thus, 

F'.  =  F.  W  =  9*n  |V9*|  W,  (29) 

l    l      e .  '   e 

and,  as  for  F.,  the  axis  of  maximal  values  of  F1.  locates  the  polar  front, 
l  i 

At  the  boundary  of  and  inside  the  search  zone  defined  by  9*+  3a 

9*    ^ 

w  =  wi =  ^€rr^ '  <30) 

9*-3a 
e 

where  x  is  a  constant  for  a  particular  front.   Thus,  W,  increases  from 
one  at  the  cold-side  boundary  to  (9  +3a)/ (9*-3a) at  the  warm-side  boundary 
of  the  search  zone,   x  values  up  to  nine  were  used  to  enhance  the  search- 
area  gradients  of  9* 

e. 
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Outside  the  search  zone, 

6*\  X 

G* 
e 

which  always  yields  a  fraction.  %   values  up  to  twenty  were  used  to  reduce 

the  gradients  outside  the  search  area  to  extremely  small  fractional  values 

As  an  example,  following  this  scheme  in  the  1200  GMT  5  February 

1964  case,  Table  17  shows  the  weight  factors  and  other  information  which 

result  for  the  specified  0*  lines.   For  cases  where  maximum  VG*  is  at  the 

e  e 

location  G*  +  3a,  F'.  does  not  change  the  front- location  estimate.   For 
e         i 

all  other  locations  of  a  VG*  maximum,  F1.  yields  a  more  consistent  front- 

e  i 

position  estimate  than  F.  or  maximum  VG*  alone.   Fig.  46  shows  the  fourth 
r  l  e 

approximation  to  the  polar  front  using  the  F'.  approach. 

Table  17.   Table  of  related  values  from  Eqs .  (27)  and  (29),  as 
applied  to  a  February  850  mb  polar  front. 


F' 
Magnitude  of      „     i 

W  ==  """~"~ 
weight  factor         F. 

20  0.31 

G  -3   33.2   33.2  W,  9  1.00 

e 


G* 
e 

(°C) 

G*n 
e 

(°C) 

24.2 

24.2 

Symbolic  wei 

ght 

factor 

W2 

Wl 

Wl 

Wl 

W2 

G      42.2   42.2  W\  9  1.31 

e 

0  +3   51.2   33.2  W,  1.69 

e 

60.2   24.2  W„  20  0.31 


Empirical  adjustment  to  the  statistical  approach  is  almost  endless 
Wind  shear  or  other  front-associated  parameters  may  be  introduced  as  com- 
ponents of  the  weighting  factor  in  Eq .  (29),  various  values  of  x  may  be 


180 


r-l 

03 

0> 

I— 1 

P> 

•  " 

en  cm 

03 

(-1 

CO 

'H    v— ' 

o 

CU 

QJ 

■H 

4-1 

c 

C       • 

U 

C 

•l-l 

o    cr 

CD 

•r4 

1—1 

•r-t    td 

E 

4J 

3 

•^ 

T3 

co     E 

C 

a 

CO 

E    o 

O 

u 

•r-l       1-1 

-o 

■^S 

4J 

X   4-i 

(D 

O 

o 

4J 

O 

TD 

j-i  -   a, 

3 

m 

CL  fa 

& 

Cl 

E 

CO 

co     E 

o 

***N 

3 

o 

4J 

& 

XJ    E 

03 

CO 

(-1     -^ 

n 

•H     X 

<f 

, — , 

+ 

i-c  « 

vD 

* 

1). 

4J       E 

ff>    CI    *(U 
,_!   w|CD     " 

>,*?       X 

^  o  b  x 

CO  CO 

2     •+!•• 

£>    »-n  cu 

Pn  cni|o>   o 

w  N 

-a  £*  x: 
H    c        o 

S     CO    X)     J-i 
O  O)     CO     CO 

/-^  T3 
O  UO  C 
O  CM  3 
CN    ^    O 

r-l  JD 


CO 


■u    en    a) 

co     CT1    C  ■)<  0)CO 

fa     O  CD 

N  [>    -C 


a-  x! 
Cu  en 


CD    -r-l 


_c         s-i    o 

O  *  £3    5-i 

CO  t>   4-1 

CU 


en 


*  ?cn 

CD    oi  u 

CU  T-l 

u 

Cu  4-1 

X  C 

01  o 


E 


o 

00 

14-1 
O 

c 
o 


f 


cj 


CO 

c 

CO 

tu 
-C 
4J 


u-i  r — 

O    CM 


x  cr 

co  w 

en  E 

•r-l  O 

r-l 

C  4-1 

O 


d 

CO 


vD    14     It 
O-      O      E 


E 
5-i 
O 
Pl,    u-i 


a 

aco 

fa     -r4 


M 


X     E 
X     O    -H 

O     U     X 
S-i     D.   CO 

a.  a,  e 

a.  co 

CO      CO  4-1 

-o    o 

C    en  O    en 

■rlrl  O  -r)    . 

en   t4  cu    X  4I 

3   4-4  en    co  ^ 


181 
employed  in  the  enhancement  scheme,  and  so  on.   A  definite  problem  is 
proper  hemispheric  combination  after  each  of  the  three  fronts,  polar, 
maritime,  and  arctic,  is  determined  separately.   Various  arithmetic 

schemes  of  combining  the  individual  F.  or  F'.  fields  were  tried.   The 

a  11 

results   were   not   satisfactory. 

Since   a  more   elegant,   as  well  as   a  more   utilitarian,   model  was 
developed    (Chapter    II),    somewhat    forced   by   the   experiences    just    recounted, 
the    illustrations    here  will  be    truncated    to   one   experimental   numerical 
plot   of   the   arctic    frontal   field    (Fig.    47),    using  a    form  of  the   enhance- 
ment  scheme   just   outlined    (Eq .    29).      The   axes    of  the   elongated   F1      isolines 
represent   the   arctic   front.      Note   that   this   approach  yields   hemispherically 
continuous    and  mostly   zonally-oriented    fronts.      Comparison  with   the 
Canadian  arctic    front   is   good,   as   might   be   expected    from  the   dependence 

of  F'      on  Harley's   values    of  ©*;    these   values   are   used   to  some   extent   by 
a  e 

the  Canadians    in   their   frontal  analysis   approach.      A  more   extensive 
analysis    of  Fig.    47   is   not  warranted   in  view  of   the   resume1    immediately 
following . 

Summary 
The  various    frontal-analysis   schemes   described    in   this  Appendix 
restrict    the    fronts    to   a   definite   number,    as   well  as    to   specific   zones 
determined  by  air -mass   statistics    (Harley,    1962).      That    is,    fronts   out- 
side  the   zone   of  0     +  3a  for  polar,   maritime,   and  arctic   fronts   are   not 
allowed.      This    is   a  very  artifical  restriction.      Thus,    final   placement    of 
a   front    is   not  dependent   solely   on   factors    involved   in   the   definition  of 
a   front   accepted  here    (see  Chapter   I)   but   rather  on  air-mass    labels.      A 
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Fig.  47.   Frontal  field,  F'   generated  for  use  in  locating  the 
850  mb  arctic  front  at  0000  GMT  5  March  1964.   The  search  zone  boundaries 
given  by  F'  =  1:   heavy  lines.   The  axis  of  F'  within  search  zone  locates 
the  estimated  position  of  the  arctic  front.   The  CAO  version  of  the  arctic 
front  is  shown  for  comparison:   xxx . 
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dynamic-synoptic   approach   to  the   problem  is   desirable.      Such   is    the 

model   finally  adopted,   as   described   in  Chapter   II. 

Further,    0*  and   0*  were   abandoned  as    front- location   parameters 
we 

in   favor   of  0*  due   to   the   relative    inaccuracy   of  numerical  moisture 
analyses.      It  was    found   that  motorboating  values   and/or  very   low  dew 
points   were    improperly  rejected    in    the    numerical   processing   procedures. 
This   has    the   effect   of  giving  excessive   0*  and   0*  values    in   the   relatively 
cold  dry  air,    thereby  reducing   the   0*  and   0*  gradients   and  causing   the 
misplacement   of  the    frontal   positions.      Additionally,   detail   inherent   in 
the  moisture    field  was   not   resolved  by  analysis    procedures    tuned   to   the 
gross   mesh   size   appropriate   to  current-day  numerical  analyses.      It  might 
be  added   that   some  hope   of   improved  and  detailed  moisture   analysis   has 
been  given  recently  through  at    least   partial   incorporation  of  satellite 
data   in   the   numerically-processed   fields.      However,    objective    frontal 
analysis   experiments   with   parameters   dependent   on  detailed  moisture 
information  are   considered   premature   even  to  this   date. 
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